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Abstract and Summary 



The theory of the strong interaction, Quantum Chromodynamics (QCD), may 
be studied using perturbation theory, the standard tool of quantum field theory, 
for high energies. Around 1 GcV, the scale of the hadronic world, however, the 
coupling constant is increased and perturbative methods fail. In this domain. 
Lattice QCD provides a systematic approach to the non-perturbative study of 
QCD. The theory is formulated on a discrete Euclidean space-time grid which 
non-perturbatively regularises QCD and allows for computer simulations of the 
theory using Monte Carlo methods. 

The discretisation of the continuum action can be done in many different 
ways. Concerning the pure gauge sector of QCD, besides the standard discreti- 
sation, the Wilson action, there are several so-called improved actions which 
reduce artifacts coming from the discretisation. The most radical scheme of im- 
proving lattice actions, based on Wilson's Renormalisation Group approach, has 
been suggested by Hasenfratz and Niedermayer, namely the creation of actions 
that are classically perfect, i.e. there are no lattice artifacts on the solutions of 
the lattice equations of motion. 

In Lattice QCD, the energy of a physical state is measured studying the decay 
of correlators of creation and annihilation operators having an overlap with the 
state under consideration. If the state is heavy, these correlators decay very 
fast in time and one has to make sure that the temporal lattice spacing is small 
enough such that the signal of the correlator can be accurately traced over a few 
time slices before it disappears in the statistical noise. At the same time, one 
has to pay attention to choose the physical volume of the lattice large enough 
such that there are no significant finite-size effects. Both these requirements 
together lead to lattices with a large number of lattice sites, which means that 
it is computationally very expensive to perform the simulations. The obvious 
way out of this dilemma is to use a smaller lattice spacing in temporal direction 
compared to the spatial directions, i.e. using anisotropic lattices. 

Anisotropic lattices have been widely used in the last few years. Studies 
comprising excited states of nucleons, heavy-quark bound states, heavy meson 
semi-leptonic decays, long range properties of the quark-antiquark potential 
as well as states composed purely of gluons (glueballs) or gluons and quarks 
(hybrids) have been performed, mainly using the standard Wilson discretisation 
of the anisotropic action or a mean-link and Symanzik improved anisotropic 
action. A classically perfect anisotropic action has still been absent. 

It is thus the goal of this work to fill this gap, presenting a way of constructing 
classically perfect anisotropic gauge actions, building up on a recent parametri- 
sation of the isotropic classically perfect action (FP action) that includes a rich 
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structure of operators as it bases on plaqucttes built from simple gauge links as 
well as from smeared links. The procedure leading to the anisotropic action is 
examined analytically on scalar fields as well as in the quadratic approximation 
for gauge fields. We then construct the ^ = 2 action valid on coarse configura- 
tions occurring typically in Monte Carlo simulations. Its properties are studied 
performing measurements of the torelon dispersion relation (which serves as a 
means of determining the renormalisation of the bare (input) anisotropy), of 
the potential between a static quark and a static antiquark, of the deconfining 
phase transition and finally determining the spectrum of low-lying glueballs in 
pure gauge theory. 

The feasibility of iterating the procedure, obtaining a classically perfect ac- 
tion with ^ = 4, is briefly checked. Furthermore, we examine properties of the 
newly created anisotropic actions (as well as of the underlying isotropic action) , 
such as autocorrelation times of different updates of the gauge configurations 
in the Monte Carlo simulation or the computational overhead of the classically 
perfect action compared to the widely-used mean-link and Symanzik improved 
action and to the standard Wilson gauge action. 

It turns out that the construction of anisotropic classically perfect gauge 
actions is feasible. The iteration of the process yielding higher anisotropics 
seems to work as well. Measuring the renormalised anisotropy using the torelon 
dispersion relation turns out to be stable and unambiguous and shows that the 
renormalisation of the anisotropy is moderate and under good control. The mea- 
surements of the static quark-antiquark potential indicate that the violations of 
rotational symmetry are small if the (spatial) lattice is not exceptionally coarse. 
This shows that the parametrisation describes accurately the full action, which 
is known to have good properties concerning rotational symmetry. The study 
of the glueball spectrum is facilitated a lot due to the anisotropic nature of the 
action, even for (rather small) ^ = 2. Reliable results, including continuum 
extrapolations, are obtained for glueball states having much larger mass than 
the highest-lying states that could have been resolved with the same amount 
of computational work using the isotropic action. However, the scaling prop- 
erties of the glueball states, i.e. the behaviour of the measured energies as the 
lattice spacing is increased, seems to be rather unfavourable compared to the 
mean-link and Symanzik improved anisotropic action. Especially, this is the 
case for the lowest-lying scalar glueball which could be caused by the presence 
of a critical end point of a line of phase transitions in the fundamental-adjoint 
coupling plane assumed to define the continuum limit of a scalar field theory. 
Our action includes in its rich structure operators transforming according to 
the adjoint representation. If their coupling (which we do not control apply- 
ing our method) lies in a certain region, the effect of the critical end-point on 
scalar quantities at certain lattice spacings (sometimes called the "scalar dip") 
may even be enhanced compared to other (more standard) discretisations with 
purely fundamental operators. Concerning the cut-off effects observed in the 
glueball simulations, we have to add that other effects, such as effects due to 
the finite size of the lattices, may also be present and that this issue requires fur- 
ther study. Furthermore, the computational overhead of the classically perfect 
anisotropic action compared to the standard Wilson action and as well com- 
pared to the mean-link and Symanzik improved action is considerable and one 
has to weigh up the pros and cons carefully, if one is about to choose a gauge 
action for large-scale simulations. 



Chapter 1 

Introduction 



This thesis covers work done in collaboration with Urs Wenger and Ferenc Nie- 
dermayer. Results already extensively discussed in Urs Wenger's PhD thesis 
that are necessary for the work presented here are briefly recapitulated. Parts 
of the results have already been published |^-|^ . 

The aim of this Introduction is to motivate the work that has been done as 
well as to present the basic knowledge necessary for the subsequent Chapters. 

1.1 Motivation 

The most systematic approach to the non-perturbative study of QCD, the theory 
of strong interactions, is Lattice QCD (LQCD). The QCD action is discretised 
by putting the theory on a space-time lattice. As QCD, LQCD needs as in- 
put the quark masses and an overall scale. Then any Green's function may be 
evaluated by taking averages of certain combinations of the lattice fields (mea- 
suring operators) on an ensemble of vacuum samples (lattice configurations). 
This allows the study of masses (spectroscopy) as well as the extraction of ma- 
trix elements. In contrast to the experiment, parameters such as quark masses, 
boundary conditions or sources may be arbitrarily varied in LQCD, which allows 
a wide range of studies to be contrasted with experimental results. 

In order to measure the energies of states in LQCD, one studies the decay 
in time of correlators of respective operators. Thus the distance between lattice 
sites has to be small enough such that the signal of the correlators may be traced 
along several time-slices of the lattice before it disappears in the statistical noise. 
On the other hand, however, the physical size of the whole lattice has to be 
large enough such that there are no significant finite-size effects. Both these 
requirements together lead to lattices with a large number of lattice sites that 
are computationally expensive. 

The obvious way out of this dilemma is the use of anisotropic lattices, whose 
spatial lattice spacing as , the distance between two neighbouring lattice sites in 
spatial direction, is larger than at, the one in temporal direction. The anisotropy 
or aspect ratio is conventionally denoted by ^ = Os/at- 

Anisotropic lattices have been widely used for studying heavy states in 
gauge theory for the last few years. For measurements in quenched QCD 
see e.g. ||^ for N* masses, |6|-H] for heavy quarkonia using relativistic heavy 



3 



4 



Chapter 1. Introduction 



quarks, for heavy quarkonia including heavy hybrid states composed of 

quarks and glue or [^,|l6| for heavy-light mesons employing NRQCD. Another 
work has studied heavy meson semi-leptonic decays using NRQCD. Stud- 
ies in pure SU{^) gauge theory have mainly examined glueballs p^[2^. The 
study of string-breaking using static quarks in the adjoint representation pH 
employed an anisotropic action to be able to have large spatial separations of 
the sources on a lattice with a moderate number of lattice sites. Other studies 
of static potentials including higher representations comprise works like |^^-^ . 
Anisotropic lattices have been also used to study (quenched) QCD at finite 
temperature ||26|-|30t] as the fine lattice in temporal direction allows for a fine 
variation of T = l/(atNt) without the need of having to increase the number 
of spatial lattice sites. A study of the autocorrelation of spatial and temporal 
S'?7(3) gauge operators on anisotropic lattices has been performed in Q. Per- 
turbative properties of improved anisotropic gauge and fermion actions, e.g. the 
dependence of the couplings gs and gt on ^ as well as the ratio between the 
renormalised and the bare anisotropy ?? = Cr/^Oj have been studied in |]3^-^. 
Matching parameters for the study of matrix elements (mainly for heavy meson 
semi-leptonic decays) have been calculated in |36|. 

The discretisation of the continuum theory leads to cut-off effects (or lattice 
artifacts) in the results obtained on the lattice. These effects grow as the lattice 
spacing a is increased (i.e. as the lattice gets more and more different from the 
continuum). As simulations using anisotropic actions are carried out on lattices 
where the spatial lattice is rather coarse — in order not to have too large cut- 
off effects it is advantageous to use improved discretisations of the continuum 
action suppressing such artifacts. Whereas it has been common to use improved 
anisotropic gauge actions for some years, above all in glueball simulations, the 
most radical concept of improving actions, the construction of classically perfect 
actions, has never been applied to this subject. It is thus the goal of this work 
to construct anisotropic classically perfect gauge actions]^ and to examine their 
properties. 

Different ways of constructing improved gauge actions are discussed in the 
remainder of this Chapter. On one hand, there are actions that are improved 
employing perturbation theory, on the other hand, the concept of renormali- 
sation group provides means of improving actions non-perturbatively. One of 
the most radical ways of improving actions is the construction of the classically 
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perfect Fixed Point action, which is described in detail in Sections 1.3.2 and 

In Chapter || we present some physical objects that may be studied in the 
context of pure lattice gauge theory and explain how we measure them in prac- 
tice. The dispersion relation of the torelon is used to determine the renormalised 
anisotropy of the parametrised action as well as to provide an estimate for the 
scale at a given coupling (3. The measurement of the critical temperature Tc 
corresponding to the deconfining phase transition of pure gauge theory yields as 
well information about the (temporal) scale of a simulation. To determine the 
lattice spacing precisely as well as to obtain information about the rotational 
invariance of the parametrised FP action, we employ the static quark-antiquark 
potential. Finally, we describe bound states of pure glue, the glueballs, and 
discuss how they are measured on the lattice. 



^In the following, we will often drop the word classically for brevity. However, it has to 
be understood that the perfect actions in this work are always classically and not quantum 
perfect. 



1.2. Perturbatively Improved Gauge Actions 
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Chapter contains a brief summary about the newly parametrised isotropic 
FP action The action is presented as well as simulations studying the 

deconfining phase transition, the static quark-antiquark potential and glueballs. 
The scaling properties of the isotropic action are discussed as well. 

In Chapter ^ we present the method used in this work to construct aniso- 
tropic classically perfect actions. Instead of repeating the whole construction 
programme as described at the end of this Chapter, we use the parametrised 
isotropic FP action and perform one single, purely spatial, blocking step to end 
up with a perfect anisotropic action with anisotropy = 2. Furthermore, we 
describe how the isotropic parametrisation using smeared links may be gener- 
alised such that it may as well account for the anisotropy between spatial and 
temporal directions. 

This approach is then tested on free scalar fields that may be treated analyt- 
ically as well as in the quadratic approximation of the FP action, as described 
in Chapter |^. It turns out that in both cases the ansatz works nicely. 

The construction of the ^ = 2 perfect action in full SU(3) gauge theory 
is presented in the next Chapter, including as well results of measurements of 
torelon energies , the static quark-antiquark potential, the critical temperature 
and the glueball spectrum. 

Chapter describes how the spatial blocking step may be repeated to obtain 
actions with larger anisotropics and presents first results about the renormalised 
anisotropy of the ^ = 4 perfect action. It turns out that the iteration of the 
spatial blocking step seems to be feasible and poses no additional difficulties. 

In Chapter || we present results of measurements of autocorrelation times as 
well as information about the overhead to be expected in simulations using the 
perfect anisotropic actions compared to the standard Wilson action. 

Conclusions and prospects are given in Chapter |^. 



1.2 Perturbatively Improved Gauge Actions 

On very smooth lattice configurations, i.e. configurations with a very small lat- 
tice spacing a close to the continuum, any discretisation of the continuum gauge 
action yields basically the same results, the leading lattice artifacts of O(a^) are 
very small. Increasing the lattice spacing a, the form of the discretisation gets 
more and more important and simple discretisations of the continuum action 
(e.g. the Wilson action) may lead to large cut-off effects. It is thus favourable 
to use more sophisticated — improved — actions that suppress the lattice ar- 
tifacts. There exist several ways of constructing improved actions, in the next 
few sections we present the most important ones that are used for lattice gauge 
theories. 

One method to construct improved actions is following Symanzik's proposal 
[^Q, killing the lattice artifacts order by order in a^, adding operators of 
corresponding dimensions that cancel the artifacts. This yields e.g. the well 



known Liischer-Weisz action |39| which has leading corrections at 0{g'^a'*^). 

A technique very important in the context of this work is mean-field (or 
tadpole) improvement which helps to calculate perturbative coefficients, as those 
appearing in Symanzik improved actions, more reliably. This is achieved by 
changing the set up for lattice perturbation theory. On the lattice, gauge fields 
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are represented in the following way: 

2 2 

U^{x) = e^-aA,^-) = 1 + lagA^ix) - ^Al{x) + ■■■, (1.1) 

where U^{x) is the gauge link at lattice site x. This form obviously gives rise 
to local quark-gluon vertices with 1, 2, ... gluons. All these vertices except 
the lowest one are lattice artifacts. Contracting the two gluons in the term 
a^g^A^{x)/2 one obtains so-called tadpole diagrams. As it was recognised by 
Parisi and later by Lepage and Mackenzie |^l[| these artifacts are in fact 
only suppressed by powers of as the ultraviolet divergences generated by the 
tadpole loops kill the suppressions. These contributions are made responsi- 
ble for the poor match between short distance quantities and their perturbative 
estimates as well as for large coefficients occurring in perturbative lattice ex- 
pansions. 

Mean-field (or tadpole) improvement has been proposed to get rid of this 
kind of artifacts. It assumes that the lattice fields can be split into an ultraviolet 
(UV)and an infrared (IR) part and that the UV part should be integrated 
out 0,0: 

^^a^A,ix) ^ g»aff(Ajf ^ ygg^agAfCa,) ^ uoU^{x). (1.2) 

It just amounts to a simple rescaling of the link variables U^{x) by an overall 
constant factor uq < 1. Instead of the common expansion parameters, the 
coupling (7^, the hopping parameter k in the fermion action and the links U one 
uses their improved counterparts g^ = g^/ufj, k = kuq and U = U/uq. There 
are two common choices for the tadpole factor uq, either the fourth root of 
the plaquette expectation value or the expectation value of the link in Landau 
gauge. Tadpole improvement may be especially useful for the construction of 
anisotropic actions because it leads to a small renormalisation of the anisotropy 
^, if this quantity is introduced perturbatively, as e.g. in Symanzik improved 
anisotropic actions. Indeed, the anisotropic action most widely used [p|P-|2l|,p^, 
^ in the last few years, initially presented in [Q, combines mean-link 
and Symanzik improvement (at tree- level). 

Other, non-perturbative methods of improving actions base on the renormal- 
isation group. This concept, as well as resulting actions are briefly introduced 
in the next section. 

1.3 The Renormalisation Group 

A Field Theory is defined over a large range of scales from small, physical ener- 
gies up to the cut-off that is sent to infinity to reach the continuum limit. The 
fields associated with very high (unphysical) scales do influence the physical 
predictions indirectly through a complex cascade process; however no physical 
question involves them directly. This makes the connection between the local 
form of the interaction and the final predictions obscure. Additionally, the large 
number of degrees of freedom present makes the problem technically difficult. 
That is why one attempts to integrate out these non-physical high momen- 
tum scales in the path integral. The method accomplishing this, taking into 
account the effect on the remaining variables exactly, is called the Renormali- 
sation Group transformation (RGT) p3|-|50[|. 



1.3. The Renormalisation Group 
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C2, . . . 




The sequence of theories defined by the repeated use of the RGT on an initial 
theory defines a flow trajectory in the space of couplings. Fixed points (FPs) 
of this transformation are theories that reproduce itself under the RGT. Since 
the correlation length of the theory scales by the scaling factor of the RGT, 
its value has to be or oo at the fixed point. In Yang-Mills theory there is a 
non-trivial fixed point (the Gaussian FP) with correlation length — oo whose 
exact location in coupling space ci, C2, . . . depends on the RGT used. There is 
one so-called "relevant" coupling whose strength increases as one is starting at 
this FP and performing RGTs. The flow along this relevant scaling field whose 
end-point is the FP is c alled the renormalised trajectory (RT). This situation 



is sketched in Figure 1.1 Simulations performed using an action which is on the 
exact RT reproduce continuum physics without any discretisation errors. 

There are different approaches to take advantage of this knowledge about 
field theories in the vicinity of a fixed point, presented in the next two sections. 



1.3.1 Renormalisation Group Improved Actions 

Iwasaki pl| has studied RGTs in perturbation theory and has obtained a gauge 
action using two loop-shapes that after a few blockings comes close to the RT. It 
shows better rotational invariance and better scaling of Tc/ \fd than the standard 
action. Today, this action is used e.g. by the CP-PACS collaboration for MC 
simulations including domain- wall fermions as it turned out that the Iwasaki 
action shows better behaviour concerning chiral symmetry. 

Gupta and Patel |Q have studied the flow of a scale \fZ RGT and approx- 
imate the action using four parameters and two loop shapes in different repre- 
sentations. Studies of lattice topology have shown that short distance artifacts 



are suppressed 154,55 



The QCD-TARO collaboration |5q, p7| has examined the flow of a scale 2 
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RGT and evaluated the couplings on blocked lattices using Schwinger-Dyson 
equations. The blocked action is reproduced by two couplings corresponding to 
the simplest Wilson plaquettes 1x1 and 1x2. 



1.3.2 The Classically Perfect FP Action 

Another, more radical, method has been suggested by Hasenfratz and Nieder- 
mayer For asymptotically free theories like SU(A^) Yang-Mills theory, the 
action at the FP where the (only) relevant coupling g = Q (corresponding to 
(3 = oo) may be determined using a saddle-point equation. The FP action is 
classically perfect, i.e. there are no lattice artifacts on the solutions of the lattice 
equations of motion. This is equivalent to an on-shell Symanzik improvement 
at tree-level to all orders in a. 

Let us consider SU(iV) pure gauge theory in d = 4 Euclidean space-time on 
a periodic lattice. The partition function reads 

Z{(3) = / dUe-'^^'^^\ (1.3) 



where dU is the invariant group measure and A{U) a discretisation of the con- 
timmm action. We perform a real space RGT: 



g-/3U'(y) ^ 



J dU exp {-piAiU) + T{U, V))}, (1.4) 



where V is the blocked link variable and T(C/, V) is the blocking kernel defining 
the explicit form of the transformation. As for asymptotically free theories the 
FP lies at /? = oo, in this limit the path integral can be calculated in the saddle- 
point approximation. This leads to an equation in classical field theory defining 
the FP action which is mapped onto itself by the RGT: 

^FP(F) = min{^FP(t/) +r(C/,y)} . (1.5) 

1.4 The Construction of a FP Action for Coarse 
Configurations 



The FP saddle point equation 1.5 derived in the last section may be studied 



analytically up to quadratic order in the vector potentials [£9|,|60(] (see Section 



5.2). However for solving the FP equation on coarse configurations one has 
to resort to numerical methods. Additionally, although the full FP action is 
local, it is described in principle by infinitely many couplings that decrease 
exponentially with the space-time separation of the two coupled variables in the 
action. In order to use a FP action in a numerical simulation it thus has to be 
approximated by some parametrisation with a finite number of couplings. 



1.4.1 The Method Used 

The construction of the FP action on rough configuration relies on the obser- 
vation that the fluctuations of the minimising configurations U in eq. |1.5| are 
reduced by a factor of 30 to 40 compared to the initial coarse fields V . That 



1.4- The Construction of a FP Action for Coarse Configurations 
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is why after a small number of inverse RGTs (calculating the minimising con- 
figurations U corresponding to coarser configurations V ^ i.e. using eq. 



1.5 



as a 

recursion relation linking two actions A, A') the configurations are so smooth 
that they can be described accurately by a simple discretisation of the gauge 
action, e.g. the Wilson action. In practise it is not feasible to perform the 
necessary number of RG steps at once due to memory and computer time lim- 
itations, we therefore resort to a cascade process involving several single RG 
steps of first inverse blocking coarse co nfig urations and then parametrising the 
action on these configurations using eq. 



1.5, This idea is displayed in Figure [L 



In the following sections we describe the starting point of the procedure, i.e. a 
simple discretisation of the action for smooth fields, as well as the proceeding 
to parametrise actions Aiy) at intermediate and final levels, using eq. 



1.5 



Ao 







RG 
< > 


C/2 


RG^ 












Vo 









Ai 



A2 



v. 



A, 



Figure 1.2: The cascade procedure used to obtain perfect actions for strongly 
fluctuating fields V2 starting with a simple action ^0 suitable for very smooth 
configurations Uq. Coarse configurations V2 (on the level of MC simulations), 
Vi, Vo; fine configurations C/2, C/i, Uq (of which the latter may be described 
with the staring action Ao)- The boxes contain configurations on the same level 
of fluctuations. The actions A describing the respective configurations are also 
indicated. Note, that the number of steps actually used may differ from the 
number displayed in this Figure. 



1.4.2 The Starting Point 

As a starting point, we use very fine configurations (corresponding to "C/q" in 
Figure that are described well by a simple action. As an example, for the 
parametrisation of the isotropic action we generate coarse configurations ("Vb") 
where the plaquette variable u(Vb(nB)) < 0.25 for all ub and thus for the 
corresponding minimised configurations u{Uo{n)) < 0.009 for all n. These fine 
configurations can be accurately described by the quadratic approximation of 
the FP action fulfilling the 0{a^) Symanzik condition, denoted by ^0 0:11] ■ 

To check how well this action ^0 is working on the minimised configurations, 
we generate a set of configurations showing fluctuations of the same magnitude 
and corresponding configurations which are identical up to one link which is 
changed. The change of the action due to this changed link may then be calcu- 
lated on one hand using the parametrisation ^0 , on the other hand we minimise 

^Note, that the average plaquette variable for configurations occurring in MC simulations 
amounts to about 1.5 ~ 2, whereas the maximum values may reach 4.5 which is the upper 
bound for SU(3) matrices. 
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all the configurations and use the r.h.s. of eq. |1.5| which rests upon the action 
on the minimised configurations whose fluctuations are again smaller by a fac- 
tor of 30-40, i.e. this result is expected to be very close to the exact one. The 
difference of these two action values is thus a good estimate of the error of the 
parametrisation. It turn s ou t that the typical error is about 1-2% only. 

We may now use eq. L5 to fix the parameters of an appropriate parametri- 
sation of the action A\ valid on configurations with u{Vq) < 0.25. 

This process is now repeated with coarser configurations (l^i) whose min- 
imised counterparts (J7i) show fluctuations u{U i) < 0.25, describing these min- 
imised configurations on the r.h.s. of eq. |l.5| by the action Ai- It turns out 
that step by step the complexity of the parametrisation (i.e. the number of free 
parameters) has to be increased in order to be able to describe the action well 
on the coarser and coarser fields. 



1.4.3 Quantities to be Fitted 

Performing the numerical minimisation of the configurations V is an expensive 
task, even for small configurations. It is thus very useful to be able to fit not 
only one quantity per configuration (the total action A(V)) but to include as 
well the derivatives of the gauge action with respect to the gauge links in a given 
colour direction a: 

Ml^, M=l,-.-,4, a=l,... ,iV2-l, (1.6) 

where Nc denotes the number of colours. These quantities are easily accessible 
and deliver V ■ {N^ — 1) • 4 additional residues for the fit per configuration. It is 
understood that these quantities are not completely independent, but still the 
increase of useful data per configuration is dramatic. 

The intermediate actions Ai, A2 etc. still fulfill the O(a^) Symanzik con- 
ditions in order to be valid as well for configurations with smaller fluctuations 
— however, in the last step yielding the final parametrisation of the FP action 
for coarse fields, this requirement is dropped to get an effective FP action de- 
scribing well the typical MC configurations using a compact set of parameters. 
Additionally, we include scale-invariant instanton solutions |]6l|-|6^ to make the 
action performing well concerning topology. 

Adjusting the weights of the different quantities entering into the fit, we 
may tune the importance of the respective quantity in the determination of the 



action, see also Section 1.4.5 



In the last step we include into the fit several sets of 10-40 configurations 
each at different values of the coupling /3 to obtain a parametrisation which is 
valid on a certain range of couplings suitable for MC simulations. 

In the following Sections we describe the steps performed to parametrise a 
perfect action at a certain level of the cascade procedure. 

1.4.4 The Non-Linear Fit 

Before we may describe the actual fitting procedure used for the parametrisation 
of the isotropic FP action as well as for the perfect anisotropic actions, let us 
briefly explain some of the main features of our parametrisation for perfect 



actions. For a more detailed account on this issue, see Section 4.3 
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The action is a mixed polynomial of traces of simple loops (plaquettes) built 
from simple gauge links as well as from (APE- like) smeared links. The parame- 
ters in the smearing enter the fitting procedure non-linearly, the parameters in 
the polynomial enter the fit linearly. This distinction is crucial for the following 
description of the actual fitting process. The main advantage of this parametri- 
sation ansatz is the rich structure including a large number of different loops 
at moderate computational cost; however, the couplings of the loops in the ac- 
tion are not independent from each other but complicated combinations of the 
non-linear and linear parameters appearing in our parametrisation. 

First we perform a full non-linear fit varying the non-linear parameters 
describing the smearing as well as the linear parameters accounting for the 
composition of the action in terms of simple and smeared plaquettes. The linear 
part of the fit is performed exactly, its computational cost is negligible. The non- 
linear part is done using a Simplex algorithm which is rather slow, however the 
danger of missing the true minimum is small. A rough estimate of the number 
N of non-linear fitting steps (where every step includes the recalculation of all 
the residues) as a function of the number n^i of non-linear parameters is 

7V = 50(20)e°-23(4)"ni_ (17) 



Due to this behaviour {N increases roughly by a factor of 10 for 10 additional 
non-linear parameters), it is not easily possible to increase the number of free 
non-linear parameters above 15 ~ 20 in the full non-linear fit. 

The non-linear fit is performed for different sized sets of non-zero non-linear 
parameters with a (fixed) large number of free linear parameters. The only 
quantities fitted at this stage are the derivatives of the action with respect to 
the coarse couplings to obtain a "basic" value of which may be used to 
estimate the quality of the fits in the next step where additional quantities are 
taken into account or the number of linear parameters is decreased. Restrictions 
such as the norm of the action or Symanzik conditions are applied already to 
the non-linear fit, however. 

The number and the values of non-linear parameters of the parametrised 
action are fixed looking at the of the fit — which in principle is rather 
delicate, as it is unknown whether e.g. a 10 % decrease of when introducing 
an additional free parameter is a sign of importance of this extension or whether 
this set is already large enough to account for the FP nature of the action. 
These questions can be answered only at the end of the whole procedure when 
the final parametrisation is checked on physical quantities. However, whether 
the parametrisation is stable, i.e. whether the number of data points is already 
large enough to account for a parametrisation that does not depend on the 
actual configurations in the fit, may be — and has to be — checked after the 
non-linear fit, measuring x^ on independent sets of configurations. 



To give an example, we copy Table 1.1 from Section 4.3.4. In this case, the 



40% decrease of between the # = 1 and the # q = 3 non- linear parameter 
set (see Section 4.3 for an explanation) seems to be vital for an appropriate 
parametrisation of the perfect action as the rcnormaliscd anisotropy of the 
first parametrisation deviates considerably from the input value ^ = 2 whereas 
the renormalisation is small for the second, more sophisticated, parametrisation. 
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Ci 


i 


niax(fc + l)sp 


niax(fc + Otm 






4 


1 


4 


4 


4 


0.0250 


1.63(2) 


2 


3 


3 


4 


4 


0.0238 




4 


3 


3 


4 


4 


0.0144 


1.912(9) 



Table 1.1: Example for the behaviour of in the full non-linear fit (for the ^ = 2 
perfect action), depending on the number (and kind) of non- linear parameters. 
Values of the (physical) renormalised anisotropy measured later are also 
given. 



1.4.5 The Linear Fit 

Having at hand the set of non- linear parameters we may proceed performing fast 
linear fitting steps. The action values of the configurations as well as additional 
data such as the action of scale-invariant instanton solutions are included at this 
stage. The stability of is an indication of the ability of the chosen parameter 
set to incorporate additional information. At the same time the number of 
linear parameters may be decreased until starts to rise considerably. Figure 



1.3 displays the change of (for the derivatives and the action values of 20 
configurations at /? = 3.5 and (3 = 4.0, respectively) if one changes the weight of 
the action value residues compared to the (fixed) weight of the derivatives. The 
value of Wact finally used is marked by the dotted line. This choice preserves 
the good value of x^ for the derivatives up to a few per cent while of the 
action values is about twice the minimum value reachable with a very large 
weight. Without inclusion of the action values into the fit, the average relative 
error of the parametrised action compared to the true one amounts to about 
10%, with the weight chosen to about 0.5% and with infinite weight it might be 
decreased to about 0.2% (destroying the good parametrisation of the derivatives 
completely). 

However, in the last step of the cascade process, when one aims at para- 
metrising the final action to be used in MC simulations, the main difficulty 
of the choice of the linear set to be used is the behaviour of the parametrised 
action in terms of the simple and smeared plaquettes of the configurations, 
determined by the linear parameters. It happens that all the data in the fit is 
properly described by the parameters, but that varying the simple or smeared 
plaquette value slightly above or below the values occurring in the fit results in 
fake negative contributions to the action. In a MC simulation, this is a very 
dangerous behaviour as this might result in configurations with an artificial 
structure, i.e. the algorithm could try to maximise the number of plaquettes 
having a negative contribution to the action etc. One thus has to study the 
linear behaviour of the action, e.g. looking at plots of the action contribution 
vs. the values of the simple plaquette u (which is bounded, < u < 4.5) and 
the smeared plaquette w (which in principle is unbounded as w is not projected 
to SU(3)). It may be that accidentally none of the parametrisations behaves 
well enough for MC simulations. This is especially the case if conditions on the 
linear parameters, such as the norm of the action, that are determined by the 
values of the non-linear parameters lead to the bad behaviour of the linear part 
of the action. In this case we include a number of constraints A{u^ w) > Q (for 
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Figure 1.3: The effect of changing Wact, the weight of the action compared to 

the weight of the derivatives, in the Unear fit for the £, = 2 perfect action. 
The values of for the derivatives of /3 = 3.5 (circles) and /? = 4.0 (boxes) 
configurations as well as for the action values of /3 = 3.5 (pluses) and /? = 4.0 
(crosses) configurations arc given. The dotted line marks the final choice of the 
relative weight Wact • The number of parameters remains fixed. 
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different pairs of {u, w) at tfie edge of the (u, i(;)-region expected to occur in MC 
simulations) into the non-linear fit in order to get a different set of non-hnear 
parameters allowing for linear parameters with a good behaviour. (Generally 
in the space of non-linear parameters several almost degenerate minima in 
exist for very different sets of parameters, so this does not result in a significant 
rise of (of course, this has to be checked). It turns out that the best way 
of introducing these constraints into the non-linear Simplex minimisation is 
adding the absolute value of negative contributions to the action at the {u, w) 
pairs multiplied with a large factor (we use 10^). "Hard" constraints with a 
step function at A{u, w) < may lead to failure of the simplex algorithm as in 
this case it is not directed towards the region of the parameter space, where the 
action is behaving as desired. 

Figures 1.4 and 1.5 show linear sets behaving badly and well. These examples 
are taken from the fit of the ^ = 4 perfect action (see Chapter 0) where it is 
necessary to force several points (u, w) in the non-linear fit to result in a positive 
contribution to the action. The first figure (standard non-linear fit) shows that 
in this case the spatial linear parameters lead to negative contributions of spatial 
plaquettes in a large region of {u, w) values whereas the temporal part of the 
action is o.k. The second figure shows the behaviour of the parameters finally 
chosen, where conditions on the action contributions are included into the non- 
linear fit. There are no dangerous regions with negative contributions to the 
action. Pilot MC runs prove that indeed no abnormal behaviour of the action 
is observed. 
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Figure 1.4: A linear set of parameters with bad behaviour (danger of generating 
artificial configurations with a lot of negatively contributing plaquettes in a MC 
run). Upper plot: Spatial parameters, lower plot: temporal parameters. 




Figure 1.5: A linear set of parameters with good behaviour. Upper plot: Spatial 
parameters, lower plot: temporal parameters. 



Chapter 2 



Physical Objects in Pure 
Lattice Gauge Theory 



There are several interesting objects existing in pure gauge theory that can be 
studied non-perturbatively on the lattice. Some of them, the ones we study 
by Monte Carlo simulation, are presented in this chapter, together with ex- 
perimental results (if present). They comprise the torelon (including non-zero 
momenta), the deconfining phase transition, the static quark-antiquark poten- 
tial and glueballs. Furthermore, we explain how these states are measured on 
the lattice and what techniques we use to determine the energies of the states 
or the critical temperature of the phase transition respectively. For the results 
of the simulations we refer to Chapter |^ for a recapitulation of results obtained 
using the isotropic FP action, to Chapter || for ^ = as/ at = 2 perfect action 
results and to Chapter for first results obtained with the ^ = 4 perfect action. 



2.1 The Torelon 



2.1.1 Torelons in the Continuum 

The torelon is a closed gluon flux tube encircling the periodic spatial boundary 
of size L of a torus. It is created and annihilated by closed gauge operators 



winding around this boundary, see Figure 2.1 



S^xT 




Figure 2.1: A creation operator for torelons. 
String models which are a good approximation for long flux tubes predict 
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its energy to behave as 

E{L)^aL + ^, (2.1) 

where a is the string tension and the second term describes the string fluctuation 
for a bosonic string [|65|-|68t. In contrast to the static quark- antiquark potential 
(see Section |2.3| ) where the string models the gluon flux between static colour 
sources, there is no Coulomb part corresponding to gluon exchange.]^ The string 
models assume no restriction on the transverse modes of the string, the trans- 
verse extension S of the volume thus has to be large enough in order to be able 
to apply their predictions. 

To learn more about QCD in a finite volume, see 



2.1.2 Torelons on the Lattice 

On the periodic lattice, a torelon associated to a point in the (x, ?/)-plane is 
created by a Polyakov loop winding around the lattice in z direction starting 
at that point. Because of the relatively small lattice sizes necessary and the 
rather simple operators that can be used, torelon measurements on the lattice 
have quite a long history ||70|-[79|]. In order to be able to compare the results to 
string models as well as not to have too large transverse momenta we choose to 
work on an S*^ X i X T lattice where S and T are large, whereas the length of 
the z-direction L is chosen to be rather short, because long torelons are more 



difhcult to measure due to their higher energies (see eq. 2.1). 

As the colour flux is wandering as it loops around the lattice, it is useful to 
employ (iteratively) APE smeared links |Q to model this, which improves the 
overlap of the operators with the ground state. One step of APE smearing acts 
on a gauge link Uj (x) as 



S,U,{x) = Vsxj(3){u,{x) + KY,{Uk{x)Uj{x + k)ul{x+j) (2.2) 

+Ul{x - k)U,{x - k)Uk{x - + j))}, 



the original link variable is replaced by itself plus a sum of the four neighbouring 
spatial staples and then projected back to SU(3). The smeared and projected 
links SiUj{x) have the same symmetry properties under gauge transformations, 
charge conjugations or reflections and permutations of the coordinate axes. The 
smearing is used iteratively allowing to measure operators on configurations of 
different smearing levels, SnU . 
We thus measure the operators 

L-l 

Tn{x,y,t) = Tr ]^ SnUz{x,y, z,t), (2.3) 

z=0 

for all X, y, t, where SnUz is the n times iteratively smeared link in the longitu- 
dinal z direction. By discrete Fourier transform, we project out the state with 



^Note, that despite the absence of the Coulomb part, extracting the string tension from 
torelon measurements is difficult. This is due to the large torelon energies for long strings 
(large L) where the relevance of the string model term 7r/3L decreases. 
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momentum p = {Px,Py) = {nx,ny){2'K/Sas): 

T„(p, t)^Yl ^-(^' V' t)e'^'"'^+Pyy^ (2.4) 



and build the correlators 

Cnn'{p,t) = V(r„(p,T)T„,(p,r + t)). (2.5) 

r 

Using these correlators of operators measured on the different smearing levels we 
obtain the energy values atE{p) employing variational methods (see Appendix 
The continuum dispersion relation 

= p2 + (2.6) 

on the lattice becomes (in temporal units) 

(2.7) 



where the components of the (transversal) lattice momentum. On an 

anisotropic lattice, this equation allows for the extraction of the renormalised 
anisotropy = as /at (measuring the "renormalisation of the speed of light") 
as well as the torelon mass mat, which in turn may be used to get an estimate 
of the scale using eq. 2.1 and known values of the string tension a. 



2.2 The Deconfining Phase Transition 

Pure lattice gauge theory is invariant under a global unitary transformation 
which transforms all temporal links U4 {x, t) of a given timeslice (fixed t) as 

Ui{x,t) — > zUi{x,t), (2.8) 

where z is an element of the center Z{Nc) of the gauge group SU(A^c)- For 
SU(3), z = exp(27riZ/3), / =0, 1, 2. 

To study this symmetry, we look at the Polyakov loop (or Wilson line) 



Nt-l 

F(.f)=Tr [] C/4(f,i), (2.9) 
t=o 

which is a loop wrapping around the lattice in temporal direction, representing 
a single static quark at finite temperature, as the temperature T is proportional 
to the inverse of the temporal extension Nt, T (x 1/Nt. 

The correlation of the Polyakov loop P{x) and its adjoint ^'^(O) is related to 
the free energy Fqq of a static quark separated from a static antiquark by the 
distance x according to 



(P(f)pt(o)) cx exp(-i^,5/r), 



(2.10) 
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where T denotes the temperature. For large separation, cluster decomposition 
requires 

(P(f)pt(o)) ^|(F(0))|2. (2.11) 

For low temperatures below the critical temperature, T < Tc, center symme- 
try holds and thus (P(0)) = which allows for increasing Fqq(x) with growing 
distance the static quarks are confined. For temperatures above Tc, center 
symmetry is broken, (P(0)) and thus the free energy of the static quarks 
has to vanish, Fqq{x) = 0, the quarks are free and form, together with the gauge 
degrees of freedom, the quark-gluon plasma. 

Figure displays the distribution of the Polyakov loop for different cou- 
plings P (for the ^ = 2 perfect action, see Chapter ^. The critical couphng on 
the lattice used has been determined to (3c = 3.118(1). It is clearly visible that 
below the critical coupling the Polyakov loop values are close to zero whereas 
above (3c they choose one of the three degenerate vacua. The coexistence of the 
two phases may be observed in the plot for P — 3.12 which is very close to (3c. 

The critical temperature Tc may be determined by the location of the peak 
in the susceptibility 

X^N!{{\P\^)-{\P\r). (2.12) 

There is strong evidence for this transition to be first order for SU(3) |]8l|-^, 
however in finite volume (on lattices used in MC calculations) the 5-peak of the 
susceptibility x is washed out and limits the precision with which one may 
determine the critical temperature Tc or, equivalently, the critical coupling Pc 
corresponding to a fixed temporal extension Lt = aNt of the lattice. 



2.2.1 Determining the Peak of the Susceptibihty 

As mentioned in the last section, the critical coupling Pc is determined by locat- 



ing the peak of the susceptibility x (eq. 2.12 ). To accomplish this, simulations 
at different values of P near the estimated critical coupling Pc have to be per- 
formed. For all these simulations a certain statistics has to be reached and it 
may thus be very time consuming to approach Pc closer and closer performing 
simulations at more and more values of /?. As the Monte Carlo simulation may 
exhibit critical slowing down, i.e. a dramatic rise of the autocorrelation times 
(see Section ^.iD near the phase transition, the procedure gets even more ex- 
pensive. A nice solution to this problem is the use of reweighting techniques 
that allow extrapolating or interpolating results obtained at a given coupling or 
temperature to other (rather close) values of these parameters. To determine 
the peak of the Polyakov loop susceptibility we employ Ferrenberg-Swendsen 
reweighting and thus present the main idea of this procedure very briefly 

(for a more extensive presentation, see Appendix C of ||l|). 

All the information about a statistical system at temperature T is contained 
in the partition function 

Z(/3)=^e-^^(-), (2.13) 

where {a} denotes the set of all possible configurations of the system, S{a-) is 
the energy for a given configuration a. The same partition function may be 



2.2. The Deconfining Phase Transition 



21 



|)=3.07 



-0.4 -0^ 0.2 0.4 0.6 

ReP 



|)=3.11 



-U.4 -U.2 



|)=3.15 



Figure 2.2: The distribution of the Polyakov loop at different couphngs /? for 
the ^ = 2 perfect action, Nt = 5, Ns = 9. The critical value of the coupling 
has been determined to f3c = 3.118(1). The system is heated up such that the 
center symmetry is broken as /3 > /3c- For clearness, only (randomly chosen) 
10% of the data are displayed. 
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written as 

Z(/3)=^W^(5)e-^^ (2.14) 
s 

where W{S) is the density of states at the energy S, also caUed spectral density 
function. This function is universal, it is the same function for every temperature 
and thus contains in principle all the information about the system at any 
temperature T . In practice, our numerical results at a given value of [3 will 
allow us to estimate W{S) only for a finite range of energies S occurring often 
enough in simulations at the given coupling. However, if we are studying a phase 
transition, due to the broad probability distribution of the examined states in 
simulations near criticality, /3 being close enough to /3c, it is possible to estimate 
the spectral density function in a large range of energies. 

Let us now consider n measurements of an observable obtained in a numerical 
simulation at coupling /3. We plot a histogram of the action values corresponding 
to the measurements and denote the number of measurements corresponding to 
one point of the histogram by N{S). Trivially, J2s^(^) ~ '^here S now 
denotes a certain range of energies binned together to build one point of the 
histogram. We can estimate the probability of a configuration to have energy 
S at coupling /3: 

Z{P) np 

Using this we have an estimate of the spectral density function W{S) of the 
system: 

W{S) « M^z(/3)e'5^. (2.16) 

Suppose now, we have performed MC runs at K different values of (3 (/3i, 
/32, ■ • ■ , Pk) where the frequencies Nk{S) have been measured. For every run 
we estimate the spectral density function 

Wk{S) « ^^ZiPk)e^>'^. (2.17) 
nk 

We know however that W{S) is universal, there should thus be a unique func- 
tion W{S). Using the measured estimates of the functions Wk{S), this unique 
function can now be constructed in an optimal way following the procedure 
by Ferrenberg and Swendsen |8^, |8^, which we do not describe here. Once 
we know this unique function we may calculate the partition function at an 
arbitrary value of /3: 

Zij3)^Y.W{S)e~^^. (2.18) 
s 

We are now ready to estimate the expectation value of any observable at an 
arbitrary coupling /? using 

{0)p = 1^ E 0{S)W{S)e-f'''. (2.19) 
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Figure ^ shows the action (energy) histograms for four values of the cou- 
phng /3 in the vicinity of the phase transition. At each couphng, 45000-55600 
measurements have been performed. Note the double peak structure which is 
prominent for the /3 — 3.03 distribution (solid line) closest to f3c = 3.032(1) and 
which starts to evolve for the (3 = 3.04 distribution (dashed line). Figure 2.4 
displays the corresponding result of a reweighting procedure for the Polyakov 
loop susceptibility. Notice, that the error of the peak position (which corre- 
sponds to the critical coupling f3c) is small, in this example the error of Pc is as 
small as A/3c = 0.001. 



< 




Figure 2.3: The action (energy) distribution at 4 different values of /3c on the 
4 X 8^ lattice for the ^ = 2 perfect action. Dotted: [3 = 3.02, solid: [3 = 3.03, 
dashed: (3 = 3.04, dot-dashed: (3 = 3.05. 



Figure 2.5 displays the MC history of the angle of the Polyakov loop arg(P) 
for (3 close to (3^. It becomes apparent that most of the information about the 
phase (and thus about the Polyakov loop susceptibility xp) is contained in the 
number of flips between the phases as well as in the time the configuration 
is staying in one of the phases. An important quantity in this respect is the 
persistence time Tp which is the MC time (number of sweeps) divided by the 
number of flip-fiops (change to another phase and another change back to the 
original phase). Because of critical slowing down, the autocorrelation time of 
the operators has to be considered as well. Near the phase transition knowledge 
about these quantities is crucial to be able to estimate the quality of the statistics 
given the number of sweeps performed. 
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3 3.02 3.04 3.06 



Figure 2.4: The Polyakov loop susceptibility on the 4x8'^ lattice for the 
^ = 2 perfect action. The solid curves are the interpolations using Ferrenberg- 
Swendsen reweighting, the dashed lines indicate the bootstrap error. 




MC time 

Figure 2.5: MC history of the Polyakov loop measurements on a 5 x 9'^ lattice at 
/? = 3.12 (/3c=3.118). The argument of the Polyakov loop P is given, arg(F) = 
0, ±27r/3 for a configuration being in the deconfined phase. 



2.3. The Static Quark- Antiquark Potential 



25 



2.3 The Static Quark- Antiquark Potential 

The potential between a quark and an antiquark, both carrying colour charge, 
may be studied even in pure gauge theory, creating and annihilating a static 
quark and a static antiquark in the gauge background. As this potential is 
rather easy to measure and as it may be compared to experimental data ob- 
tained by studying bound states of a quark and an antiquark, the heavy mesons, 
it provides an excellent mean of setting the scale of lattice simulations in the 
absence of dynamical quarks (i.e. fermionic theories in the quenched approxi- 
mation or pure Yang-Mills theory). The scale is needed if lattice results, such 
as masses, should be converted into physical units. In this Section, we will thus 



present sources of experimental data coming from light mesons in Section 2.3.1 



and from heavy mesons in Section 2.3.2 as well as how the potential may be 
examined in pure lattice gauge theory (Sections ^.3.3 , 2.3.5) and how it is used 



to set the scale in Section 2.3.4 



2.3.1 Light Mesons 

The elementary constituents of hadronic matter, the quarks, carry colour charge 
which enables them to build bound states, held together by the strong force. 
The simplest such state consists of a quark and a corresponding antiquark. Since 
the early sixties it has been noticed that these states, the mesons, as well as 
bound states of three quarks, the baryons, having mass m and spin J group 
themselves into almost linear, so-called Regge trajectories JSSf-pOt, i.e. if one 
plots w? vs. J the connection is linear. This can be observed up to spins as 
high as J = 11/2. The behaviour is thus mainly 

J(m) = a(0) +a'm2, (2.20) 
where a(0) is known as the Regge intersect and 

a' = (2.21) 

271(7 

is the Regge slope with cr, the "string tension" . The explanation for this name 
stems from a simple model that explains the linearity of the Regge trajectories: 
Imagine a rotating string of length 2d with a constant energy density per unit 
length a. If this string spans between (almost) massless quarks these quarks 
will move at (almost) the speed of light c = 1 with respect to the center of mass. 
The velocity at a distance r of the centre is thus given by v{r) — r/d. From this 
we may easily calculate the energy stored in the string, 

drcr 

= TTda, (2.22) 

Jo ^l-i;^(r) 

and the angular momentum, 

Jo VI -^^'W 2 2^(7 ' ■ ' 



which is the relation between Regge slope and string tension in eqs. ^.20 , ^.21 



The experiment predicts values of 470 MeV < ^/a < 480 MeV for trajectories 
starting with a pseudo-scalar (tt or K) and values around 430 MeV for the other 
trajectories. 
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2.3.2 Heavy Mesons 

After the discovery of the J /ip meson consisting of a charm quark and a charm 
anti-quark in e+e^ annihilations it was suggested to treat such states buih out of 
heavy quarks non-relativisticaUy ]9l|] . Because of the analogy to positronium in 
electrodynamics these states have been named quarkonia. They are composed 
of charm or bottom quarks, the top quark does not appear as a constituent 
because its weak decay rate t bW^ is large (see @). The best studied 
states are cc and bh, less is known about bound states of a 6 and a c quark. The 
string tension appearing in the Regge slope may be determined experimentally, 
e.g. through the mass splitting between different states of the same constituents, 
e.g. IP — IS for bb which is Xb ~ T. However, these measurements effectively 
probe a range of about 0.2 ^ 1 fm where it is still difficult to extract the string 
tension a unambiguously. 

For sufficiently heavy quarks the characteristic time scale of the relative 
movement of the quarks is much larger than the one of the gluonic degrees of 
freedom. In this case we may apply the adiabatic (or Born-Oppenheimer) ap- 
proximation and describe the effect of gluons and sea quarks by an effective 
instantaneous interaction potential between the heavy quarks. In this approx- 
imation quarkonia are the positronium of QCD. However, unlike in QED, it is 
not possible to calculate the potential perturbatively and predict the spectrum 
but we have to solve the inverse problem, namely guess the form of the potential 
from the observed spectrum (and decay rates). The Cornell potential ]9^ , 

V{R) = VQ + -+aR (2.24) 
K 

describes the interaction reasonably well in the energy range probed. Further- 
more, there is not enough structure in the measured potential to require a more 
complex parametrisation. Recent results for the parameters a and a (Vq is 
an unphysical constant) are a ~ —0.51, ^/a ~ 412 MeV and a w —0.52, 
^/a « 427 MeV |9^]. However in the range 0.2 fm< r < 1 fm which is probed 
by quarkonia splittings these parametrisations do not differ significantly even 
from earlier values a « —0.25, ^/a w 455 MeV |9^] because the higher value 
of the Coulomb coefficient is compensated for by a smaller slope a. Still, it is 
very interesting to see the compatibility of these values with the estimate for 
the string tension a from Regge trajectories of light mesons. 

For an extensive review about QCD potentials and their study on the lattice 
see ref. Q. 

2.3.3 The Static Quark- Antiquark Potential on the Lat- 
tice 

In lattice gauge theory, the rectangular Wilson loop W(i?, T), having extensions 
R and T in spatial and temporal directions respectively, starting at the point 
{x,t), creates at time t a static pair of a quark q and an antiquark q sitting 
at X + R and x and annihilates it again at time t + T. The corresponding 
static potential between the colour sources is related to the expectation value 
of >V(i?,r) as 



{W{R,T)) (X exp{-V{R)T) 



(2.25) 



2.3. The Static Quark- Antiquark Potential 
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for large T. The potential can thus be calculated in principle by 

y(i?)= hm -iln(W(i?,r)). (2.26) 

1 — >oo 1 

In pure lattice gauge theory the potential is confining the quarks and from strong 
coupling expansion it follows that the Wilson loop respects an area law 

(W(i?, T)) cx exp{-<7RT) for i? -> oo, (2.27) 

where a is again the string tension between the static quark sources. 

Due to the asymptotic freedom of QCD, perturbation theory is reliable 
at short distances. So for small R, it predicts a Coulomb-like interaction 
V{R) cx a/R between the quark and the antiquark. The most simple ansatz 
to describe both these properties of the gg-potential is the Cornell potential, 
eq. [2.24^ mentioned in the last section. All our measurements of the Wilson 
loop will be fitted to this ansatz. 



2.3.4 Setting the Scale Using the Static qq Potential 



Using the string tension a whose value is in principle known experimentally, 
one may set the scale corresponding to the coupling (3 used in the simulation. 
However this is plagued by several problems: firstly, it is not very easy to 
determine the long-range quantity a on the lattice due to a poor signal/noise 
ratio for large values of R, the value of a is reached asymptotically which leads to 
the demand for large (and thus computationally expensive) lattices, and finally, 
the string tension a is not well defined in full QCD because if the energy of 
the string is large enough it may break and create a light quark-antiquark pair 
("string breaking"). Due to these problems, the hadronic scale rg - the so-called 
Sommer scale - has been introduced through the force F{R) between the 
two quarks having an intermediate distance 0.2 bn < r < 1.0 fm. There are 
several advantages of the hadronic scale: Firstly, it can be measured much better 
on the lattice than the string tension a, secondly, in this intermediate distance 
there is reliable information about the physical scale from phenomenological 
potential models ||9S , 100 |. The scale is defined like 



dV_ 
dR 



(2.28) 



where originally c — 1.65 has been chosen. The corresponding experimental 
value from bottomonium phenomenology 101, 102| is r^^ = (394 ± 20) MeV 
and thus ro ~ 0.50 fm. However, it may be useful to measure the hadronic 
scale at different distances corresponding to other values of c, that is why we 
have collected values for c and from high-precision measurements of the static 



potential performed with the Wilson action [103-105| listed in Table 2.1 



2.3.5 Measuring the Static Potential on the Lattice 

The correlation functions of the strings in MC measurements are always con- 
taminated by high-momentum fiuctuations. The common way to reduce these 
excited-state contaminations is smoothing the spatial links employing APE 
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0.662(1) 


0.89 


1.00 


1.65 


1.65(1) 


4.00 


2.04(2) 


6.00 



Ta bic 2. 1: Parameter values for the determination of the hadronic scale through 
eq. |2^. 



smearing (see Section 2.1.2). The Wilson loop operator Wnm{R^T) is con- 
structed as a product of iteratively smeared spatial links on smearing level n on 
time slice t and smearing level m on time slice t-\-T connecting two spatial points 
of distance R and unsmeared link products of length T in temporal direction 
(see Figure 2.6). The Wilson loops easiest to measure are the ones parallel to 
the spatial axes — in this case the choice of the shortest paths is unique and the 
number of paths which could be taken into account and have to be measured on 
the lattice is small. However, taking into account separations R not parallel to 
the spatial axes as well has got the advantages that the discretisation of the sep- 
aration \R\ is not as coarse-grained, i.e. quantities like tq that rely on differences 
on the lattice can be determined more accurately and more reliably, furthermore 
the energies of off-axis separated quark-antiquark pairs allow to estimate effects 
due to violations of rotational invariance. This is why in some runs we include 
in the measurements separations that are multiples of the lattice vectors (1,0,0), 
(1,1,0), (1,1,1), (2,1,0), (2,1,1), (2,2,1) (and lattice rotations). The spatial paths 
are constructed such that all the shortest paths to the vectors are calculated 
first, out of these initial six vectors the longer separations are constructed. 

Averaging the operators over the three spatial directions as well as over the 
whole lattice we get the zero momentum correlation matrix that can be analyzed 
using the variational method described in Appendix 



Uwwvwvwvv 

R 



Figure 2.6: The smeared Wilson loop operator Wnm{R,T). 

This method yields projectors to the ground state of the string for each R 
which can be used in a correlated fit to the phenomenological ansatz for the 
potential, eq. 2.24, where R= \R\. The hadronic scale Vc/a is then determined 
by local fits to V{R) including 2 ~ 4 separations and converted to r^/a using 
Table 2.1 if necessary. Depending on the lattice spacing this can be done for 
several values of c yielding information about systematic uncertainties. 
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2.4 Glueballs 

2.4.1 Introduction 

The particles mediating the strong interaction of QCD, the gluons, carry colour 
charge and thus interact with each other, unlike, e.g., their counterpart in elec- 
tromagnetism, the photons, which have zero electric charge. This fact has also 
been established experimentally, studying 4-jet events in e~^e~ annihilation. 
The spectrum of QCD may thus contain bound states of (mainly) gluons, called 
glueballs. These states are described by the quantum numbers J denoting the 
(integral) spin, P denoting the eigenvalue ±1 of the state under parity and C, 
denoting the eigenvalue ±1 under charge conjugation. Thus the eigenstates of 
the Hamiltonian corresponding to glueball states are labeled | J''^'^). 

Due to their non-perturbative nature, glueballs can be theoretically studied 
most reliably doing numerical simulations of QCD on a space-time lattice. As 
(pure) glueballs are composed entirely of glue, the study using the quenched 
approximation of QCD (where fermions, if present, are infinitely heavy) makes 
sense .0 

Experimentally, there is evidence for exotic glueballs or hybrid particles con- 
sisting of quarks and gluonic excitations. These states (often called "oddballs" ) 
have exotic (or "odd") quantum numbers (e.g., 0"' or 1 '") and cannot be 
explained by simple bound states of purely quarks. Due to this, they cannot 
mix with mesons and are thus particularly interesting to study in pure gauge 



theory. However, these states are found to be high lying [|18| , |106| - |108|| (above 
twice the mass of the lightest glueball observed) . This may lead to ambiguities 
in the analysis of the experimental results because these states may mix with 
bound states of hghter particles. 

The lighter glueballs with conventional quantum numbers are difficult to 
distinguish from the dense background of conventional meson states observed 
in experiments. They are expected to be created in "gluon-rich" processes, 
such as radiative J/ip decays, central production (two hadrons passing each 
other "nearly untouched" without valence quark exchange) or pp-annihilation. 
Currently, there is an ongoing debate whether light glueballs (above all the 
scalar O"*""*" which is the lightest state in pure lattice gauge theory with a mass 
of about 1.6 GeV) have been observed experimentally at about the mass that is 
predicted by quenched simulations on the lat tice, whether the lightest glueball is 
much lighter (below 1 GeV) and ver y bro ad |l09| , or whether glueballs have not 
been observed at all in experiments |llC| . There are mainly two reasons for this 
uncertainty. On one hand, the experimental data seem not yet to be accurate 
and complete enough, despite large efforts in the last years, driven by the lattice 
results; on the other hand, lattice simulations with high statistics, measuring 
glueball states, have been performed only in the quenched approximation, where 
the quarks are infinitely heavy and thus static. Decreasing the sea (dynamical) 
quark mass (finally down to the physical value) will allow to track the glueball 
states as sea quark effects are increased. It may turn out, that indeed the 
glueball mass is lighter than the one measured in pure gauge theory (for partially 



^The question, whether pure Yang-Mills Theory has a mass gap is itself a very important 
problem. To prove this rigorously is considered one of the seven Millennium Problems, for- 
mulated by the Clay Mathematics Institute (www.claymath.org); the mathematical proof is 
awarded $1 million. 
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quenched results possibly indicating this see |111-113|). It may even happen 
that by "switching on" the sea quarks the scalar glueball acquires a very large 
width and thus decays (almost) instantaneously to qq states, i.e. it ceases to exist 
physically. However, the (partially) unquenched results are rather indecisive yet. 

Following the explanation that is dominant at this time, the observed low 
lying scalar mesons /o(1370), /o(1500) and /o(1710) are all mixtures of glue 
and mesonic components as dd, uu and ss. The proposed decompositions of the 



wave functions into these contributions [114-121| are rather different, however 
they share some robust common features. Studies in pure gauge theory may 
help separating the pure gauge part from the mesonic part of these states. 

2.4.2 Glueballs on the Lattice 
States 

Glueballs in the continuum are rotationally invariant and have a certain (inte- 
gral) spin J. On the lattice, the rotational 0(3) symmetry is broken, only its dis- 
crete cubic subgroup Oh survives the discretisation. Therefore, the eigenstates of 
the transfer matrix are classified according to the five irreducible representations 
of Oh'- ^1, A2, E, Ti, T2 with dimensions 1, 1, 2, 3, 3 respectively. Their trans- 
formation properties may be described by polynomials in the components j/, 
z of an 0(3) vector as follows: Ai ~ {!}, A2 ~ {xyz}, E ^ {x^ — z^, — z^}, 
Ti ^ {x, y, z}, T2 ~ {xy, xz, yz}. Generally, an 0(3) representation with spin J 
splits into several representations of the cubic group. Since Oh is a subgroup of 
0(3), any representation Dj with spin J in the continuum induces a so-called 
subduced representation Dj | Oh on the lattice. This subduced representation 
no longer has to be irreducible but is a direct sum of irreducible representations 
TP of Oh-. 

Dj i o,, = © © . . . . (2.29) 



Table 2.2 lists the subduced representations of Dj for J = 1, . . . ,6. The spin 
J = 2 state for example splits up into the 2-dimensional representation E and 
the 3-dimensional representation T2. Approaching the continuum, rotational 
symmetry is expected to be restored and thus, as a consequence the mass split- 
ting of these two states will disappear and the two representations form together 
the 5 states of a spin J = 2 object. 

Operators 

Pure glue physical states on the lattice are created and annihilated applying 
gauge invariant operators to the pure gauge vacuum. In our simulations, we 
use space-like Wilson loops in the fundamental representation of SU(3). Since 
we do not aim at measuring non-zero momentum glueballs we consider only 
translationally invariant operators, i.e. operators averaged in space. 

It is computationally feasible to measure Wilson loop operators up to length 
8. The composition of the irreducible representations r^*-^ of the cubic group in 



terms of these 22 loop shapes has been done already in Ref. Il23]. Figure 2.7 dis- 



plays all Wilson loop shapes up to length 8 together with the numbering which 



will be used in the forthcoming sections. Table 2.3 lists the number of orienta- 



tions of the different loop shapes. Tables 2.4, 2.5 list the irreducible contents 
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Table 2.2: The composition of the subduced representation Dj I Oh in terms 
of the irreducible representations of the cubic group Oh- 



of of the C = + and C = — representations of the cubic group, respectively. 
Certain operators contribute to the two- or three-dimensional representations 
with two or three different polarisations, in analogy to different magnetic quan- 
tum numbers m for a given angular momentum / in the 0(3) group. Measuring 
all these polarisations may suppress statistical noise more than just increasing 
statistics since the different polarisations of a loop shape are expected to be anti- 
correlated. Examining the measured correlators of different polarisations of the 
same representation, it turns out that in most cases one of the polarisations is 
measured very well, whereas the others exhibit a bad signal/noise ratio. 



loop shape # 
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12 
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96 



Table 2.3: Dimension d of the representation of Oh 'x C on the loop shapes, i.e. 
number of different orientations. 



The accuracy of the measurements of different operators may differ consid- 
erably; firstly, because shapes with larger multiplicity show smaller statistical 
noise, secondly, because operators on different smearing schemes may exhibit 
less or more fluctuations. To give an impression of the differences between var- 
ious operators, Tables ^ display correlator "lifetimes" obtained from the 
glueball measurements at coupling f3 = 3.3, with a temporal lattice spacing 
at « 0.08 fm « (2.4 GeV)^^. The correlator lifetime Tcorr is (conventionally) 
defined such that at time t — Tcorr the relative (bootstrap) error of the corre- 
lator is 25%. The relative error is linearly interpolated between t < Tcorr and 
t -I- 1 > Tcorr- Table 2.6 compares the shapes contributing to different glueball 
representations, measured on smearing level 1 (3 times smeared) at /3 = 3.3, 
Table 2.7 compares the same quantities measured on smearing level 3 (9 times 
smeared). For shapes which contribute with more than one projection to a given 
representation, the best measured projection is displayed — usually the other 
polarisations are measured much worse due to the anti-correlation. 

The operators measured on different smearing levels are again used in the 
variational method described in Appendix ^ The scalar representation A'^'^ 
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Figure 2.7: All space-like Wilson loop shapes up to length • 



2.4- Glueballs 



33 



loop 


1 


2 


E++ 


T++ 
1 


T++ 


A7+ 
1 


2 


E-+ 


T-+ 
1 


T-+ 

2 




1 





1 

























1 


1 


2 























#3 

IT'-' 


1 





1 





1 











1 


1 


#4 


1 











1 



















1 





1 

























1 





1 





1 











1 


1 


tt ' 


1 





1 





1 











1 


1 


#8 


1 


1 


2 


3 


3 


1 


1 


2 


3 


3 


#9 


1 





1 


1 


2 


1 





1 


1 


2 


#10 


1 


1 


2 























#11 


1 


1 


2 


1 


1 











2 


2 


#12 


1 


1 


2 


1 


1 

















#13 


1 





1 




















1 


#14 


1 





1 


1 


2 

















#15 


1 





1 























#16 


1 


1 


2 























#17 


1 





1 





1 

















#18 


1 





1 


1 


2 


1 





1 


1 


2 


#19 


1 





1 





1 

















#20 


1 





1 


1 


2 


1 





1 


1 


2 


#21 


1 





1 





1 











1 


1 


#22 


1 


1 


2 


3 


3 


1 


1 


2 


3 


3 



Tabic 2.4: Irreducible contents of the C-parity plus representations of the sym- 
metry group of the cube on Wilson loops up to length eight. 
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Table 2.5: Irreducible contents of the C-parity minus representations of the 
symmetry group of the cube on Wilson loops up to length eight. 
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Table 2.6: "Lifetimes" of glueball correlators at /3 = 3.5 (first smearing scheme, 
3 times smeared); see text. 



picks up a vacuum expectation value due to having the same quantum numbers 
as the vacuum. The standard procedure is to measure these vacuum expectation 
values and subtract them from the correlation matrix elements as 

^^^■W ^jfY.^O,^^)0,{T + t)) - (0.)(0,), (2.30) 

* T 

however during the analysis of our measurements, it turns out to be better to 
treat the vacuum on the same footing as the other states in the vacuum channel, 
i.e. as a state having mass 0. We thus just cut out the vacuum state obtained 
from solving the initial generalised eigenvalue problem (see Appendix ^ and 
perform the fit on the remaining states. 
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2.2 


3.0 


2.7 


#15 


6.0 




4.8 














2.2 




#16 


5.1 


1.4 


2.2 


















#17 


5.1 




2.6 




3.1 














#18 


5.9 




4.7 


2.5 


6.1 


4.2 


3.4 


4.0 





2.1 


0.7 


#19 


6.2 




5.1 




3.1 











2.3 




#20 


6.0 




5.0 


2.5 


6.1 


4.2 


3.4 


4.0 




2.2 


0.7 


#21 


5.9 




4.7 




5.9 






3.9 




1.6 


0.4 


#22 


6.1 


3.0 


5.1 


2.6 


6.1 


4.2 


3.5 


4.0 


2.3 


3.0 


2.2 



Tabic 2.7: "Lifetimes" of glueball correlators at /3 = 3.5 (third smearing scheme, 
9 times smeared); see text. 
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Why Anisotropic Lattices? 

Due to the rather high mass of glueball states and large vacuura fluctuations of 
the operators the signal to noise ratio decreases very rapidly as the temporal 
separation of the source and the sink is increased. This demands a small lattice 
spacing such that the signal can be traced over several slices increasing the 
computational cost drastically because the physical volume still has to be kept 
large enough (La > 1.2 fm) in order to avoid severe finite-volume effects. 

The common way to ensure good resolution of the signal on a lattice of 
moderate size is using anisotropic lattices. 

The history of glueball measurements on the lattice employing anisotropic 
gauge actions is not very long, despite this the most reliable results from pure 
gluodynamics stem from anisotropic simulations. The first study, using a tree- 
level and tadpole improved action has been performed by Morningstar and 
Peardon in 1996 ]123|] . This was followed by additional studies of the same 
group p^,pj|, p^ , p^ . In 1997, the Kentucky Glueball Collaboration stud- 



ied glueball matrix elements on anisotropic lattices [125|. A recent study of the 
glueball spectrum, usi ng t he same action as Morningstar and Peardon, has been 
performed by Liu ||2C , 12£ ] . 



Chapter 3 



Recapitulation: Properties 
of the Isotropic Action 



As the anisotropic perfect action is based on the isotropic FP action paramet- 
rised using "fat" hnks it is useful to recapitulate briefly the properties and the 
results of measurements of physical quantities with the isotropic action . 

The measurements that have been performed using the isotropic action com- 
prise the critical temperature Tc of the deconfining phase transition, the static 
quark-antiquark potential as well as the lower lying part of the glueball spec- 
trum. These quantities and techniques used to measure them are presented in 
detail in Chapter ||. 



This Chapter is organised as follows: In Section 3.1 the action is presented, 



Section 3.2 contains results about the deconfining phase transition, in Section 
3.3 the results about the static qq potential are summarised, in Section 3.4 we 



report about scaling tests of the FP action, and finally, in Section 3.5 we present 
the results about glueball spectroscopy. For additional information, such as run 
parameters and more detailed results, consult 



3.1 The Isotropic Action 

The isotropic FP action has been constructed using the block transformation 



introduced in ||60[ (see Section 12) which has been optimised for a short in- 
teraction range of the FP action as well as for improved rotational invariance 
compared to the standard Swendsen-blocking which uses long staples. 

The action has been parametrised using the parametrisation presented in 



Section 4.3.1 , where the non-linear parameters ry, Ci have been chosen to be 
constant (no x^(n)-dependence), i = 1, 2, 3 and the linear parameters pki in 
the mixed polynomial of the standard and smeared staples u, w are non-zero 



for < k + I < 4:. The values of the parameters are given in Appendix DA 



3.2 The Deconfining Phase Transition 

The properties of the deconfining phase transition in pure Yang-Mills theory are 



described in Section 2.2. On the isotropic lattice we have measured the critical 
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couplings (3c corresponding to temporal extensions Nt = 2,3,4 of the lattice. 



3.2.1 Details of the Simulations Performed 

To determine the susceptibility X; we measure, in the equilibrated system, the 
Polyakov loops averaged over the lattice 



(3,1) 



where P{x) is defined in eq. as well as the action values of the corresponding 
configurations. The measurement is performed after each sweep of updates as 
it is computationally inexpensive. Both values are stored for later use in the 



spectral-density reweighting procedure, described in Section 2.2.1 



We have performed a large number of simulations on lattices with temporal 
extensions A't = 2, 3, 4 at 3 to 6 different (3 values near the estimated critical 
couphngs l3c{Nt). Additionally, various spatial extensions Ns/Nt — 2. 5... 5 
have been explored in order to be able to examine the finite size scaling of Pc- 
The configurations have been generated by alternating Metropolis and pseudo- 
over-relaxation steps (see Sections |8.1.2 ). 

The critical coupling, i.e. the location of the peak of the susceptibility, is 
determined employing the spectral density reweighting method which allows to 
calculate observables away from the actual (3 values at which the simulations 
have been performed 



_ 37|,|127|,|12§] (see Section |2. 2. l| ). 
To get the value of (3c for infinite spatial volume we use the finite-size scaling 
law for first order phase transitions: 



(^ciNt,Ns)=pciNt,^)--h 



(3.2) 



where h denotes a universal quantity in principle independent of Nt 



3.2.2 Results 



In Table we display the val ues of (3c{Vs) together with the infinite volume 
extrapolations according to Eq. 3_^. Studying the deconfining phase transition 
provides accurate information about the scale at couplings /3c. This has been 



used in scaling tests (see Section 3^) where the critical temperature cx 1/a 
is compared to the hadronic scale tq cx a whose determination for the isotropic 
action is described in the following Section. 



3.3 The Static Quark- Antiquark Potential 

The static quark-antiquark potential as described in Section ^.3| is used to set 
the scale of our simulations (determining the lattice spacing a) as well as to 
examine the scaling properties of the FP action, i.e. measuring lattice artifacts. 

3.3.1 Details of the Simulation 

Simulations with the isotropic FP action including measurements of the static 
quark-antiquark potential have been performed at six different values of three 
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Ns 


(ic{Nt = 2) 


f}c[Nt = 3) 


l3c{Nt = 4) 


6 


2.3552(24) 






8 


2.3585(12) 


2.6826(23) 




10 


2.3593(7) 


2.6816(12) 


2.9119(31) 


12 




2.6803(10) 


2.9173(20) 


14 






2.9222(20) 


oo 


2.3606(13) 


2.6796(18) 


2.9273(35) 


h 


0.14(9) 


-0.05(7) 


0.25(9) 



Table 3.1: Results of the critical couplings /3c from the peak location of the 
Polyakov loop susceptibility and the corresponding infinite volume limit ob- 
tained according to eq. 3.2. The finite size scaling constant h is also given. 



of them corresponding to the critical couplings indicated in Section 3.2. The 
updates have been performed using alternating Metropolis and pseudo-over- 
relaxation sweeps. Based on earlier observations in [103,104, 129[ the spatial 
extent of the lattice has always been chosen to be at least (1.5 fm)'^. 

The measurements have been performed as described in Section 2.3.5, for 
measuring the potential we have chosen = 0.2 in the APE smearing of the 
operators (see Section 2.1.2) and used smearing levels n = 0, 1, 2, 3, 4. 



3.3.2 Results 

It turns out that on coarse lattices with Os > 0.15 fm systematic ambiguities 
occur determining on different fitting ranges or using different (reasonably 



chosen) values of c, see Section 2.3.4. The main reason for this is the difficulty 
emerging to define the derivative dV/dR having at hand discrete values of R. 
Estimates for these systematic ambiguities are displayed together with the re- 



sults for ro/a in Table B.2, A possible solution to this problem is including 
off-axis separations of the quark and the anti-quark which makes the measure- 
ments rather expensive. However, violations of rotational symmetry show up 
in this case and allow the estimation of systematic errors in ro if the result is 
changed beyond statistical uncertainty. 

The full static qq potentials aV{R) measured on lattices with different /3, 
i.e. different lattice spacings a, may be compared directly if V{R) is measured 
in terms of tq and the unphysical constant roVb is eliminated by subtracting 



e.g. roy(ro) for each lattice spacing. This yields Figure 3.1 and it turns out that 
all the measurements agree excellently. Also, the result of |130| , represented by 
the dotted line, is hardly distinguishable from our measurements. 



3.3.3 The Interpolating Formula 

The information about the scale obtained in the potential measurements may 
be nicely collected in an (ad-hoc) interpolating formula that connects the bare 
coupling 13 to the lattice spacing (in units of the hadronic scale rg) cl/tq. Our 
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fi 


Nr 


ro/a 


3.400 




4.833(39)(;22) 


3.150 




3.717(23)(+J^) 


2.927 


4 


2.969(14) J 


2.860 




2.740(10)(+3^) 


2.680 


3 


2.237(7)(;jJ) 


2.361 


2 


1.500(5)(1?^) 



Table 3.2: The hadronic scale rp/a determined from local fits to the potential. 
The first error denotes the statistical error and the second is the estimate of 
systematic ambiguities (see text). 
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aP=3.150 








o |3=2.927 








□ P=2.860 


/ 






+ P=2.680 








X p=2.361 









1 2 3 4 5 



Figure 3.1: The static quark- antiquark potential V{R) expressed in terms of 
the hadronic scale rg. To match the curves aX r — tq the unphysical constant 
raV{ra) is subtracted for eac h val ue of (3. The dashed line is a fit of all the data 
to the Cornell potential, eq. ^.24 , the dotted line stems from rcf. [130|. 
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fit (analogous to the ones in |104 , 105|) yields 



ln(a/ro) = -1.1622(24) - 1.0848(95)(/3 - 3) + 0.156(17)(/3 - 3)^ (3.3) 
which describes very well the data points in the range 2.361 < /3 < 3.4. 

3.4 Scaling Properties 

The parametrised FP action is assumed to reduce lattice artifacts consider- 
ably. In order to check this assumption one may study the scaling properties of 
renormalisation group invariant quantities, measuring them on lattices of dif- 
ferent lattice spacings. Dimensionless ratios or products of physical quantities 
which are RG invariant are constant if there are no scaling violations (due to 
lattice artifacts), it is therefore very convenient to study such observables. 

Having at hand the physical quantities tq, the hadronic scale, and a, the 
string tension, from the potential measurements as well as the critical temper- 
ature Tc from the study of the deconfining phase transition, we may build up 
the dimensionless combinations Tc/y^, tqTc and r^^/a. 

In this section we present the results for the FP action and compare them 
to results obtained for the Wilson action and different other improved actions 
as far as they are available. 

3.4.1 Te/v^ 

The results for Tc/i/a, the critical temperature of the deconfining phase tran- 
sition in terms of the string tension, measured on lattices of different [3 are 



compared to the continuum extrapolations for other gauge actions in Table 3.3 



The Wilson action data stem from [131| where results from |132] have been used. 



From the same reference we state the value for the 1x2 tree-level improved ac- 
tion. "DBW2" stands for "doubly blocked from Wils on in two coupling space" , 
the results of this action stem from |l33|. Iwasaki |l34|] has performed mea- 



surements with the RG improved (Iwasaki) action [135|. Finally, the results 
obtained using a tree level and tadpole improved action are from [|l36| . The 
continuum extrapolations for all these different actions stem from a systematic 



reanalysis |137] of the data 



These values, together with those for finite lattice spacings, are displayed in 



Figure 3.2. Our measurements are all compatible with the continuum extrapo- 
lation of the Wilson data and we observe no scaling violations of the FP action 
within statistical errors even on very coarse lattices corresponding to Nt—2. 
The continuum results obtained using the RG improved action as well as the 
tree-level and tadpole improved action lying considerably higher (about 3a for 
the latter) may be explained by ambiguities occurring with the measurement of 
the string tension a. Depending on the separations r of the static quarks taken 



into account performing the fit to the Cornell potential, eq. 2.24, the obtained 
value of the string tension may vary considerably (together with variations of 
the Coulomb parameter a), such that determinations by different groups may 
differ systematically if the methods do not coincide. Especially, Bliss et al. |136{ , 
having studied the treelevel and tadpole improved action, determined the string 
tension a using results coming from torelon measurements and employing the 



string formula, eq. 2.1 
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action 


/3 




FP action 


2.927 


0.624(7) 




2.680 


0.622(8) 




2.361 


0.628(11) 



Wilson |31| 



1x2 tree-level improved [131 

DBW2 ITa 



RG improved [134| 



tree-level/tadpole improved [136 



oo 

CO 

oo 
oo 
oo 



0.630(5) 

0.634(8) 

0.627(12) 

0.651(12) 

0.659(8) 



Table 3.3: Results of the deconfining temperature in units of the string tension 
obtained with the FP action and continuum values for different other actions. 



0.70 



'q 0.65 



0.60 



0.00 0.05 0.10 



O RG improved action 

□ DBW2 action 

1x2 tree level improved 

<Wilson action 

■ FP action 

X tree level/tadpole improved 



1/N 



0^15 



0.20 0.25 0.30 



Figure 3.2: The dimensionless quantity Tc/^/ct vs. for different actions. 
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3.4.2 roTc 



Whereas it was easy to find data for the quantity Tc/^/a presented in the last 
section, precise determinations of tq /a are missing in the hterature except for 
the Wilson action 104 , 105 1 and thus we are not able to compare our results 
for tqTc to a large number of other actions. Our results compared to those 
obtained with the Wilson action are displayed in Table The Wilson results 
for the critical couplings (3c{Nt), Nt = 4, 6, 8, 12 have been taken from |131| and 
the corresponding values of r^/a have been calculated using the interpolating 
formula from |105|. The continuum value has been extrapolated using a linear 



fit in the leading correction term on Nt = 4, 6, 8. Figure 3.3 displays the 



same quantities; including the systematic ambiguities on coarse lattices (dashed 
lines), see Section 3.3.2. 



A^T Wilson action FP action 



2 




0.750(3) 


3 




0.746(3) 


4 


0.719(2) 


0.742(4) 


6 


0.739(3) 




8 


0.745(3) 




12 


0.746(4) 




00 


0.750(5) 





Table 3.4: Results for the critical temperature in terms o f the hadronic scale, 
r^Tc, from measurements with the Wilson action |l05, 131|, and the FP action. 



Compared to the Wilson action, which shows scaling violation of about 4 % 
at TVt = 4 the FP action performs much better in this scaling test exhibiting 
scaling within 1 % even on very coarse lattices corresponding to A^'t = 3 and 2. 
However this statement is softened due to the large systematic uncertainties at 
larger lattice spacings. 



3.4.3 rov^ 



The dimensionless product r^y/a contains two (physical) quantities obtained in 
the fit of the static quark-antiquark potential. However its scaling properties are 
not trivial as ro probes the intermediate range and a the long range properties 
of the potential. We are t hus using the values of ro/a obtained from local 
fits to the potential (Table 3^) , whereas the valu es o f cra^ are determined 
from global fits including large values of R. Table 3.5 collects the resulting 
values of r^^fa together with results from the Wilson action calculated using 



the interpolating formula for rp/a from |105| and the parametrisation for the 
string tension from J104|. The continuum extrapolation of the Wilso n data has 
been performed in |137[ . The same data is displayed in Figure 3.4. The error 
bars are purely statistical and dominated by the errors of the string tension — 
the systematic ambiguities in ro may thus be neglected in this case. 

Again, the Wilson action shows scaling violations of about 4 % for A''^ = 4 
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0.77 



0.76 



0.75 



0.74 



0.73 



0.72 



O Wilson action 
■ FP action 



0.71 



0.00 0.05 



0.10 0.15 
1/N' 



0.20 0.25 0.30 



Figure 3.3: The dimcnsionle ss quanti ty roTc vs. l/N^ for the FP action com- 
pared to the Wilson action 1 105, 131 1, the dashed error bars indicate the sys- 
tematic error from ambiguities in calculating the force F{R) — V'{R) on coarse 
lattices. 



Wilson action 




FP action 






(3 


roV^ 


5.6925 


1.148(12) 


2.361 


1.194(21) 


5.8941 


1.170(19) 


2.680 


1.196(15) 


6.0624 


1.183(13) 


2.860 


1.190(23) 


6.3380 


1.185(11) 


2.927 


1.191(12) 






3.150 


1.185(16) 






3.400 


1.198(12) 


oo 


1.197(11) 


oo 


1.193(10) 



Table 3.5: for the Wilson pM , 



105 



137 1 and the FP action. 
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1.25 



1.20 



O 1 



1.15 



O 



O Wilson action 
■ FP action 



1.10 



0.00 0.10 0.20 0.30 



0.40 



0.50 



Figure 3.4: Scaling behaviour of rQ^/a for the Wilson action [104, 105, 137| 
(empty circles) and the FP action (filled squares). 
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whereas the FP action exhibits no artifacts even on lattices as coarse as f3 — 
2.361 (corresponding to criticality for Nt — 2). 



3.5 Glueball Spectroscopy 

The physical nature of glueballs and how they can be studied on the lattice is 



described in Section 2.4. Numerical simulations using the methods described 
there have been performed using the isotropic FP action and will be described 
in this Section. 



3.5.1 Details of the Simulation 



The simulations have been performed on three different lattice spacings in the 
range 0.10 fm < a < 0.18 fm and volumes between (1.4 fm)'^ and (1.8 fm)'^. The 
updates have been performed using alternating Metropolis and pseudo-over- 
relaxation sweeps (se e Sec tions 8.1.1 , 8.1.2 ). For the listing of all the simulation 
parameters see Table |3.6|. 





lattice 


ro/a 


a [fm] 


# sweeps / measurements 


2.86 


10^ 


2.740(10) 


0.18 


32000 / 8000 


3.15 


124 


3.717(23) 


0.13 


101000 / 10100 


3.40 


144 


4.833(39) 


0.10 


86520 / 14420 



Table 3.6: Run parameters for the glueball measurements using the isotropic 
FP action. 



In the first simulation, performed on a coarse lattice with /3 = 2.86, only the 
mass of the lowest-lying A^^ ch ann el has been determined measuring the loop 
shapes 2, 4, 7, 10, 18 (see Figure 2/7) which had been shown to be measured well 
and having a large overlap with the state in a short pilot run. For the other two 
simulations all 22 Wilson loop shapes up to length 8 are measured and projected 
to all 20 irreducible glueball representations. To kill ultraviolet fluctuations and 
to enhance the overlap with the glueball states the measurements are performed 
on five smearing levels SnU, n = 2, 4, . . . ,10 with smearing parameter As = 0.2 
(see Section 2.1.2| ). 

It turns out that there are operators which are measured notoriously bad, 
i.e. that have a small signal/noise ratio. As there are representations with con- 
tributions of as many as 145 operators (the representation on five smearing 
levels) this may destabilise the fitting procedure described in Appendix That 
is why we remove badly measured operators right from the start of the analysis. 

For the scalar glueball, the vacuum expectation values are not subtracted 
from the operators, but the vacuum state is treated like the other states in the 
vacuum channel, see Section 2.4.2. 



3.5.2 Results 

The full (and rather lengthy) lists of results are contained in ||^, we only display 
the masses of the lowest-lying and best measured states in units of the hadronic 
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scale ro compared to the results of other groups in Tables 3.7, The results 



denoted by UKQCD, Teper and GFll have been obtained using the Wilson 
action whereas the results of Morningstar and Peardon (M&P) and Liu have 
been measured using a tree-level/tadpole improved anisotropic action. 



Collab. 



romo++ 



rom2++ 



year 



UKQCD 1 107] 



Teper [|137 



GFll |138| 



M&P [18 1 
Liu I20II 



4.05(16) 
4.35(11) 
4.33(10) 
4.21(15) 
4.23(22) 



5.84(18) 

6.18(21) 

6.04(18) 

5.85(8) 

5.85(23) 



1993 
1998 
1999 
1999 
2000 



FP action 



4.12(21) 



[5.96(24)] 



2000 



Table 3.7: Comparison of the two lowest glueball masses in units of rg. Our 
2++ value is not extrapolated to the continuum but is the mass obtained at a 
lattice spacing a = 0.10 fm. 



Collab. roTOQ-+ rom2-+ rQmi+- year 



Teper 
M&P 


13 
18 


7 


5.94(68) 
6.33(13) 


8.42(78) 
7.55(11) 


7.84(62) 
7.18(11) 


1998 
1999 


FP action 


[6.74(42)] 


[8.00(35)] 


[7.93(78)] 


2000 



Table 3.8: Comparison of glueball masses in units of tq. Values in brackets 
denote masses obtained at a lattice spacing a = 0.10 fm and are not extrapolated 
to the continuum. 



The results for the sc alar (0^^) glueball are displayed together with the ones 
of other groups in Figure [s^ , the results for the 2++ glueball which is composed 
of the two representations and T^^ in Figure 3.6. 
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X Wilson action 

C treelevel/tadpole improved 

anisotropic action 
• FP action 



^ 0.05 0.1 0.15 0.2 

OA)' 

Figure 3.5: Mass estimates for the 0++ glueball. Results from simulations of 

the Wilson action (crosses) and a tree Icvcl/tadpolc improved anisotropic action 
(empty circles) are shown together with the results obtained with the FP action 
(filled circles). 
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7 



6 



E 




X Wilson action 

O treelevel/tadpole improved 

anisotropic action 
• FP action 



4 I . . . . , . . . . , . . . . , . . . . 1 

0.05 0.1 0.15 0.2 

Figure 3.6: Mass estimates for the 2+"'" glueball. Results from simulations 
of the Wilson action (crosses) and a tree level/tadpole improved anisotropic 
action (empty symbols) are shown together with the results obtained with the 
FP acticm (filled symbols). Squares and circles denote the £■++ and T2"^ mass 
estimates, respectively. 



Chapter 4 



The Construction of Perfect 
Anisotropic Actions 

In this Chapter we describe our concept of generating anisotropic perfect actions 
without having to repeat the whole lengthy procedure described in Section 1.4, 
making use of the parametrisation of the isotropic FP action [|l]-D . In Section 



4.1 we describe two methods one might think of, in the following we concentrate 



on the one of the two ways which is presented more deeply in Section 4.2. The 



parametrisation ansatz, for the isotropic as well as for the anisotropic case is 



finally presented in Section 4.3 



4.1 Methods 



The method for constructing FP actions (see Section 1.3.2 ) may be straight- 



forwardly generalised to anisotropic lattices in t wo different ways, illustrated, 



together with the isotropic method, in Figure |4.1| : Firstly, one might repeat the 
procedure used to generate the isotropic FP action (see Section 1.4), starting 
with an anisotropic action A^^ which behaves well on very small fluctuations 
and doing the same cascade process involving blocking steps and parametrisa- 
tions as before. Of course, the parametrisation has to be adapted such that it 
may account for the space/time asymmetry. 

Secondly, one may start with an isotropic parametrised FP action on rather 
large fluctuations (e.g. with the action Af^'^, see Appendix ^ and perform a 
small number of purely spatial blocking steps. If a scale 2 block transformation 
is used, the resulting actions will have anisotropics ^ = 2, 4, 8, .... This second 
approach may be performed much faster than the first one because the number 
of minimisation/fitting steps is smaller and one does not have to construct and 
check a new (anisotropic) starting action Aq\ It is because of this we concentrate 
in this work on the spatial blocking approach described in the next section. 



4.2 Spatial Blocking 

To perform the blocking only spatially, doubling the lattice spacing in spatial 
direction while leaving the temporal lattice spacing unchanged, we modify the 
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* ♦ t 



* ♦ * 



Figure 4.1: Different methods of creating a FP action, left: isotropic action; cen- 
ter: anisotropic action, using an anisotropic starting action; right: anisotropic 
action, performing purely spatial blocking steps in the end. 
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blocking kernel T{U, V) of the symmetric RGT 3 transformation ||6^], appearing 
in cq. 



T{U, V)=Y^ (AA-(nB) - ^Re Tr[V,{nB)QlinB)]) , (4.1) 

nB,fJ. 

with the normalisation 

Af^ins) = ^max^^ {^Re Tr[WQl{ns)]} (4.2) 

and the blocked link 

Q^(ns) = W^i2nB)Wp{2nB + A), (4.3) 

where Wp(n) denotes a smeared (fuzzy) fine link, which is constructed as follows. 
In addition to the simple staples we build as well diagonal staples by first going 
in the planar or spatial diagonal directions orthogonal to the direction fi of 
the final smeared link, then doing a step in direction fi and finally returning 
along the corresponding diagonal to n + jl. To be specific, let us first build 
matrices ) connecting the sites n and n' where n' is a neighbouring 

site orthogonal to the smeared link W^(n) to be built: 

iy(o)(n,n) = 1, (4.4a) 
W'-^Hn,n + i') = U^{n) , (4.4b) 

W^^Hn,n + {^ + p) = ^iU4n)Upin + i^) + Upin)U,in + p)) , (4.4c) 

W'^^\n,n + i) + p + \) ^ i {Ut,{n)Up{n + i')Ux{n + O + p) + perms.) . 

(4.4d) 

Here i/, p and A go over all (positive and negative) directions different from 



p and from each other. In eqs. 1. 4c, 4.4c] the sum is taken over all shortest 



paths leading to the endpoint n' of the corresponding diagonal. The fuzzy link 



operator Wp,{n) is then constructed by a modified Swendsen smearing |139-143 



W^.{n) = 5] ^c™VF('")(n,n')C/^(n')VF(")(n' + /i,n + /i). (4.5) 

rn—O n' 

The coefficients Cm are free parameters subject to the constraint: 

Co + 6ci + 12c2 + 8c3 = 1. (4.6) 

This normalisation condition ensures that for a trivial field configuration, Wp{n) 
is equal to the unit matrix. These parameters Cm have been optimised for a 
short ranged quadratic FP action j6^] and take on the values ci=0.07, C2=0.016, 
C3=0.008 and thus co==0.324. 

To modify the isotropic definition, in order to end up with a purely spatial 
blocking, we set 

Qi{nB) = Wi{2nB) (4.7) 

thus doing only a smearing and no blocking in temporal direction. (In the 
case of spatial blocking the expression "2nB" stands for (2n]j, 2n^, 2n^, n^).) 
Additionally, we have the freedom to choose different values of k for spatial (k^) 
and temporal {nt) links (see also Appendix^. 
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4.3 The Parametrisation 

As already stated, the FP-type actions, described in principle by an infinite num- 
ber of couplings, have to be parametrised in order to be useful for numerical 
simulations. We have shown that a parametrisation of the isotropic FP gauge 
action including APE-like smearing behaves much better and is much more 
flexible compared to common parametrisations using traces of closed loops with 
comparable computational cost . Thus the starting point for the parametri- 
sation of the anisotropic perfect gauge action is our isotropic parametrisation, 
presented in Section 4.3. l[ . We extend this ansatz to suit anisotropic actions in 



Section 1.3.2. The different extensions will be compared in Section 4.3.4 



4.3.1 The Isotropic Parametrisation 

The isotropic action has been parametrised using mixed polynomials of traces 
of simple loops (plaquettes) built from single gauge links as well as from (APE- 
like) smeared links. In this Section, the parametrisation will be described in 
detail. 

We introduce the notation S'^'^'' (n) for the sum of two staples of gauge links 
connecting two lattice sites in direction fi lying in the /ii/-plane: 

Sl^\n) = U,{n)U^{n + v)Ul{n + A) + Ul{n ~ v)U^{7i - 0)U,{n - + fi). 

(4.8) 

Our parametrisation may depend on local fluctuations measured by x^(ji) 
which is defined like 

x^in)=ReTT{Ql{n)Ulin)), (4.9) 
with the symmetrically smeared link 

This parameter is negative, —4.5 < x^{n) < 0, it vanishes for trivial gauge 
configurations: x^{n) = 0. 

To build a plaquette in a plane from smeared links it is convenient to 
introduce asymmetrically smeared links. First define]^ 



(^) 



Sl'^ + 77(a;,)5M j - (^1 + ^vix,)^ U,. (4.11) 

Out of these sums of matrices connecting two neighboring points n, n + jl, we 
build the asymmetrically smeared links 

l^M = u,, + ci(a;^)Q|r) + C2(x,0QMt/tQM + . . . , (4.12) 

where r]{x), Ci{x) are polynomials with free coefficients (determined later by the 
fit to the FP action): 

r^(x) = r]^°^ + r]''^'^x + rj^^'^x^ + ■■■ , (4.13) 
c,{x) = cf^ + cf ^ X -I- cf^a;2 -)-... . (4.14) 



^The argument n is suppressed in the following. 
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Of course, these asymmetrically smeared links, built out of a large number 
of paths connecting the neighbouring lattice sites are no longer elements of the 
SU(3) gauge group. They might be projected back to SU(3), however this task 
increases the computational cost in actual numerical simulations, and addition- 
ally our studies have shown that projection reduces the degrees of freedom in 
defining the action, we are thus using the smeared links wl^^ as they are. 

We are now able to build a smeared plaquette like 

w,,,^ReTi{l~Wll), (4.15) 
as well as the ordinary plaquette 

u^, = Re Tr(l - Uf,), (4.16) 



where 



and 



Wll{n) = Wl^Hn)Wi'^\n + + v)Wi^^'^\n) (4.17) 



Ufi{n) = U^{n)UAn + (i)Ul{n + v)Ul{n). (4.18) 
Finally, the action is built out of these plaquette variables as 

-^[f^] = ]^ E E Epfc'^^^''(")'^^^''(")'' (4.19) 

^ n fj,<u k,l 

where the coefficients pki are again free parameters defined by the fit to the FP 
action. 



For the actual fitting procedure (see Sections 1.4.4 , 1.4.5 ) it is very important 
to note that the parameters pki enter linearly, whereas the parameters Ci as well 
as rj enter non-linearly. 

4.3.2 The Anisotropic Extensions 

The isotropic parametrisation presented in the last section does not distinguish 
between spatial and temporal directions on the lattice. It is thus not able to 
describe anisotropic actions at all. However there are some straightforward ways 
of extending it to parametrise anisotropic perfect actions. 



Firstly, the coefficients pui in eq. 4.19 may be chosen differently depending 
on the orientation of the plaquette ^jlv, i.e., for G {12,13,23} (spatial 
plaquettes) and p^™ for ^iv G {14,24,34} (temporal plaquettes). 



Secondly, the parameter ry entering in eq. 4.11 describing the asymmetry 



between different staples contributing to a smeared link wji^^ depending on their 
orientation may be generalised to distinguish between spatial and temporal links 
contributing to spatial and temporal smeared plaquettes W^^: 

Q?' = ^(E^f'+^i^P^+^35f')-(i + ^^i + iry3)C/„ (4.20) 



l(4) 



= i( E s^^'^ + msl"^) - (1 + lm)u„ (4.21) 



^(E si^^+v,si'^)-{l + ^v,)U„ (4.22) 

(4.23) 
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*■ 



vi m m m 

Figure 4.2: Asymmetry in the construction of the smeared hnks matrices Q^'^''. 
From left to right: spatial staple contributing to a spatial smeared link in a 
spatial plaquette (rji), temporal staple contributing to a temporal smeared link 

(772) 7 temporal staple contributing to a spatial smeared link in a spatial plaquette 

(773) , temporal staple contributing to a spatial smeared link (774). 




Cil Ci2 CiS 



Figure 4.3: Asymmetry in the construction of the smeared links Wj),'^ out of the 
matrices Q^l''^ ■ From left to right: spatial finally contributing to a spatial 
smeared plaquette (cn), spatial finally contributing to a temporal smeared 
plaquette (€,2), temporal Q4 (0^3). 



where the anisotropic parameters rji, . . . ,774 may be again polynomials in the 



local fluctuation x^. These situations are depicted in Figure 4.2 



Finally, also the construction of the smeared links W^i^'^ out of the matri- 
ces Q^i^\ described by the parameters q in eq. 4.12 is generalised such that 
these parameters are chosen differently for temporal links (always contribut- 
ing to temporal plaquettes), spatial links contributing to spatial plaquettes and 
spatial links contributing to temporal plaquettes, introducing qi, Ci2 and ^3 
respectively, see Figure 

It is not a priori clear whether all the three extensions are necessary at the 
same time. Results from fits to different extended parametrisations are given in 



Section 4.3.4 



4.3.3 The Normalisation Condition 

Our anisotropic action shall have the correct normalisation and anisotropy in 
the continuum limit. This can be accomplished demanding the following two 
normalisation conditions to be exactly fulfilled: 



Poi+-Pfo + 2poiCi3+PoiCi3?72 + 2p^lcn773 -|-p^*iCi27?4 = ^, (4.24) 
Pol+Pio + 2po'iCii + 2p^*iCi2 + 2p^lcii77i = ^ (4.25) 
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where all parameters 77, c denote the constant (O"' order) terms in the polyno- 
mials in X. The way how these normalisation conditions are derived is presented 
in Appendix |^. 



4.3.4 Comparing Different Sets of Parameters 

During the construction of the perfect ^ = 2 action (see Chapter ^ we check the 
quality of the fits to different anisotropic extensions of the APE parametrisation. 
As it turns out, it is indispensable to have the full parameter set including four 
different asymmetry parameters 77 as well as distinguishing between different 



plaquette orientations (cii, Ci2, ca) in building up W^^. Table 4.1 lists the 
results of fits to the ^ = 2 perfect action as well as measured anisotropics 
using the torelon dispersion relation (see Section |6.2| ). It is clear from these 
results that setting all Cy, j = 1,2,3 equal leads to a large renormalisation of 
the anisotropy; as well, setting rya = 774 = in order to have a positive definite 
transfer matrix decreases the quality of the fit considerably. We thus choose 
to use at the same time anisotropic linear parameters (p^^ and p^™) as well as 
non-linear parameters 771, 772, rya, 774 and cn, Ci2, Ci3 to parametrise our final 
anisotropic perfect gauge actions. 

The proceeding and the results of the parametrisation of anisotropic perfect 
actions are presented in Chapter || for the anisotropy ^ = 2 and Chapter for 
e = 4. 
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0.0238 
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4 


0.0144 


1.912(9) 



Table 4.1: Comparison of the accuracy of the fit on derivatives Xd '^^ 20 con- 
figurations at /3 = 3.5 and the measured renormalised anisotropy (using 
the torelon dispersion relation) at /3 = 3.3 for different choices of the set of 
non-linear parameters. 



Chapter 5 

Scalar Fields and 
Perturbative Tests 



The methods for creating anisotropic actions presented in Chapter ^ are first 
examined for the case of d = 2 free scalar field theory. In this theory, the renor- 
malisation group transformation and thus the calculation of the FP action may 
be performed analytically. In Section 5.2 we study the quadratic approximation 
to the FP action in d = 4 gauge theory. 

It turns out that the ansatz works in both theories and that the results look 
very similar to what is obtained for the isotropic FP action. 



5.1 The Free Scalar Field in d=2 



The FP action for free scalar fields has been studied earlier |144], a recent study 
also involves anisotropic lattices and massive fields [145|. 



5.1.1 Anisotropic Blocking out of the Continuum 

A simple way to get a FP action for the free scalar field in d = 2 is to block it 
out of the continuum choosing different lattice spacings Os and at in temporal 
and spatial direction, respectively. The FP propagator is derived in Appendix 
0. The result is 

with a 27r-periodic function Tl{q) defining the blocking. The momentum q = 
[qi, q2) is expressed in units of the lattice spacing a and runs in the first Brillouin 
zone ] — 7r/a, 7r/a]. The parameter k is due to the generalisation of the blocking 
from S functions to a Gaussian (k — > oo corresponds to the S function), i.e. it is 
a sort of a smearing parameter (with extensive smearing for small values of n). 
For isotropic lattices it can be shown analytically that k = 6 yields ultra-locality 
for d = I. The blocking function n(fc) is a dimensionless quantity, while the 
propagator and k carry the dimension m^^ and respectively. It is purely a 
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matter of convention whether we express (or regulate) these quantities in terms 
of as or at, yielding different ^-factors at different places: 

where Ks and the propagator are now expressed in terms of and 

if they are expressed in terms of at. The values KsC^ — Kfol = & are optimal 
for locality in spatial and temporal direction, respectively, and we obtain ultra- 
locality for the Id projection onto the corresponding axis. As, k, ^ m? we get 

The corresponding FP action is defined as: 

S^'' ^\Y.p{r)Mn+r (5.4) 

n,r 

on the scalar field 0r, where 

p{q)^D-\q). (5.5) 

If we require the action to go over to the form of classical field theory for 
small momenta q, the coefficient of the O(g^) term in the action should be 1: 

oo oc 

-5(9)= E E P(n)e-'''" ^ forq«l. (5.6) 

ni — — oo n2 — — oo 

Expanding exp(— igyi) and comparing the 0{q^) terms we obtain the sum rules: 

^ nlp{ns,nt) = -j, (5.7) 
J2ntP{ns,nt) = -2e, (5.8) 

if we express our quantities symmetrically in a^at. 

Requiring the spectrum to be massless we obtain another (trivial) sum rule: 

Y,Pins,nt) = 0. (5.9) 

(5.10) 

Calculating the couplings in position space p(r) we can either choose k = 6/^ 
which yields optimised spatial locality and ultra-locality at the projection to 
the temporal axis, or k = 6^ for optimal temporal locality and ultra-locality 
after projecting to the spatial axis. The resulting magnitudes of the couplings, 
|p(rs,0)| and |p(0,rt)| are displayed in Figures 5.1, 5.3, 5.5 for anisotropics ^ = 2, 
4, 10, where the separations and rt are given in lattice units. Note that the 
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physical separation in temporal direction is a factor ^ smaller compared to the 
spatial one. 

By truncating the action, we introduce 0{af) and 0{af) artifacts. As for 
our anisotropic actions as ^ at we expect the temporal artifacts to be much 
smaller than the spatial ones and therefore truncate the action to nearest- 
neighbour couplings in temporal direction ( "nnt" ) which is only possible without 
completely distorting the sum rules, the spectrum etc. if the action has been 
optimised for temporal localit y choosin g k = 6^. This truncation scheme yields, 
employing the sum rules eqs. 



p{ns,l) ^ p{ns,-l) = ^^njp^^{ns,nt), 



(5.11) 



p(n„0) = ^/P(n„nO-5]ntVK,"0- (5-12) 



To estimate the physical quality of the truncated perfect actions we display 
the deviations of the exact spectrum E{p) ~ p/^ in the first Brillouin zone for 
different truncations and different values of k in Figures 5^ , 5.4, 5.6. 



Comparing the spectra of the symmetrically truncated actions (5 x 5) it turns 
out that indeed the temporally optimised couplings {k — 6^) yield much better 
spectra than the spatially optimised couplings — for ^ = 10 the spectrum is even 
completely off in the latter case. As well, truncating these actions to nearest- 
neighbour in time does only mildly decrease the quality of the spectra. The 
larger the anisotropy gets, the more important is good locality of the temporal 
couplings. 



5.1.2 Spatial Blocking 

For the {d — 2) scalar field, we may perform purely spatial blockings (cf. Section 
4.2) analytically and check whether (in this model) we end up with the perfect 
anisotropic action, eq. We start from the isotropic FP propagator which we 
know analytically for a blocking that averages the continuum field isotropically 
around lattice sites (without overlapping): 



sin'(g,./2) 



n 



1 

3k ' 



(5.13) 



and perform a single spatial {p 
spacing by a factor of two: 



1) RG step, increasing the spatial lattice 



l',=0 



(note the normalisation of the k factor). Inserting expression eq. 5.13 and 
performing the substitution 171/2 + ttI[ + 27rZi — > qi/2 + nli we get 



<^(9i,92) 



1 



sin2(f )sin2(f ) 



1 

2.^ (f + nhY + (92 + 27r/2)2 (f + 7r/i)2(f + nl^y 



I — — 00 

1 

3k ' 



(5.15) 
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Figure 5.1: Couplings of the ^ = 2 free scalar field in d = 2, optimised for spatial 
locality {k = 3, circles: on the spatial axis, squares: on the temporal axis) and 
temporal locality {k = 12, pluses: on the spatial axis, crosses: on the temporal 
axis) . 
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Figure 5.2: Deviations from the exact dispersion relation of the ^ = 2 perfect 
action for the d = 2 free scalar field. Circles: k ~ 3, 5 x 5 couplings, squares: 
K = 3, 5x nnt couplings, pluses: k = 12, 5 x 5 couplings, crosses: k = 12, 
5x nnt couplings. 
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Figure 5.3: Couplings of the ^ = 4 free scalar field in d = 2, optimised for 
spatial locality {k = 1.5, circles: on the spatial axis, squares: on the temporal 
axis) and temporal locality {k = 24, pluses: on the spatial axis, crosses: on the 
temporal axis). 
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Figure 5.4: Deviations from the exact dispersion relation of the ^ = 4 perfect 
action for the d = 2 free scalar field. Circles: k = 1.5, 5x5 couplings, squares: 
K = 1.5, 5x nnt couplings, pluses: k = 24, 5x5 couplings, crosses: k = 24, 
5x nnt couplings. 
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Figure 5.5: Couplings of the ^ = 10 free scalar field in d = 2, optimised for 
spatial locality {k = 0.6, circles: on the spatial axis, squares: on the temporal 
axis) and temporal locality (k = 60, pluses: on the spatial axis, crosses: on the 
temporal axis). 
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Figure 5.6: Deviations from the exact dispersion relation of the ^ = 10 perfect 

action for the d = 2 free scalar field. Circles: k = 0.6, 5x5 couplings, squares: 
K = 0.6, 5x nut couplings, pluses: k = 60, 5x5 couplings, crosses: k = 60, 
5x nnt couplings. 
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which is the anisotropic C = 2 perfect propagator for this model, compare eq. 5.3 



(Performing another single temporal RG step on this propagator gives back the 
initial FP propagator D^^ , which proves that the normalisation of both terms 
has been chosen correctly.) Of course this discussion of a single spatial RG 
step in the d — 2 scalar theory cannot be generalised straightforwardly to gauge 
theory as the RG step cannot be treated analytically there, however these results 
indicate that this approach may be feasible as well for gauge theory. 



5.2 Gauge Fields in the Quadratic Approxima- 
tion 

We will now turn to SU(3) gauge fields in d = 4 dimensions. On small fluctua- 
tions we may approximate the FP gauge action quadratically: 

^^"^ = ^ E P,Ar)TT[A,in + r)AAn)] + OiA') (5.16) 

^ n,r 

neglecting terms of order 3 and higher in the gauge potential. In momentum 
space the action reads 

= E P,Ak)TT[A,i-k)A4k)] + 0{A'). (5.17) 

The calculation of the couplings Pfj.u{r) and the propagator D^^{r) is de- 
scribed in Appendix 



5.2.1 Physical Quantities 

Having at hand the quadratic couplings p^i, and the propagator D^i, we may 
easily calculate physical quantities, as e.g. the spectrum of the action or the tree 
level potential as described in the next sections 



The Spectrum 

The quadratic FP action, eq. |S.16 , is defined on the lattice. Nevertheless 
it describes transverse massless gluons with the exact relativistic spectrum 
E(k) = |fc| for ki — — OO..CX). This is a fundamental difference to perturba- 
tively (e.g. Symanzik) improved actions where the spectrum shows O(a^) or 
0(0*) lattice artifacts. 

The spectrum is determined by looking at the singularities of the propagator 
Dfii,{k): We fix the momentum k and determine the corresponding energy k^. 
This may be done analytically using exactly the same arguments as for the 
isotropic FP action (see ||5^) with the result that the energy behaves as k^ — 
—i\k\, where k ^ q + 2t:1, which is the exact continuum relation realised by a 
tower of poles in D^,^ for any given q, qi e (0, 27r). 



The Perturbative Potential 

Another physical quantity that can be determined using the quadratic couplings 
is the static quark-anti-quark potential up to tree level. 
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At the quadratic level the potential in a finite continuum box of size is 
proportional to the Coulomb potential: 

Vco„t(:E) = ^ Ve^''^-i-, where ki = ^k. (5.18) 

On the lattice, the corresponding relation reads 

Mat(r) e''''^44(fc, = 0), (5.19) 

where I?^iy(fc) is the propagator corresponding to the quadratic part of the lattice 
action. As this quantity is unchanged from the isotropic case (except a global 



factor 1/^) the result is exactly the same as in [|60[. Figure 5/7 summarises the 
results for the anisotropic actions. As V ^ 1/^ the dev iations are accordingly 
smaller for larger values of the anisotropy ^. Figure compares the deviation 
to the continuum result of the perfect action to the deviation present with the 
Wilson action (for the isotropic action; the picture is exactly the same for ^ 7^ 1) 
— note the difference in the scale of the deviations. It is self-evident that the 
perfect action heavily suppresses discretisation errors. The deviations that are 
still present are on one hand due to the RGT transformation not averaging 
exactly spherical, on the other hand due to remaining quantum fluctuations on 
the fine lattice. 
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Figure 5.7: The deviations to the exact tree-level quark-anti-quark potential of 
the quadratic approximation to the perfect anisotropic action for ^ = 1 (circles), 
^ = 2 (squares), ^ = 3 (triangles), <^ = 4 (crosses). 
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Figure 5.8: The deviations to the exact tree-level quark-anti-quark potential of 
the Wilson action (crosses) and the FP action (circles). 



Chapter 6 

The C =2 Perfect Action 



The main topic of this work is the construction and the apphcation of a classi- 
cally perfect gauge action for ^ = 2 anisotropic lattices. These tasks and their 
results are presented in this Chapter. The construction follows the proceeding 
presented in Section 1.4 using the spatial blocking method and the parametri- 
sation described in Chapter ^. 

The anisotropy of the resulting action is measured using the torelon disper- 
sion relation (see Section the (spatial) scale is determined by the static 
quark-antiquark potential (see Section |2.3| ) , where for some values of the cou- 
pling P, off-axis separations of the quarks are considered as well, to be able to 
estimate violations of rotational symmetry by the parametrised action. Another 
way of determining the (temporal) scale is the study of the deconfining phase 
transition (see Section 2.2), measuring the critical couplings /3crit for different 
spatial extensions Nt of the lattice. 

Finally, the glueball spectrum is measured performing simulations at three 
values of (3 using the ^ = 2 action. 

On one hand, the study of the classically perfect action with rather small 
anisotropy ^ = 2 is a feasibility test of the ansatz as it employs only one single 
spatial blocking step and one single parametrisation of the resulting action. On 
the other hand, to determine accurately the continuum limit, one should perform 
simulations on fine lattices as well — using actions with large anisotropics, this 
may be more expensive than necessary, at least when the states to be studied are 
not extraordinarily heavy. The feasibility of additional spatial blocking steps, 
followed by parametrisations is briefly (and incompletely) checked in Chapter 





6.1 Construction 

We construct the ^ = 2 classically perfect anisotropic action starting with the 
parametrised isotropic action performing one spatial blocking step as de- 
scribed in Chapter ^ The quadratic approximation suggests nt — ^^Hs = 
for keeping the block transformation as close to the isotropic case as possible. 
Furthermore, spatial locality should stay at its optimum for the "isotropic" 
value of Ks = 8.8, thus we choose Kt = 35.2. Indeed, it turns out that varying 
Kt away from nt = 35.2, keeping Ks = 8.8 fixed, makes it much more difhcult to 
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parametrise the resulting perfect action such that the renormaHsed anisotropy 
(see Section |6^) stays close to the input value ^ = 2. 

Examining this blocking, we observe that coarse configur ation s obtained in 
a MC run with the ad-hoc anisotropic action (see Appendix pl).3| ) are mapped 
to fine configurations that are close to isotropic (concerning spatial and tem- 
poral plaquctte values and expectation values of the Landau gauge-fixed link 
variables) for input anisotropy ^ad-hoc ~ 3.2. 

We thus generate 20 configurations each at /3ad-hoc = 2.5, 3.0, 3.5, 4.0 using 
the ad-hoc anisotropic action^ with Cad-hoc=3.2. These configurations are spa- 
tially blocked using the isotropic action Af^^ (see Appendix p.2| ) to describe 
the isotropic minimised configuration on the r.h.s. of eq. 1.5. 



T o construct the action we follow the procedure described in Sections |1.4.3 
f .4.5| . Several non-linear fits are performed for the derivatives of 20 configura- 



tions at P — 3.5, using different sets of parameters. Concerning the non-linear 
parameters, we find out that it is necessary to keep all asymmetry parameters 
r]i non-zero. Moreover, it is crucial to distinguish the Ci parameters for different 



plaquette orientations, see Table 4.1 and Section 4.3.4. The parameters Cj are 



chosen to be non-zero for i = 1, 2, 3 as adding additional free C4 parameters does 
not improve significantly. 

Having fixed the non- linear parameters 77 and Ci, we include the derivatives 
of 20 configurations at /3 = 4.0 and the action values of all the configurations 
at /3 = 3.5, 4.0 in the linear fit, keeping the non-linear parameters fixed. The 
resulting values of and the linear behaviour of the actions suggest a linear 
set with max(/c -I- ^)sp = 4, max(A: -I- ^)tm = 3 where the action values are 
included with weights Wact — 0.018 relative to the derivatives (see also Figure 



1.3). Using this value, the good parametrisation of the derivatives obtained in 



the full nonlinear fit is preserved, at the same time the mean error of the action 
value due to the parametrisation is as small as 0.35%. 

It is checked that the values of are still acceptable on sets of configurations 
down to /3 = 3.0, besides that they do not increase significantly for configura- 
tions independent of the ones used in the fitting procedure. The behaviour of 
the linear parameters of the action as a function of the s imp le plaquette values 
u and smeared plaquette values w is displa yed in Figure |6.l| , the parameters of 
the resulting action are listed in Appendix D.4. 



6.2 The Renormalised Anisotropy 

The renormalised anisotropy of the action presented above is measured us- 



ing the torelon dispersion relation as described in Section 2T. The main ad- 
vantage of this approach is obtaining the renormalised anisotropy as well as 
an estimate for the scale a^, at performing computationally rather inexpensive 
measurements. 



^Note, that these values of /3ad-hoci used in this section, may not be directly compared to 
the values of /3perf corresponding to the parametrised classically perfect action, used in the 
following sections; but rather the fluctuations (plaquette values, expectation values of Landau 
gauge-fixed links etc.) of the gauge configurations produced in MC runs using the perfect 
action should lie approximately in the same range as the fluctuations of the initial coarse 
configurations produced with the ad-hoc action. 




Figure 6.1: The behaviour of the linear parameters pki as a function of the 
values of the standard plaquette (u) and the smeared plaquette (w) for the 
parametrisation of the classically perfect ^ = 2 action. 
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6.2.1 Details of the Simulations and Results 



The parameters of the simulations carried out are cohected in Table 6.1 





X LxT 


# sweeps / measurements 


3.00 


8^ X 4 X 20 


72000 / 14400 


3.15" 


8^ X 5 X 20 


54000 / 10800 


3.15^ 


52 X 8 X 20 


47600 / 9520 


3.30 


122 X 6 X 3Q 


35600 / 7120 


3.50'^ 


122 X 8 X 24 


39400 / 7880 


3.50^ 


142 X 6 X 30 


36800 / 7360 



Table 6.1: Run parameters for the torclon measurements using the ^ = 2 perfect 
action. The lattice extensions in torelon direction L, the extension in the two 
transversal spatial directions S as well as the temporal extension L are given. 



The Polyakov line around the short spatial direction is measured on APE- 
smeared configurations of smearing levels 3, 6, 9, 12, 15 for /? > 3.3 and 2 , 4, 6, 
8, 10 for /? < 3.15 with smearing parameter As — 0.1 (see Section 2.1. 2|) . The 



measured energies E{p'^) determined using variational methods (see Appendix 
^ depending on the lattice momentum are given, together with the number 
of operators used in the variational method, the fit ranges and the values of 
per degree of freedom, x^/A^df, in Tables G.l, p.2| in Appendix ^ The two 



simulations denoted by 3.15'' and 3.50'' are used to check the stability of the 
method. The first one is carried out on a lattice with small transversal spatial 
size (the torelon string measured is even longer), the latter is measuring the 
energies of a rather short string. 

Figures |6.2[ -]6?^ display the dispersion relations at the different values of the 
coupling [3 including all values of p at which the energies could be determined 
on the lattices considered. Note, that the energies of the p = torelons (the 
torelon masses) may be hard to determine because the effective masses do not 
reach a plateau within the temporal extent of the lattice; however, using the 
non-zero momentum energies, the masses may still be accurately determined. 

To determine the renormalised anisotropy as well as the torelon mass in 
units of the temporal lattice spacing mat we employ fully correlated fits over 
different ranges of p^ and choose the values depending on x^/^df and the 
precision of the dispersion relation data in the respective range. The results are 



given in Table 3.2 



Using the finite size relation for the torelon mass corresponding to eq. 2.1 



MriLa,) = (a + -J^)Las (6.1) 

it is possible to calculate the string tension a and thus the spatial scale Os . The 
string model predicts D = — | for long strings. To estimate the scale using 
the results of the torelon simulation, we use this string picture value and stay 
aware that the estimate gets worse for short strings. To obtain the hadronic 
scale To and the lattice spacings we employ ro^/CT=l-193(10) from Q and use 
the definition rg = 0.50 fm. The results are collected in Table |6.3|. 
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p 


fit range 




TUT at 




3.00 


0..5 


1.903(81) 


1.324(98) 


0.39 


3.15" 


0..4 


1.966(39) 


0.700(35) 


0.83 


3.15'' 


0..5 


2.022(104) 


1.262(84) 


0.62 


3.30 


1..9 


1.912(9) 


0.311(5) 


0.94 


3.50" 


1..10 


1.836(9) 


0.149(10) 


0.94 


3.50'' 


1..8 


1.826(16) 


0.208(16) 


1.61 


Table 6.2: 


Results of the torelon simulations 


using the ^ 


— 2 perfect action. 


The fit ran 


ge in is 


given in units of {2tt/SY . 





The scale is measured accurately using the static quark- antiquark potential 



(see Section 6.3) such that we can compare our torelon estimates to the true 



value. This is done in Table 6.4 and it turns out that the deviation of the torelon 
estimates for the scale ro/flg to the much more reliable potential values docs not 
exceed about 6% if the length of the string is not too short (say Lag > 1 fm). 
For short strings, the string picture does not apply, we may therefore not expect 
good results in this case. But even for very long strings, it is not possible to 
obtain very accurate information about the (hadronic) scale. Firstly, long strings 
are very heavy and thus difficult to measure (comparable to the very long range 
region of the static quark-antiquark potential), especially when the anisotropy 
^ is not very large; and secondly, using the string tension a to set the scale is 
plagued by ambiguities. Depending on the analysis of the respectiv e pote ntial 
measurements, rQ^/a may vary as much as about 4% (see Section |3.4.3 ) and 



additionally, the effective string tension at separations around r = 1 fm does 
not correspond to a occurring with the torelon strings. 

There seems to be no problem if the transverse volume is rather small as in 
the /3 = 3.15'' simulation, except of course the large momenta occurring which 
make the determination of the energies more difficult. As well, the estimates 
of the renormalised anisotropy do not show significant deviations neither for 
small transverse volume nor for short strings. 



a 






as [fm] 


at [fm] 


3.00 


0.834(20) 


1.43(5) 


0.350(12) 


0.184(14) 


3.15" 


0.563(8) 


2.12(5) 


0.236(6) 


0.120(5) 


3.15'' 


0.579(18) 


2.06(8) 


0.243(9) 


0.120(11) 


3.30 


0.3578(21) 


3.33(5) 


0.150(2) 


0.079(1) 


3.50" 


0.225(5) 


5.31(16) 


0.094(3) 


0.051(2) 


3.50'' 


0.304(7) 


3.92(12) 


0.128(4) 


0.070(3) 



Table 6.3: Estimates of the scale determined from torelon results (see Table 



3.2), see text. 



Finally, we may conclude that determining the renormalised anisotropy using 
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p 


Las [fni] 


ro/ag (Torclon) 


j'o/fls (Potential) 


rel. error 


3.00 


1.48 


1.43(5) 


1.353(16) 


+5.7% 


3.15" 


1.23 


2.12(5) 


2.038(2) 


+4.0% 


3.15^ 


1.96 


2.06(8) 


2.038(2) 


+1.1% 


3.30 


0.95 


3.33(5) 


3.154(8) 


+5.6% 


3.50" 


0.82 


5.31(16) 


4.906(21) 


+8.2% 


3.50^ 


0.61 


3.92(12) 


4.906(21) 


-20.0% 



Table 6.4: The estimates of the scale obtained from torelon measurements (see. 



Table 6.3) compare d to the scale measured using the static quark- antiquark 



potential (see Table x6 ) . The length of the torelon string Lag in physical units 
is given as well. 




1 2 3 

(a p)" in units of (2JI/8)' 



Figure 6.2: Torelon dispersion relation for /? = 3.0. The straight line is the 
correlated fit to E'^{p) = rn^ + in the range — 0..5. 
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(a p) in units of (27C/8)* 



Figure 6.3: Torelon dispersion relation for (3 = 3.15". The straight hne is the 
correlated fit to E'^{p) — rn^ + p"^ in the range — 0..4. 
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Figure 6.4: Torelon dispersion relation for /3 = 3.15*' on a small transverse 
volume. The straight line is the correlated fit to E'^{p) = +p'^ in the range 
p2 = 0..5. 
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(a p)^ in units of (llUuf 



Figure 6.5: Torelon dispersion relation for (3 = 3.3. The straight line is the 
correlated fit to E'^{p) = m\ in the range = 1..9. 




Figure 6.6: Torelon dispersion relation for (3 = 3.5°. The straight line is the 
correlated fit to E^{p) = m\ in the range = 1..10. 
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6 8 10 12 
(a p)" in units of (2JC/14)" 



Figure 6.7: Torelon dispersion relation for (3 = 3.5^ for a short torelon string. 
The straight Hne is the correlated fit to E'^{p) = rn^ in the range = 1..8. 



the torelon dispersion relation is a stable and apparently sensible procedure. 
There are no manifest problems and it works for coarse and fine lattices as well. 
Checks employing different lattice extensions indicate that the determination 
of the renormalised anisotropy is very stable. The estimate of the lattice scale, 
however, is not very accurate, mainly due to the use of the string picture relation, 



eq. 3T; the (systematic) error of the torelon scale is about 5% for sensible torelon 
lengths L > 1 f m with a tendency of underestimating the lattice spacing a^. 

Choosing the lattice size, one has to pay attention on one hand to the length 
of the torelon string, i.e. the extension of the (short) longitudinal direction: if 
the string is too short, there is no useful estimate for the scale; if it is much 
too long, it will be difhcult to measure due to its high energy. On the other 
hand, the transversal direction in principle may be short in order to make the 
simulation cheaper what concerns computational effort, however, if it is too 
short, the momenta are large, and again the measurement becomes difficult. 
Finally, the temporal extension should be chosen large enough such that the 
torelon mass (which may be rather light) may be determined directly. 

6.3 The Scale 

To set the scale of the action at a given /?-value we measure the quark-antiquark 



potential (see Section 2.3) and determine the Sommer scale rg ~ 0.50 fm. For 



the anisotropic action, we still employ the definition of rg: 

r''V'{T)\r^ro = 1-65, (6.2) 
however V{r) is now measured in units of the temporal lattice spacing and thus: 

-a + ^CT( — )2 = 1.65, (6.3) 
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where a denotes the fitted dimensionless value a — aagat of the string tension. 
The renormalised anisotropy ^ thus has to be known before the spatial scale 
may be determined. 

The static quark-antiquark potential is measured on lattices of different 
scales. The measurements &% (3 — 3.00, 3.30 include the determination of the 
potential between quarks that are off-axis separated along the lattice vectors 
(1,0,0), (1,1,0), (1,1,1), (2,1,0), (2,1,1), (2,2,1) (and lattice rotations) in order 
to estimate violations of rotational symmetry. The rest of the measurements 
include only on-axis separations due to the large computational cost (concern- 
ing speed and memory) of th e of f-axis measurement. The parameters of the 
simulations are given in Table [S.Sj 



/3 


off- axis scp. 


5*3 X T 


# sweeps / measurements 


3.00 


yes 


8^ X 16 


42000 / 2800 


3.15 


no 


10^ X 20 


39800 / 3980 


3.30 


yes 


10^ X 20 


27000 / 1800 


3.50 


no 


12^ X 24 


36400 / 3640 



Table 6.5: Run parameters for the measurements of the static quark-antiquark 
potential using the ^ = 2 perfect action. 



T he va lues of the off-axis potential atV{r) a t (3 — 3.30 are collected in 
Table G.4 in Appendix ^ and displayed in Figure 6^, the values of t he o ff-axis 



The 



G.5 



potential on the coarse lattice at /? = 3.00 are collected in Table p. 5 
potential values atV(r) of the on-axis simulations are collected in Table 

The parameters a and a (the string tension) in the Cornell potential, eq. 2.24 , 
are determined using global fits, the hadronic scale tq is determined performing 
local fits. The fit ranges and results are given in Table |6.6| . In the global fit to 
the off-axis potential at /3 = 3.30 we exclude the badly measured separations 
2(2,1,0) and 2(2,1,1). 



p gl. 


fit 


a 




loc. 


fit 






3.00 1 - 


2\/6 


-0.201(18) 


0.3595(87) 


^/6^ 


~3 


1.353(16) 


1.119(51) 


3.15 1- 


~ 5 


-0.1503(12) 


0.1654(6) 


1 - 


3 


2.038(2) 


1.162(15) 


3.30 ^/2 


- 5 


-0.1539(7) 


0.0683(3) 


2%/2- 


2V3 


3.154(9) 


1.140(8) 


3.50 1- 


- 6 


-0.1478(4) 


0.0320(3) 


2- 


4 


4.906(21) 


1.189(14) 



Table 6.6: Results of the measurements of the static quark-antiquark potential 
using the ^ = 2 perfect action. The parameters a and a of the Cornell potential 
are given together with the global fit range chosen to determine them. The 
hadronic scale rg is given in spatial lattice units together with the local fit 
range used to determine that quantity. Note, that c = 6.00 has been chosen 
for (3 = 3.0 and c = 0.89 has been chosen for j3 — 3.5 (see Section 2.3.4). 
Additionally, the dimensionless quantity tq^/u is given. 



In principle, it is possible to determine the renormalised anisotropy using 
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Figure 6.8: Ground state of the static quark-antiquark potential at (3 
Including a global fit to the Cornell potential in the range r — \/2..5. 



3.3. 



the first excited state of the potential together with the ground state. As can 
be seen from Figure |6.9| however, for the off-axis potential at /? = 3.3, the 
energy values of the first excited state have large errors. The result for the 
separation at(F*(ro) — F(ro)=0. 555(47) thus shows a large error. Comparing 
this value to ra{V*[ro) - V{ro)) « 3.25(5) of |l4|l one obtains ^ = 5.86(59) 
and thus = as/at=l. 86(19) for [3 = 3.3 which is in agreement with the 
value determined using the torelon dispersion relation, ^^=1. 912(9). A similar 
determination for the on-axis potential at /3 = 3.50 is even more difficult due 
to the small number of data. The result is ^r=1. 66(32) which does again agree 
with the torelon result, — 1.836(9). Due to the large errors (caused by the 
choice of (3 and the operators used) this method of determining the renormalised 
anisotropy ^r is not suitable in the context of this work. However, for the 
determination of the scale as well as of the anisotropy on lattices with finer 
temporal lattice spacing (above all with higher anisotropics) this way might be 
feasible at least if off- axis separations are included in the measurement. 

The measurements at /? = 3.0, 3.3 including a large number of off-axis separa- 
tions provide information about the deviations from rotational invariance. This 
issue has been addressed in |60| for the renormalisation group used in this work. 
It turned out that violations of rotational symmetry caused by the blocking are 
small (see also Section 5.2.1 ). That this is still true f or t he parametrisation em- 
ployed is shown by the off-axis potential, see Figure 6.S. The devia tions on the 
coarse lattice at /3 = 3.00 are significantly larger (see Figure |6.10 ) which is no 
surprise as this coupling corresponds to a spatial lattice spacing as ~ 0.37 fm. 

The static qq potential is an effective test of scaling. Expressing the po- 
tential measurements performed at different couplings (3 in the RG invariant, 
dimensionless ratios r/rp and ri^V , subtracting the unphysical constant raV{rQ) 
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Figure 6.9: Static quark antiquark potential at (3 — 3.3 including the first 
excited state V*{r) at /3 — 3.3. The ground state values are again fitted to the 
ansatz V{r) = Vq + a/r + ar, whereas the values of the excited state are fitted 
to the (ad-hoc) ansatz V*{r) ~ A/r + B + Cr + Dr^ . The value of ro where the 
determination of the anisotropy is done, using the gap, is marked by the dotted 
line. 
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should lead to potentials lying exactly on top of each other. Deviations indicate 



either scaling violations or ambiguities in the determination of tq. Figure 6.10 



which includes all potential measurements shows that the different curves can 
hardly be distinguished from each other, except some energies determined on 
the coarse lattice at /? = 3.00 that deviate notably from the curve which is a fit 
to the Cornell potential including all the measurem ents that have been included 
into the global fits of the single f3 values, see Table |6.6| . The results for the two 
physical parameters are a = —0.27799(5) and crrQ = l. 3690(3). 
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Figure 6.10: The static quark- antiquark potential measurements expressed in 
RG invariant, dimensionless units r/vQ, roV. The unphysical constant roV^(ro) 
has been subtracted such that all the curves exactly coincide at r — ro- The 
dashed line is a global fit to the Cornell potential, eq. 2.24. 



6.4 The Critical Temperature 

Another interesting quantity to study with the ^ = 2 perfect action is the critical 



temperature of the deconfining phase transition, see Section 2.2. This quantity 
contains mainly information about the temporal scale of the lattice at a given 
coupling /3 as Tc = \/{Ntat). Comparing this information to quantities obtained 
from the measurements of torelons or the static quark-antiquark potential offers 
a lot of possibilities of scaling (and other) tests. Because measurements of these 
quantities at the determined values of /3crit are yet absent and because interpo- 
lations in (3 are very difficult, because of the renormalisation of the anisotropy 
^ which is another parameter entering, these tests are a future project. 

In this first study, we decide not to examine the finite-size scaling of /?crit 
but to choose Lg/Lt = ^ ■ Ng/Nt « 3.5 ~ 4. Obviously, the anisotropic nature 
of the action makes the computational effort of the simulations for exploring 



80 



Chapter 6. The (, — 2 Perfect Action 



the deconfining phase transition smaller. The effects of the finite volume are 
(rather conservatively) estimated by using the finite-size scaling relation for 
Pcrit, eq. The parameter \h\ is set to the largest value \h\ = 0.25 that has 
been observed with the isotropic action, the value of ^r, used to estimate the 
spatial lattice volume, is set to the lowest value appearing in the whole range 
of (3 values considered. These finite- volume errors are given together with the 
run- t ime p arameters and results of the ^ = 2 simulations in Table 6/7. Figures 
3.11 - 3.15 show that the determination of Pc is rather precise, partly due to the 
reweighting technique. As well, it shall be understood that in order to make the 
results even more precise, a finite-size extrapolation is indispensable. 




2.8 2.82 



2.9 2.92 



Figure 6.11: The Polyakov loop susceptibility on the 3x6'^ lattice for the 
^ = 2 perfect action. The solid curves are the interpolations using Ferrenberg- 
Swendsen reweighting, the dashed lines indicate the bootstrap error. 
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lattice size 


P 


sweeps 


Tp 


Tint 


XL 


/3c 


3 X 6^ 


2.80 


45000 




70 


1.45(5) 


2.863(1)(5) 




2.85 


45000 


410 


259 


11.3(20) 






2.87 


45000 




312 


21.6(28) 






2.90 


50000 




100 


2.34(10) 




4 X 8^ 


3.00 


45000 




81 


1.63(5) 


3.032(1)(5) 




3.02 


45000 


870 


102 


3.97(26) 






3.03 


45000 


830 


372 


43.2(22) 






3.04 


45000 




188 


10.9(18) 






3.05 


55600 




144 


3.94(21) 




5 X 9^ 


3.07 


44000 




117 


0.84(4) 


3.118(1)(6) 




3.10 


45000 


560 


154 


4.38(63) 






3.11 


20500 


640 


198 


15.0(27) 






3.12 


44500 


3700 


278 


33.2(19) 






3.15 


20000 




110 


3.26(18) 




6 X 11^ 


3.17 


45000 


1440 


121 


7.4(11) 


3.181(1)(6) 




3.18 


45000 


1750 


201 


30.0(15) 






3.185 


45000 




232 


23.2(24) 






3.19 


45000 




123 


10.4(14) 






3.20 


45000 




100 


5.62(36) 




7 X 133 


3.22 


32400 


700 


87 


6.7(5) 


3.236(1)(6) 




3.23 


44500 


750 


129 


18.4(11) 






3.24 


29000 




149 


20.7(16) 






3.25 


22400 




83 


9.2(10) 





Table 6.7: Run-time parameters of ^ = 2 Tg-simulations. The number of sweeps, 
the persistence time Tp (if definable), the integrated autocorrelation time Tint 
and the Polyakov loop susceptibility Xh arc given for all values of /? used in the 
final reweighting procedure, as well as the resulting /3c . The systematic error 
in the second brackets is a (rather conservative) estimate of the finite volume 
effects and has to be considered if the given value should be an estimate of the 
infinite-volume critical coupling. 
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3 3.02 3.04 3.06 



P 



Figure 6.12: The Polyakov loop susceptibility on the 4x8^ lattice for the 
^ = 2 perfect action. The solid curves are the interpolations using Ferrenberg- 
Swendsen reweighting, the dashed lines indicate the bootstrap error. 




Figure 6.13: The Polyakov loop susceptibility on the 5x9'^ lattice for the 
£, — 2 perfect action. The solid curves are the interpolations using Ferrenberg- 
Swendsen reweighting, the dashed lines indicate the bootstrap error. 
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Figure 6.14: The Polyakov loop susceptibility on the 6 x 11^ lattice for the 
^ = 2 perfect action. The solid curves are the interpolations using Ferrenberg- 
Swendsen reweighting, the dashed lines indicate the bootstrap error. 




Figure 6.15: The Polyakov loop susceptibility on the 7 x 13^ lattice for the 
^ = 2 perfect action. The solid curves are the interpolations using Ferrenberg- 
Swendsen reweighting, the dashed lines indicate the bootstrap error. 
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6.5 Glueballs 

Finally, our ^ = 2 classically perfect action is used to determine the glueball 



spectrum in pure lattice gauge theory (cf. Section 2.4). As glueballs are rather 
heavy, mc ^1-5 GeV, to study their masses the use of anisotropic lattices is 
particularly useful. 

Due to mixing with other states existing in pure gauge theory on a peri- 
odic lattice, such as torelons or glueball pairs, and due to the breaking of the 
continuum rotational symmetry by the lattice which leads to glueball states 
transforming under the cubic group (see Section [2.4| ), the identification and the 
continuum spin assignment of single glueball states are additional vital questions 
which will be addressed in Sections |6.5.3 , 6.5.5| . Finally, the resulting energies 



of all the states that can be measured and assigned a value of continuum spin 
to, are given in units of the hadronic scale as well as in MeV for convenience. 

6.5.1 Simulation Parameters 

To allow a continuum extrapolation, at least for the lighter states, we de- 
cide to perform simulations at three different lattice spacings in the range 
0.10 fm < fls < 0.25 fm in volumes between (1.2 fm)^ and (2.0 fm)^. The 



simulation parameters are given in Table 3.S. The gauge fields are updated 



by performing compound sweeps consisting of 4 pseudo over-relaxation and 1 



Metropolis sweep (see Sections 3.1.1, 3.1.2), and after every 2 compound sweeps 



we measure all 22 loop shapes up to length 8 (see Figure 2.7). The measurement 
is performed on APE-smeared configurations (see Section 2.1.2) to spatially en- 
large the operators, thus improving the overlap with the glueball states, and 
to reduce unphysical high- momentum fluctuations. On the coarse lattice at 
P = 3.15 we use smearing levels 2, 4, 6, 8, 10, on the other lattices we use levels 
3,6,9,12,15, always with smearing parameter As = 0.1. The measured oper- 
ators are projected into all the irreducible representations of the cubic group. 
We measure all possible polarisations to suppress statistical noise (see Section 

EH)- 
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S^xT 
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3.15 


0.2408(5) 


8^ X 16 


(1.96 fm)3 


241000 / 24100 


3.30 


0.1005(6) 


10^ X 20 


(1.59 fm)3 


99000 / 9900 


3.50 


0.0415(4) 


123 X 24 


(1.22 fm)3 


115000 / 11500 



Table 6.8: Run parameters for the glueball simulations using the ^ = 2 perfect 
action. 



6.5.2 Determination of the Energies 

The masses of the lowest states and the first excited states (if possible) in all 
the representations are determined using the variational methods and techniques 
described in Appendi x |^. The results in units of the temporal lattice spacing 
at are listed in Tables G.6 - G.8 in Appendix]^ together with the number N of 
operators entering into the fitting procedure, the time slices on which the initial 
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generalised eigenvalue problem is solved (usually 1/2), the number of operators 
M in the final fit and the corresponding value of x^/-^df- These results are then 
multiplied by ■ tq/ to obtain the glueball masses in units of the hadronic 
scale ro, listed in Table 6.£. 

In order to obtain reliable estimates of the glueball masses from the varia- 
tional method, one has to pay attention that among the large number of op- 
erators (up to 145 on five smearing levels) there are no correlators entering 
the process that are measured exceptionally bad, i.e. with large errors even at 
small time separation, as this may destabilise the method. To filter out such 
data, we calculate the "correlator lifetimes" defined in Section 2.4.2 and drop 
all correlators whose lifetime is smaller than a given minimum. Usually, this is 
repeated with different minimal values, yielding different sized sets of operators 
(2 ~ 4 sets per representation). The O(IOO) mass determinations using differ- 
ent parameter sets (see Appendix ^ are then carried out on each set of initial 
operators, and it turns out that generally the variational method used is stable, 
provided that the badly measured operators are absent from the beginning and 
provided that no important operators are missing. Patently, having such a win- 
dow of the number of initial operators may be difficult if there is a rather small 
number of operators. 

For some heavy states, such as the PC = glueballs, whose masses could 

only be determined from the /3 = 3.50 measurements, we have to resort to using 
the correlators at to = and ti = 1 for the solution of the first generalised 
eigenvalue problem, eq. E.3, yielding the lowest lying dual vectors n = 
1, . . . , M, that are used in the subsequent optimisation. This has the advantage 
that the correlators from separation to = are positive definite by definition 
(which is not true ioi to > I), thus rendering the second generalised eigenvalue 
problem, eq. E.8, well defined even if all initial operators are kept, M = N. 
However, the correlators C{t = 0) are under suspicion of containing rather little 
physical information about the correlation lengths (which correspond to the 
decay of the correlators), actually they almost only provide information about 
the relative normalisation of operators. Obviously, operators with a large signal 
at t = might be preferable because their signal/noise ratio is probably rather 
better on average, however in general they do not correspond to the operators 
having the largest overlap with the lowest-lying states. It turns out that using 
to = 0, ti = 1 one has to pay much more attention in choosing the fit range 
t = ^min • • ■ ^max because at small t the contamination of the lowest masses due 
to higher states is much more significant. 

During the analysis, it turns out that the mass determinations for the scalar 
representation Af~^ are rather challenging. The ground state receives a larger 
relative error than comparable states where the operators are measured about 
as well; for the first excited state it is very hard to obtain a stable determination, 
and the error of its mass may turn out to be huge (e.g. at /3 = 3.50). Probably, 
it is the underlying vacuum, having the same quantum numbers as the glueball 



state, that is responsible for these troubles. As already specified in Section 2.4.2 



we treat the vacuum just like another state in this representation, such that 
the glueball ground state is effectively a first excited state and the glueball first 
excited state is effectively a second excited state of the representation. Attempts 
of using other ways of getting rid of the vacuum such as the usual subtraction or 
the subtraction of large t correlators (that are assumed to contain solely noise, 
however correlated to the noise at lower t) or even more sophisticated methods 
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(like solving a generalised eigenvalue problem using large t correlators to dig 
out the vacuum state) do not at all succeed in improving the situation. 



6.5.3 Identification of Glueballs 

Besides single glueballs, the spectrum of pure gauge theory on a lattice with 
periodic boundary conditions also contains states consisting of several glueballs, 
torelons or mixed states of glueballs and torelons. Although we expect that 
the operators used to measure glueball energies couple most strongly to the 
single glueball states, other states with similar masses and compatible quantum 
numbers might mix with them and distort the result, i.e. the energy determined 
by the analysis of the correlation matrix may be dominated by a multi-glueball or 
torelon state with smaller energy than the single-glueball state to be measured. 

In principle, there are several means of determining the nature of a state 
that has been measured. Firstly, the simulation may be repeated on lattices of 
different size, keeping the lattice spacing fixed. As multi-glueball and torelon 
states show a finite-volume behaviour very different from single glueballs, such 
states stand out and may be dropped from further analysis. Secondly, one may 
measure additional operators that couple strongly to torelon or multi-glueball 
states. Including and excluding these operators in the variational method and 
studying the coefficients obtained from the optimisation, the mixing strength is 
determined and states which do not mix considerably with any of the additional 
operators may be safely considered to be single-glueball states. Finally, using the 
mass estimates for the low-lying glueballs, one may determine the approximate 
locations of the lowest-lying multi-glueball states. Also, the minimal energy of 
mixing torelons may be estimated using the string formula, eq. 



Torelons 

The first two of the methods mentioned above require additional work and com- 
puter time and go beyond the scope of this work. The last method, however, 
can be done rather easily. Let us first calculate the energy of the lowest-lying 
torelon states that may interfere with our measurements. Single torelons (see 



Section 2.1) transform non-trivially under Z3 symmetry operations; our opera- 
tors, closed Wilson loops, however, are invariant under these transformations. 
This means that they cannot create single torelon states, however the creation 
of two torelons of opposite center charge is possible. We assume that they do 
not interact considerably and that our lattice extension is rather large and thus 
we use the simple formula E2T ~ 2aL to estimate the energy of a torelon pair 
with momentum zero, where a is the string tension (measured by the static qq 



potential, see Section 3.3) and L is the spatial extent of the lattice. Table [6.10 



lists the minimum energies of torelon pairs to be expected on our lattices. 

Note, that a state composed of two opposite center charge torelons and with 
total zero momentum is symmetric under charge conjugation. The operators 
for C = — states therefore do not create such torelon pairs. 



Multi-Glueball States 

To estimate the lowest energies of multi-glueball states present on our lattices, 
we follow the method used by Morningstar and Peardon fl^ , described in detail 
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Tabic 6.9: Final glucball mass estimates in terms of the hadronic scale tq. 
mcro from the measurements using the perfect ^ = 2 action. The continuum 
spin assignment J is given as well. 
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a 


r^E2T 


3.15 


10.60 


3.30 


8.24 


3.50 


6.92 



Table 6.10: The minimum energies of momentum zero torelon pairs on the 
lattices used in the glueball simulations, given in units of r^^ . 



in Appendix ^. We assume that the energy of multi-glueball states is approx- 
imately given by the sum of the energies of the individual glueballs, i.e. that 
there is no substantial interaction. It is therefore clear that the lowest-lying 
multi-glueball states with zero total momentum are the two-glueball states, 
with energy 



E2G^ \Jf + mi + +mi, (6.4) 

where mi denotes the rest mass of the first glueball with momentum p and 
m2 denotes the rest mass of the second glueball with momentum —p. The 
masses are taken from our determinations at the respective values of /3, the 
(lattice) momenta are chosen such that the energy E2G of the glueball pair 
contributing to some representation r^*-^ is minimised. Note, that glueball pairs 
contribute to different representations of the cubic group, not only depending on 
the representations according to which the single glueballs transform, but also 
on the single glueball momentum p. It happens therefore, that a pair of rather 
light glueballs with non- vanishing momenta p has a smaller energy than the pairs 
of glueballs at rest contributing to the same representation F^'-^. The details of 
the analysis are collected in Appendix ^ the resulting minimal energies in units 
of ro of glueball pairs in all the repre sentations where states can be measured 
on our lattices are collected in Tables |a9| - |G.ll| . 



The lower bounds of the multi-glueball energy region as well as the lower 
bound for torelon pairs are indicated in Figures |6.16 - 3.1^ . It turns out that 



all of the excited states measured and even some of the lowest states in a given 
representation could be affected by torelon pairs or multi-glueball states. Of 
course, it is not at all ruled out that these states are indeed single glueballs, 
however one has to be very careful with the interpretation and has to keep in 
mind that these issues require further study. 



6.5.4 The Continuum Limit 

In order to remove discretisation errors from glueball masses romc that have 
been obtained on lattices with finite lattice spacing, we have to perform the con- 
tinuum hmit, i.e. to extrapolate the results to — 0. Having obtained results 
at three values of the coupling /3, one of them corresponding to rather coarse 
lattices, and having no accurate information (e.g. from perturbation theory) 
about the behaviour of the energies depending on the lattice spacing (about the 
form of the curve used in the extrapolation) this is not an easy task and rather 
ambiguous. 
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Figure 6.16: The results of the glueball measurements at /3 = 3.15 converted 
into units of tq. The lower bounds of the energies of multi-glueball states (solid 
lines) or torelon pairs (dotted line) are given as well. 




Figure 6.17: The results of the glueball measurements at /3 = 3.30 converted 
into units of rg. The lower bounds of the energies of multi-glueball states (solid 
lines) or torelon pairs (dotted line) are given as well. 
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Figure 6.18: The results of the glucbaU measurements ai P — 3.50 converted 
into units of tq. The lower bounds of the energies of multi-glueball states (solid 
lines) or torelon pairs (dotted line) are given as well. 



We perform two different kinds of continuum limits. Firstly, we extrapolate 
the masses in terms of the hadronic scale tq, tqitiq = mQat ■ ^ • rQ/a^, secondly, 
we extrapolate mass ratios rnGi/mcs of different glueball representations. Gen- 
erally, the procedure we resort to, is the following: We include all the three 
energies or mass ratios, obtained at /3 — 3.15, 3.30, 3.50, into the fit and use the 
form 

romc \a,^ roma \a,=o +C2 [ —] ■ (6.5) 

Additionally, we perform fits including only the results from the two finer lattices 
at l3 — 3.30, 3.50 using a constant (if this fits the results reasonably well) or the 
form stated above. This value of the continuum mass (or mass ratio) is then 
compared to the result of the fit of all the three data points, which is only 
accepted if the two results coincide within their errors. 



Masses in Terms of Tq 

For the extrapolation of masses in terms of the hadronic scale ro, this procedure 
works very weh for the representations , E-+, E++ , E++* , and T+"^ 
where we may always fit the results from the two finer lattices to a constant 
which agrees with the result of the extrapolation using eq. 3.5 on all three data 
points. For the A^^ representation, we obtain the best results fitting the masses 
to a constant, for the fit to two data points as well as to three data points. 

Much more difficult is the extrapolation for the scalar glueball and its excited 
state, A^"'" and . It turns out that the fit has a very bad and describes 
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the data very badly if eq. 6.5 is used. However, we find good behaviour if we 
use 

romo \a,^ romc \a,=o +C6 { — ] (6.6) 

for the ground state as well as for the excited state energies and thus decide 
to use this (ad-hoc) fitting form for the scalar glueball. The results of the 



continuum extrapolations for various glueball states are collected in Table G.14 
In Figures |6.19 , 6.20| the curves are shown together with the measured values 



for the PC = ++ and the PC = — h sectors where extrapolations are possible 
with the present data. 

Other representations where energies can be determined at several values 
of /3, such as , T^^ or T^^ do not yield consistent results using the 
procedure presented above and we decide to abstain from performing continuum 
extrapolations in terms of romc for these states. 

Mass Ratios 

The other method of extrapolating to the continuum, namely using mass ra- 
tios mGi/"2G2 between different glueball representations, has got the advantage 
that uncertainties present in the determination of ro/og or the renormalised 
anisotropy do not enter into the continuum extrapolation. Furthermore, we 
may assume that several effects such as finite size scaling are more similar for 
different glueball states measured on the same lattice as opposed to completely 
different quantities such as the static quark-antiquark potential or the torelon 
dispersion relation measured on different sized lattices that are used to obtain 
ro/ttg and respectively. 

The common way of extrapolating glueball mass ratios to the continuum is 
using ratios mG/mQ++ (see Table G!.12| ), i.e. ratios to the scalar glueball mass. 



which is the lightest mass present. Doing this, we observe that the extrapola- 
tion is difficult, the goodness of the fit is rather bad and the errors are large. 
Furthermore, the resulting continuum mass ratios are not in good agreement to 
ratios of the continuum masses coming from the extrapolation of the masses in 
terms of tq. This is not a surprise, because the behaviour of the scalar glueball 
as a function of the lattice spacing is rather difficult, i.e. there are large cut-off 
effects that cannot be described too easily. Additionally, the errors of the mass 
estimates of the scalar glueball are quite large. Due to these reasons, we decide 
to use the well measured masses of the tensor T^~^ representation, which show 



small errors and seem to scale rather well. These ratios are listed in Table G.13 



The procedure of the extrapolation is the same as for the case of romc de- 



scribed in the previous section; the unambiguous results are listed in Table G.15 . 
Again, for the A^^ states, as well as for the E^~^, the fit to a constant works 
out best; the T^^ glueball can only be reasonably extrapolated using a fitting 



form analogous to eq. 6.6. Note, that this way we obtain continuum results for 
the representations , T.^~^ and T^'^* which is not possible performing the 
extrapolation using masses in terms of rg. Furthermore, all the other results 
listed are in agreement to the ratios of the continuum masses romc within one 
standard deviation. 

In the following, we will use the finite lattice spacing results from the simu- 
lation on the finest lattice at /? = 3.50 for the energies of states that cannot be 
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extrapolated to the continuum in one of the two ways, e.g. to study degeneracies 
or to draw a (rather quahtative) picture of the low-lying glueball spectrum. In 
tables listing continuum results, these values not extrapolated to = will be 
stated in square brackets. 




Figure 6.19: Mass estimates of the PC = ++ glueballs in terms of the hadronic 
scale ro against (os/ro)^. The curves are the continuum limit extrapolations of 
the form roTOG + C2(as/ro)^. Circles: A'^'^ , boxes: A^"*", diamonds: E^'^ , up- 
ward triangles: T^^ , downward triangles: T^^; solid symbols: ground states, 
open symbols: first excited states. Note that the continuum results of the rep- 
resentations T^~^ and T^~^ have been obtained using a fit of glueball mass 
ratios. 



6.5.5 Continuum Spin Identification 

Once the extrapolation of the measured glueball energies to the continuum, 
as = 0, has been performed, what remains is the assignment of continuum spin 



J to the different representations of the cubic group T. From Table 2.2 we know 
from which lattice representations F^*-^ a continuum glueball with quantum 
numbers J^*^ may obtain contributions, where the assignment of parity P and 
charge conjugation C is simply 1:1. Degeneracies (in the continuum limit) 
between several representations F'^*-^ contributing to the same continuum state 
J-^'-^ are a strong indication for the correctness of the assignment of all the lattice 
states involved to the same continuum J^'~^ glueball. Furthermore, we make the 
assumption that the mass of the glueballs increases with the spin getting larger. 
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Figure 6.20: Mass estimates of the PC = — + glueballs in terms of the hadronic 
scale rg against (os/ro)^. The curves are the continuum limit extrapolations 
of the form romo + C2(as/ro)^. Circles: , diamonds: E ^, downward 

triangles: T2~^', solid symbols: ground states, open symbols: first excited states. 
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In the PC = ++ sector, we observe a single low-lying state, Af~^, which 
is thus assigned J = 0. The E~^~^ and T^~^ states are degenerate to a very 
high precision as it should happen if these representations correspond to the 
five polarisations of a J = 2 glueball. The excited state A^^ again has no 
degenerate partner and is assigned to an excited state of the continuum J — 
glueball. The excited states and T^'^ again turn out to be degenerate 

and so they are assigned to an excitation of the J ~ 2 glueball. The same 
situation is met with the representations and T^^ which are thus assigned 
J = 3. The first excitation T^'^* being part of the J = 2 state, the missing 3 
polarisations should come from the second excitation of T^^ which cannot be 
measured by our simulations. 

In the PC = — h sector, the situation of the A^^ and the E ^ and Tj"^ is 
very similar to the one of their partners in the PC = ++ sector; these states are 
thus assigned J = and J = 2 respectively. The remaining three, excited states 
have no apparent degeneracies, there is thus no safe assignment of continuum 
spin to these states. The facts that the excited state A^^ is rather light and 
comes with no degenerate partner, raises the presumption that it is the excited 
J = 0-+ gluebaU. 

In the other two sectors, continuum extrapolations of our measured energies 
are not possible, except for the representation using mass ratios. However, 
looking at the degenerate representations in the PC = ++, — h sectors, it is 
noticeable that the degeneracies are apparent even for finite lattice spacing at 

P = 3.50,3.30. Assuming that this behaviour persists for PC = H — , , at 

least for energies not too high, we may assign continuum spin even to some of 
the remaining states. 

In the PC = H — sector, we notice the lowest-lying state which has no 
degenerate partner which suggests a J = 1 interpretation. Next, there is the 
A2~ which is likely to correspond to J — 3. 

Finally, in the PC = sector, there are different possible scenarios. The 

almost exact degeneracy of and suggests them contributing to the 

J = 3 continuum state. The rather heavier state could still be degenerate 
with the two latter states (note the high mass and the finite lattice spacing!) 
and carry the remaining polarisation. Admittedly, it could as well correspond 
to an excited J = 3 state or even to J = 6. The final A^^ state is very heavy 

and indicates that glueballs with even spin and PC = have large energies 

as the representation Ai contributes solely to even numbered spin states up to 
J = 8! 

The rather well supported continuum spin assign ments pre s ented above are 
given together with the glueball masses in Tables 6^, G.14 , |G.15 , the final 
results on the masses of these states are collected, together with the masses in 



MeV, in Table 3.11. Figure 6.21 provides a rather qualitative picture of the 



low- lying glueball spectrum as obtained from our ^ = 2 simulations. 



6.5.6 Discussion 

In order to discuss the results obtained from the glueball measurements, pre- 
sented in the previous Sections, let us first list the sources of errors that may 
be present. Firstly, there are finite size effects; our finest lattice (/3 — 3.50) 
has got a spatial volume of (1.22 fm)'^ which is rather small compared to other 
volumes used in glueball measurements. The magnitude of finite-size effects de- 
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Figure 6.21: The mass spectrum of glueballs in SU(3) pure gauge theory as 
obtained from our simulations using the ^ = 2 classically perfect action. The 
masses are given in units of the hadronic scale tq. The dashed boxes are used 
for states that could not be extrapolated to the continuum, instead the mass 
obtained from the simulation on the finest lattice at /3 = 3.50 is given. The 
most probable assignment of continuum spin is indicated as well. 
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Table 6.11: Final results for the masses of continuum glueballs with quantum 
numbers J^'^ , obtained from the lattice representations F^*-^. For the conversion 
to MeV, w 0.5 fm w 395 McV has been used. Quantities in brackets have 
not been extrapolated to the continuum but denote results from the simulation 
on the finest lattice at /? = 3.50 (corresponding to Us = 0.102 fm). 



pends largely on the quantity studied, more precisely on the form of the glueball 
wavefunctions (above all their extension). Due to the three volumes employed 
having such different size, finite-size effects do not only systematically shift the 
lattice results but they can also make reliable continuum extrapolations much 
more difficult or even impossible. 

Secondly, the mass of the scalar glueball Af~^ (above all the ground state) 
shows very large cut-off effects. These could be due to the presence of a criti- 
cal end point of a line of phase transitions in the fundamental- adjoint coupling 
plane assumed to define the continuum limit of a scalar field theory. Our action 
includes in its rich structure operators transforming in the adjoint representa- 
tion. If their net coupling (which we do not control during the construction and 
the parametrisation of the action) lies in a certain region, the effect of the crit- 
ical end-point on scalar quantities at certain lattice spacings (sometimes called 
the "scalar dip") may even be enhanced compared to other (more standard) 
discretisations with purely fundamental operators. 

Furthermore, our analysis shows that states consisting of several glueballs 
or torelon pairs could mix with most of the higher-lying states occurring in 
our measurements. This makes the interpretation of the measured states more 
difficult, additionally, systematic effects (e.g. different mass shifts at different 
values of /3 because of the different lattice sizes or due to the presence of op- 
erators having larger overlap with unwanted states) may again complicate the 
continuum extrapolation. 

Finally, as always, there is the statistical error. Generally, it does not intro- 
duce bias and it can be estimated quite reliably. However, glueball correlators 
are notorious for large statistical noise and this noise may hamper the reliable 
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variational analysis of states that in principle could still be measured as for their 
energies. 

There are several possibilities of improving the measurements. Firstly, and 
most simple, improving the statistics may help. This can be seen e.g. looking at 
the results of the simulation on the coarsest lattice at f3 — 3.15. Although the 
lattice is coarser, more energies can be determined than on the finer P = 3.3- 
lattice, mainly due to the larger statistics. Additionally, a lot of problems 
stated previously are solved (or at least kept under control) if for each value 
of /3, we determine the glueball energies on several lattices of different physical 
size. Studying the finite-size scaling creates the possibility of separating single 
glueball from other (multi-glueball, torelon pair) states and it provides estimates 
of finite lattice spacing infinite-volume results to be used in a safe continuum 
extrapolation. Another way of improving the measurements is systematically 
studying the operators employed to create and annihilate single glueball states 
in order to increase their overlap with the single glueballs to be measured and 
to decrease their overlap with all the other states. T his i ncludes the study of 
other smearing techniques, such as e.g. Teper fuzzing [147|. Conversely, the in- 
troduction of operators coupling most strongly to the unwanted multi-glueball 
or torelon-pair states allows for the calculation of mixing strengths in order 
to exclude all the unwanted states from further analysis. Another technically 
rather easy but computationally expensive way of improving the results is per- 
forming simulations on additional values of /3. This allows for a more detailed 
study of the continuum limit, especially in the case of non-perturbatively im- 
proved actions, where there is no clearcut information about the way how the 
continuum is approached. Naturally, one has to take care not to go too deep 
into the strong coupling region as there is not much information about the con- 
tinuum; performing simulations on two additional values of /3 corresponding to 
(os/ro)^ « 0.16,0.3 would probably be a good choice in our case. But keep in 
mind that for each additional lattice spacing the measurements of the torelon 
dispersion relation as well as of the static-quark potential have to be repeated. 

There is one error which is not cured by the measures recommended in the 
last paragraph, namely the effects on the scalar states, coming from the "scalar 
dip" . Morningstar and Peardon improve the situation adding a two-plaquette 
adjoint term with a negative coefficient which results in an approach to the con- 
tinuum on a trajectory always far away of the dangerous "dip" region [|l9|. In 
principle, the classically perfect action could be treated the same way: Extract 
all the operators with adjoint contributions present in the parametrisation, de- 
termine their sign and add another (non-linear) constraint to the fit, namely 
that the action of all these operators together corresponds to an adjoint oper- 
ator with a certain (negative) coupling. Because of the large freedom in the 
fit, the inclusion of this single criterion should not impair the quality of the 
parametrisation considerably. 

Table 3.12 compares the newly obtained continuum glueball masses with 
the result of the isotropic FP action as well as with results obtained by other 
collaborations. There is reasonable agreement between the different determi- 
nations of the 2++, 0++ , ^, 2 ^ glueball masses. The situation for the 
scalar O"*""*" and for the 1^ state is different, however. The scalar glueball mass 
comes out about 5 standard deviations below the average mass of the other 
collaborations if only the usual statistical errors are considered. Thinking of 
all the possible sources of errors stated above, which are (partly) also present 
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CoUab. 



ror7io++ 



rom2++ 



UKQCD [[07 1 
Teper [|l37| 
GFU (l3|] 
M&P @ 
Liu M 



4.05(16) 
4.35(11) 
4.33(10) 
4.21(15) 
4.23(22) 



5.84(18) 

6.18(21) 

6.04(18) 

5.85(8) 

5.85(23) 



FP action 
AICP 



4.12(21) 
3.64(10) 



[5.96(24)] 
6.15(15) 



Table 6.12: Comparison of the lowest-lying glueball masses in units of tq. Values 
in brackets denote masses obtained at a lattice spacing a = 0.10 fm and are not 
extrapolated to the continuum. 



Collab. 
Teper [ ]l37 | 

M&P mi 



5.94(68) 
6.33(13) 



7.86(50) 
6.50(51) 



8.42(78) 
7.55(11) 



romi+- 

7.84(62) 

7.18(11) 



FP action 
AICP 



[6.74(42)] 
6.46(18) 



6.90(38) 



[8.00(35)] 
8.75(26) 



[7.93(78)] 
9.45(71) 



Table 6.13: Comparison of glueball masses in units of tq. Values in brackets 
denote masses obtained at a lattice spacing a = 0.10 fm and are not extrapolated 
to the continuum. 
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in the analyses of the other groups, we tend to explain this discrepancy with 
underestimated (or disregarded) errors. Concerning the 1"' glueball, looking 
at the masses determined on the lattices (T^ in Table |6^ , we observe that the 
masses determined on the coarser (and larger) lattices at (3 = 3.15, 3.30 coincide 
with the (continuum) results obtained by other groups and do not even show a 
significant discrepancy, despite the considerable difference between the lattice 
spacings. The mass determined on the fine (and small) lattice at (3 = 3.50 
however, is much higher. This tendency is even amplified by the continuum ex- 
trapolation. We suspect that this state exhibits strong finite-size effects pushing 
up its mass in small volumes. 




Figure 6.22: Lattice results for the scalar (0++) glueball mass, together with 
results of continuum extrapolations obtained from simulations employing dif- 
ferent actions. The different continuum values for the Wilson action stem from 
different groups. 



Figures 3.22 and [6.23 compare our measurements of the A^^ , and 



T2"^ states to measurements obtained by other groups, using different actions, 
as well as to the isotropic FP action results. Concerning the scalar glueball, 
the artifacts at moderate lattice spacings corresponding to (as/rg)^ < 0.15 
{as < 0.19 fm) seem to be rather smaller for the isotropic FP as well as for the 
anisotropic classically perfect action compared to the tadpole and treelevel im- 
proved action (M&P) and certainly much smaller than the artifacts of the Wilson 
action. The mass obtained on the coarsest lattice at f3 — 3.15 (corresponding 
to as ~ 0.24 fm) using the classically perfect anisotropic action (AICP), how- 
ever, exhibits large cut-off effects of about 30% (compared to about 20% for the 
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Figure 6.23: Lattice results for the tensor (2++) glucball mass, together with re- 
sults of continuum extrapolations obtained from simulations employing different 
actions. The continuum extrapolations of the E^^ and T^~^ representations are 
averaged in order to get the mass of the single continuum glueball. The different 
continuum values for the Wilson action stem from different groups. 
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M&P action), probably due to the scalar dip. Concerning the tensor glucball, 
the picture is not clear, mainly due to the considerable statistical errors of the 
measurements, except of the ones obtained by using the tadpole and treelevel 
improved anisotropic action. 

These considerations show again, that in order to obtain fool-proof results 
measuring glueball masses (and other quantities) on the lattice, one has to 
make sure to measure everything that is necessary, taking into account a certain 
number of lattices with different couplings /3, lattices with different volumes, 
measuring the renormalised anisotropy for all the values of /? considered and 
measuring additional operators having large overlaps with states that could 
mix. 

Being the first application of the FP action technique to anisotropic gauge 
theory, it goes beyond the scope of this work to check all the possible sources 
of errors in the glueball mass determinations. Observing the difficulties arising 
with the determination of the glueball spectrum already in pure gauge theory, 
one can estimate that reliable determinations in the context of more and more 
unquenched Lattice QCD (including lighter and lighter dynamical fermions un- 
til, finally, one ends up with quarks having their natural mass) will be a very 
tough job. However, this will yield information about the fate of glueballs in 
the presence of quarks and possibly solve the big puzzle concerning the presence 
(or absence) of glueballs in nature. It will thus be really worth the big effort. 



Chapter 7 



Repeating the Spatial 
Blocking Step 



As noted in Chapter |j, the spatial blocking step used to obtain a classically 
perfect ^ = 2 action may be repeated straightforwardly to obtain actions with 
higher anisotropies. Using another spatial scale 2 blocking step we create a 
classically perfect ^ = 4 gauge action and parametrise it the same way as the 
^ = 2 action. Furthermore, we measure the renormalised anisotropy and show 
that it is indeed feasible to construct such actions to be used in MC simulations 
for heavy states such as glueballs. The classically perfect = 4 action is not 
yet examined thoroughly, however the parametrisation is ready (see Appendix 



D.6) and the analyses of the ^ = 2 action (see Chapter u) may be repeated for 



7.1 Perfect Actions for Larger Anisotropies 

Each application of the spatial blocking step (slightly) renormalises the aniso- 
tropy of the action; the renormalised anisotropies of the final actions thus have 
to be measured in the end. However, this is necessary anyhow if one is interested 
in comparing results (other than e.g. mass ratios, where the renormalisation of 
the anisotropy cancels out) to other collaborations or to the experiment. 

In order to be able to repeat the spatial blocking step getting from ^ = 2 to 
^ = 4 we construct a ^ = 2 action which is valid on minimised ^ = 4 configura- 
tions including into the fit configurations at different values of /3 covering a large 
range. The S, = 2 action presented in Chapter ^ is not suitable for this task 
as it is dedicated to be used solely on largely fluctuating configurations around 
(3 — 3.0. The construction and the parameters of the intermediate £, = 2 action 



are described in Appendix D.5 



Once the fine action to be used on the r.h.s. of the renormalisation group 



equation, eq. 1.5, is ready, we may proceed the same way as for the ^ = 2 action. 
We perform the non-linear fit on 2 ^ad-hoc = 6 configurations each at /3ad-hoc = 
4.0, 3.5, 3.0, 2.5, 2.0. It turns out that in order to obtain parametrisations that 



are free of dangerous "traps" in the u-w plane (see Section 1.4.5) we have to 
include the condition A{u,w) > at ten points {u,w) into the non-linear fit. 
Studying the values of as well as the linear behaviour of the parametrisations 
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we decide to use a set with inax(fc + Z)sp = 3, niax(fc + /)tni — 2. The parameters 
are given in Appendix D.6, the behaviour of the hncar parameters is displayed 



in Figure 1.5 



7.2 The Renormalised Anisotropy 

To check whether the construction of the classicahy perfect ^ = 4 action worked 
well and whether the parametrisation reproduces the input anisotropy, we mea- 
sure the renormalised anisotropy at one value oi (5 — 3.0, using the torelon 



dispersion relation, following the method described in Sections 2.1 and 3.2, The 
simulation parameters are collected in Table [7.1| . 





5^ X L X T 


# sweeps / measurements 


3.0 


8^ X 4 X 32 


54000 / 10800 



Table 7.1: Run parameters for the torelon measurements using the ^ = 4 perfect 
action. The lattice extension in torelon direction L, the extension in the two 
transversal spatial directions S as well as the temporal extension L are given. 



Due to the (expected) coarser spatial lattice spacing we do not need as many 
smearing steps as for the ^ = 2 action and thus perform our measurements on 
configurations of smearing levels 2, 4, 6, 8, 10 keeping As = 0.1 fixed. 



Figure |7.1| displays the dispersion relation obtained, in Table G.IC the mea- 
sured energies E{p'^) determined using variational methods (see Appendix ^ 
depending on the lattice momentum are given, together with the number of 
operators used in the variational method, the fit ranges and the values of per 
degree of freedom, x^/-^df- The renormalised anisotropy as well as the torelon 
mass are determin ed in exactly the same way as for the ^ = 2 action. The results 
are listed in Table |7.2| . Furthermore, we may again evaluate an estimate of the 
lattice scale using the string picture relation, eq. yielding y/aus — 0.532(16), 
and thus ro/ag — 2.241(86) which corresponds to Os = 0.22 fm, at — 0.060 fm. 
Having a torelon of length Lag ~ 0.89 fm the error of the scale estimate is 



expected to be of about 5-10% (see Table 6.4) 





fit range 




Tot a* 




3.0 


1..5 


3.71(8) 


0.235(13) 


0.82 



Table 7.2: Results of the torelon simulations using the = 4" perfect action. 
The fit range in is given in units of (27r/iS')^. 



Concerning the application of the spatial blocking method presented in this 
work to yield classically perfect gauge actions for higher anisotropics, we may 
conclude that the construction including the blocking as well as the parametri- 
sation works in exactly the same way as the construction of the ^ = 2 action 
to which the main body of this work is dedicated. The resulting action shows 
a renormalisation of the anisotropy of about 7% at the value oi (3 = 3.00 where 
this has been measured. Calculations of physical observables such as the static 



104 



Chapter 7. Properties of the Actions 





0.5 
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Figure 7.1: Torclon dispersion relation for /? = 3.0. The straight line is the 
correlated fit to E'^{p) = m\ in the range = 1..5. 

quark-antiquark potential or glueball masses using the ^ = 4 action are still 
absent. 



Chapter 8 

Properties of the Perfect 
Anisotropic Actions 

In this Chapter, we examine the autocorrelation times and the average change 
of the gauge hnks for different updates in the MC simulation for anisotropics 
^ — 1,2,4. Furthermore, the computational overhead of the parametrisation of 
the perfect action is compared to the standard Wilson action as well as to the 
widely-used Symanzik and tadpole improved anisotropic action. 



8.1 Updating Algorithms 

There are different methods of updating Monte Carlo (MC) configurations. The 
two methods we use for our rather complex parametrised gauge actions are the 
plain local Metropolis algorithm as well as pseudo over-relaxation (POR), both 
described in the next two sections. 

In order to define the run parameters of MC simulations one should know 
about the properties of the updates, in order to have configurations that have 
small autocorrelation, to be sure that the algorithm does not introduce any bias 
and to estimate the optimal MC time between measurements of observables. 



8.1.1 Local Metropolis Update 

A very general method for updating MC ensembles satisfying detailed balance 



has been introduced by Metropolis et al. |148|. We update link by link, first 
proposing a new link U' e SU(3). The change [/ ^ C/' is accepted with prob- 
ability 1 if the action of the proposed configuration is smaller than the initial 
action, with probability e"^"^ if the action increases with the update. In our 
updating procedure we use proposed links V that are products of the initial 
link \J and random SU(3) matrices whose fluctuations are chosen such that the 
acceptance rate of the Metropolis sweep stays between 0.3 and 0.7. 

There is the possibility of performing multi-hit updates, i.e. using "coarser" 
trial matrices that lead to larger changes and compensating for the lower accep- 
tance by trying several proposed links in a row. 
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8.1.2 Pseudo Over-Relaxation Update 

The method of over-relaxation bases on the fact that for simple actions (as 
e.g. the Wilson gauge action) it is possible to perform large changes in the 
update of a link without changing the action 5, i.e. to perform micro-canonical 
updates. For SU(3) gauge theory this may be done executing such updates 
subsequently on different SU(2) subgroups. 

For our complex parametrisation, exact micro-canonical steps are not easy 
to implement. That is why we propose new links that would not change the 
corresponding part of the anisotropic Wilson action and correct for the change 
in the perfect action by adding an accept- /reject step, i.e. our pseudo over- 
relaxation (POR) algorithm is a variant of the Metropolis algorithm presented 
in the last Section, where the generation of the proposed links is altered. How- 
ever, one has to be careful about this procedure because all the configurations 
that are reached by subsequent POR updates lie on the same -4wiison — const, 
hyperplane in configuration space. In order to keep ergodicity, POR steps have 
to be mixed with another kind of (ergodic) updates, e.g. with Metropolis steps. 
Additionally, due to the similarities of the Wilson action and the FP action, dur- 
ing the thermalisation, when a certain action level shall be reached, it may be 
favourable to use rather less POR steps than later on, generating configurations 
at the same coupling [3. 



8.2 The Autocorrelation 

The autocorrelation is a measure of correlation between an observable with 
the very same observable at a later (MC) time. Near phase transitions, the 
autocorrelations in MC simulations of observables are proportional to a power 
n of the lattice correlation length : 

~ {e'-'T, (8.1) 

where theoretically n = 2 for local stochastic (Metropolis) and n = for cluster 



or over-relaxation updates. Eq. 8.1 shows as well that one has to expect crit- 
ical slowing down, i.e. a dramatic increase of the autocorrelation time, if one 
simulates near a phase transition with divergent correlation length ^ using local 
stochastic updates only. 

On the anisotropic lattice there are two different correlation lengths (in lat- 
tice units) ^^^'^ in spatial and ^J*** in temporal direction, where = at/as — 
(where this latter ^ denotes the anisotropy). 

Operators that couple strongly to ^J**', i.e. operators living in the spatial do- 
main, are thus expected to have autocorrelation times rising with the anisotropy 
for simulations that are performed not too far away from phase transitions. 

8.2.1 The Integrated Autocorrelation Time 

Let us first introduce the (unnormalised) autocorrelation function of an operator 
O, Cooit): 



Coo{t) = {0{T)0{T + t))-{0)\ 



(8.2) 
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As well, we define the normalised autocorrelation function p{t): 

Typically, p{t) decays exponentially, p{t) ~ exp(— |t|/r) for \t\ large enough 
and we thus define the exponential autocorrelation time Tcxp 

Toxp.o = lim , I ^ — 7-T7, (8.4) 

which is the "lifetime" of the correlation. This definition is vivid, however what 
is more important for the statistics of MC runs is the integrated autocorrela- 
tion time Tint which is describing the loss of statistics due to the subsequent 
configurations not being completely independent: 

var(O)«-(2ri„t,o)^^(0), (8.5) 
n 

i.e. the true variance of the estimated observable O is a factor 2ri„t,o larger than 
if the configurations were completely independent of each other. 

There are different possibilities of estimating Tint • Firstly, there is an inequa- 
tion between the two autocorrelation times [149-151| which holds for algorithms 
that satisfy detailed balance: 



The exponential autocorrelation time Tcxp can be obtained performing fits 
to p{r). For a direct estimate of Ti nt on may use the same procedure as used by 
the QCD-TARO collaboration in p5^ ], where 

N 

^int = ^p(0)+^p(t)^^, (8.7) 

t=l 

where N is determined such that Tint is maximised, at the same time iV < 10 % 

of the total sample. 

Another method has been proposed by Sokal in |15C, 151|: 

1 ^ R{t) 



with 



n-\t\ 
\t\ ^ 



R{t) = zr^,Y.^o^-om+\t\-o), (8.9) 



where n is the number of measurements performed (extent of "MC time" ) . The 
parameter m defines a "window" in which the measurements are taken into 
account to calculate the integrated autocorrelation time. The corresponding 
variance of the estimate is 

, , 2(2m + l) . 
var(Tint) « ^ '-tI, (8.10) 
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for r ^ TO <C 71. There are several possibilities choosing to; generally, large m 
leads to a small bias of the estimated Tint, however, the variance is growing with 
increasing to. One possibility is to fix to = CTjnt with a factor c depending on 
the behaviour of p{t) (c « 4 for purely exponential decay, c « 6 for asymptotic 
or pre-asymptotic decays or even c « 10 20 for exceptionally slow decays). 
However, analogous to the procedure described before, we choose to such that 
the estimate for Tint is maximised. It turns out, that there is always a long 
plateau of maximal Tint and that the bias corresponding to the to chosen is small. 
That is why we display the resulting values of the integrated autocorrelation 
time together with their standard deviation. 

In practice, the effect of autocorrelations of the operators between measure- 
ments close in the MC trajectory is taken into account by binning subsequent 
measurements into blocks of length Tmc > Tint , where Tint should be the largest 
integrated autocorrelation time of the operators measured in the simulation. 
These binned measurements are independent from each other and yield the cor- 
rect variance of the observables, if they are used e.g. in a bootstrap procedure. 

8.3 The Change in the Update 

A measure for the change of the configuration in a sweep is the following 

(At7),pd - E - Tr(f/,[7;t)], (8.11) 

^ ^— 1 

which is the average change per gauge link in one update; A^upd is the number of 
update attempts, U and U' denote the initial and the updated link respectively 
(of course U' = U for rejected steps). It is this change that is expected to be 
considerably larger for (pseudo) over-relaxation updates compared to standard 
Metropolis steps. 

The measure defined above may also be used as a measure for autocorrela- 
tions of the most local objects, i.e. the gauge links itself, if it is used to compare 
an initial configuration ([/) to the configurations at MC time t = 1,2,3,... 
(f/'). 



8.4 Results 

To study the autocorrelation of the two update algorithms used, local Metropolis 
and POR, we perform MC runs (« 100000 sweeps each) on 6^ x 6^ lattices with 
Os/t'o ~ 0.3 using different combinations of Metropolis (M) and POR steps: 
1 M, 1 M + 1 POR, 1 M + 2 POR, 1 M -)- 4 POR; there is no sense in 
using pure POR updates as they are not ergodic and must not be used alone 
in actual simulations. The gauge operators measured are the simple plaquette, 
either purely spatial [p^s) or going into temporal direction (pst), as well as the 
1x2 rectangle: purely spatial (rgg), extension 1 in temporal direction (rgti) and 
extension 2 in temporal direction (rst2)- 

First, the simulations are performed for the £, — 2 perfect action, the results 



are listed in Table 8.1. Firstly, it turns out that the statistics necessary is indeed 



very large {n w lOOOTjnt): in the case of pure Metropolis updates where this 
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condition is not fulfilled for all the operators, the statistical errors are of about 
30%. Secondly, one may notice that the inequation between the exponential and 
the integrated autocorrelation time, eq. ^.6| , holds — thus, if one is interested in 
an upper bound of Tint, e.g. to choose the bin size in a simulation, it is enough to 
estimate Tcxp (for the measured operators) which is rather cheap as this can be 
done with much less statistics up to errors of about 10-20% which normally do 
not matter. The combination of standard Metropolis and POR updates which 
turns out to be optimal for the decorrelation of small operators, using the ^ = 2 
action, is 1:1. 



op. 


1 M 


1 M H 


- 1 POR 


1 M H 


h 2 POR 


1 M H 


h 4 POR 


^cxp ^int 




Tint 


'^exp 


Tint 


'^cxp 


Tint 


Pss 


395 125(28) 


45 


31(4) 


59 


38(5) 


67 


42(5) 


Pst 


227 71(16) 


36 


32(4) 


52 


49(7) 


66 


61(8) 




392 172(39) 


54 


34(5) 


64 


38(5) 


72 


39(4) 


rsti 


279 73(17) 


39 


25(3) 


54 


37(5) 


67 


45(6) 


rst2 


280 126(29) 


44 


30(4) 


58 


40(5) 


70 


46(5) 



Table 8.1: Autocorrelation times for simple gauge operators using different up- 
dates with the ^ = 2 perfect action. The integrated autocorrelation time is 
calculated using the method described by Sokal, the parameter m chosen cor- 
responds to the maximum of the estimated values of Tint. 



e-1 e=2 e=4 



op. 


^cxp 


Tint 


^exp 


Tint 


'^cxp 


Tint 


Pss 


17 


15(6) 


45 


31(4) 


17 


8(2) 


Pst 


18 


13(3) 


36 


32(4) 


17 


18(4) 




21 


13(6) 


54 


34(5) 


19 


10(4) 


Tstl 


21 


10(5) 


39 


25(3) 


18 


12(2) 


rst2 


21 


11(4) 


44 


30(4) 


18 


13(3) 



Table 8.2: Autocorrelation times for simple gauge operators using alternating 
Metropolis and POR updates for actions with different anisotropies ^ = 1,2, 4. 
The integrated autocorrelation time is calculated using the method described by 
Sokal, the parameter m chosen corresponds to the maximum of the estimated 
values of Tint. 



The integrated autocorrelation time listed is calculated using the method 
proposed by Sokal, the method of QCD-TARO yields results in excellent agree- 
ment, however it seems as if the Sokal method does a better job in distinguish- 
ing between small contributions to the integrated autocorrelation time and the 
noise, the method is more stable and the window used to determine Tint has 
a more reasonable size. Furthermore, because most of the values of the QCD- 
TARO method lie below the values of the Sokal method (however, within one 
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standard deviation of the Sokal results), it is more conservative to use the Sokal 
resuhs. 

FinaUy, we compare the autocorrelation times for the different anisotropic 
actions looking at the 1 Metropolis + 1 FOR update, which turned out to be 
the most efficient for ^ = 2, see Table 3.2 and Figure |8.2| . Rather surprisingly, 
it turns out that the average autocorrelation times (of these small operators) 
are by a factor k, 2.5 larger for the = 2 action compared to the isotropic and 
the ^ = 4 action. It is not observed, that the spatial operators Pss and r^s have 
larger autocorrelation times than the temporal operators, which is probably due 
to being far from the critical region near the deconfining phase transition. 




# POR steps 



Figure 8.1: Integrated autocorrelation times of different operators using different 
update algorithms (1 Metropolis step + n POR steps, n = 0, 1, 2, 4) for the 
— 2 action. The autocorrelation times displayed by the open symbols as well 
as the errors have been determined by the method of Sokal, the crosses display 
the corresponding results obtained with the method of QCD-TARO. 



Table S.3 lists the average change per link for Metropolis and POR updates 
using the perfect ^ = 2 action at different values of the coupling (3. It turns out 
that the average change of POR updates is about a factor of 10 larger than in 
the case of Metropolis updates; the normalised values are changing slowly with 
the lattice scale. 

Comparing configurations separated by N MC updates using the measure 
eq. 8.11, we obtain the following (obvious) behaviour: 

(AC/)„pd(Af)-3-ao-(l-(7)^, (8.12) 

where (AC/)upd = 3 indicates that the configurations are completely indepen- 
dent on the level of gauge links. The factor q determines the speed of the decor- 
relation. The results, obtained using only 20 starting configurations for each 
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Figure 8.2: Integrated autocorrelation times of different operators using 1 
Metropolis + 1 OR updates for the ^ — 1,2,4 actions. The autocorrelation 
times displayed by the open symbols as well as the errors have been determined 
by the method of Sokal, the crosses display the corresponding results obtained 
with the method of QCD-TARO. 





/? 


(□> 


1 / A 7- T-\ Metropolis 
■{Syv^^ Jupd 


{□) V^"^ )upd 


1 


2.86 


1.917 


0.0652 


0.472 




3.05 


1.846 


0.0527 


0.486 




3.40 


1.735 


0.0660 


0.506 


2 


3.00 


1.885 


0.0546 


0.613 




3.15 


1.702 


0.0552 


0.584 




3.30 


1.584 


0.0568 


0.577 




3.50 


1.487 


0.0572 


0.574 


4 


3.00 


1.750 


0.0400 


0.328 



Table 8.3: The average change per link for Metropolis and pseudo over- 
relaxation updates for the classically perfect actions with anisotropies ^ = 1,2,4, 
normalised by the plaquette expectation value (□) = {Nc — Re Tr L/p'). 
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anisotropy and update and performing 100 sweeps each are displayed in Table 
3.4. It turns out that these results are already completely stable. Obviously, 



these numbers disclose no secret and it is clear that all the effects apparent in the 



autocorrelation times, Tables 8.1, 8.2, stem from larger range correlations. Note 
as well, that the measure, eq. 3.11 is not gauge invariant, i.e. the gauge trans- 
form of a configuration may show to be completely independent, (A?7)upd = 3, 
of the initial configuration. 





1 M 


1 POR 


1 M + 1 POR 


1 M + 2 POR 


1 M + 4 POR 




q 


q 


q 


q 


q 


1 


0.015 


0.183 


0.102 


0.128 


0.152 


2 


0.012 


0.174 


0.094 


0.121 


0.143 


4 


0.008 


0.064 


0.037 


0.046 


0.054 



Table 8.4: The speed of the link decorrelation q (see text) for different updates 
at different anisotropics. 



8.5 The Computational Overhead 

One very important quantity of an improved action is its computational over- 
head compared to the standard (i.e. Wilson) action. 

To measure the effective gain or loss of an action is a very difficult task that 
depends very strongly on the application (spectroscopy, thermodynamics etc.), 
on the lattice sizes and lattice spacings used and as well on the architecture of 
the computers (scalar, vector, parallel and a lot of more subtle details) used. 

Here, we will thus simply give raw ratios between the anisotropic perfect 
action, the widely used Symanzik and tadpole improved action and the Wilson 
action. It is clear that our parametrisation which is quite sophisticated has 
a large "raw" overhead, mainly due to the large number of operators that are 
affected if one link is changed in the MC update. Presently, our parametrisation 
is probably the most complex one that may still be considered to be utilised in 
numerical simulations of pure gauge theory. Although, in principle, the fit to the 
exact perfect action could be done much better, e.g. employing double smearing 
or including diagonal links into the smearing, such a procedure would surely 
make it unsuitable for MC simulations. Of course, if the gauge action is used 
together with (quenched or even dynamical) fermions, the computational cost 
of the gauge part is rather negligible, even better parametrisations of perfect 
actions might be imaginable in this case. 



Table 8.5 compares the perfect anisotropic action to the Symanzik/tadpole- 
improved action of and the Wilson action on lattices of size S'* on an Alpha 
workstation. The numbers refer to the pure Metropolis update of each link of 
the configuration ("1 sweep"). 



8.5. The Computational Overhead 
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action 


^sweep [s] 


^act /^Wilson 


Wilson 


0.095 


1 


Symanzik 


0.53 


5.6 


Perfect 


6.3 


66 



Table 8.5: The updating time (1 sweep on a 8"* lattice) and the computational 
overhead of the perfect anisotropic action, compared to the Symanzik and tad- 
pole improved action B3] and the standard Wilson action. 



Chapter 9 

Conclusions and Prospects 



In this work, wc have presented the construction and paramctrisation of a clas- 
sically perfect anisotropic SU(3) gauge action based on the Fixed Point Ac- 
tion technique. The recently parametrised isotropic FP action as well as its 
paramctrisation ansatz using mixed polynomials of plaqucttcs built from single 
gauge links as well as from smeared links, have been the starting point. Perform- 
ing one and two purely spatial blocking (RG) steps, respectively, starting from 
the isotropic action, and adapting the paramctrisation to distinguish between 
spatial and temporal directions, we have obtained parametrised classically per- 
fect actions with anisotropies ^ = 2 and 4. That the method works has been 
checked as well in the context of scalar fields and in the quadratic approxima- 
tion for SU(3) gauge theory. The ^ = 2 action has been tested extensively in 
measurements of the torelon dispersion relation, of the static quark- antiquark 
potential, of the deconfining phase transition as well as of the low-lying glue- 
ball spectrum of pure gauge theory. For the ^ = 4 action we have checked 
that the renormalisation of the anisotropy remains small and that it exhibits 
no pathologies in Monte Carlo simulations. For both actions (and addition- 
ally the isotropic FP action) autocorrelation times have been measured and the 
computational overhead, compared to the widely-used anisotropic tadpole and 
tree-level improved action as well as to the standard Wilson gauge action, has 
been determined. 

The results of the torelon measurements show that the renormalisation of 

the anisotropy (due to quantum corrections and paramctrisation artifacts) is 
small (below 10% for ^ = 2 and 4). The rotational invariance of the action has 
been examined by measuring the static quark-antiquark potential, including 
separations corresponding to a large set of lattice vectors, again with good 
results. The glueball measurements, including some rather heavy states, confirm 
that the use of anisotropic lattices facilitates spectroscopy when heavy states 
arc present. Compared to the isotropic simulations, with the same amount 
of computational work it is possible to resolve states with considerably larger 
energy, allowing reliable continuum extrapolations from larger ranges of the 
lattice spacing. The mass of the scalar ghieball with quantum numbers J^*^ = 
0++ is measured to be considerably lighter (1440(40) MeV) compared to masses 
obtained by other groups (around 1670 MeV); furthermore it shows large cut-off 
effects when it is measured on a lattice with spatial lattice spacing « 0.25 fm, 
which could be caused by a large sensitivity of the action to a critical point 
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in the fundamental-adjoint coupling plane, due to adjoint terms in the action. 
However, the sources of possible errors in measurements of the glueball spectrum 
are manifold and there is no ultimate statement about the goodness of the 
parametrised classically perfect actions. 

During our studies, we have noticed that the examination of scaling proper- 
ties of lattice gauge actions is a very delicate problem: Quantities which can be 
reliably measured and which are not very sensitive to systematic factors (such 
as the volume, other states mixing with the observed ones, the method used 
for extracting the mass etc.), such as the critical temperature or the hadronic 
scale ro (for moderate lattice spacings) generally exhibit rather small cut-off 
effects and thus demand very large statistics, if the differences between the ac- 
tions should be explicitly identified. On the other hand, some glueball states 
show large cut-off effects which make the glueball spectrum another interesting 
object to study scaling violations. As pointed out in this work, however, there 
are a lot of systematic factors making the extraction of the pure lattice artifacts 
difficult. The task of comparing different actions is even more difficult if the 
measurements and analyses are performed independently by different groups 
which introduces additional systematic discrepancies and ambiguities. 

Let us therefore conclude that despite the rich paramctrisation, the corre- 
spondence to the classically perfect FP action and the approximate scale in- 
variance of classical instanton solutions there is no conclusive evidence of the 
parametrised classically perfect actions (isotropic as well as anisotropic) behav- 
ing significantly better than other improved actions (such as the tadpole and 
tree-level improved anisotropic gauge action) — however, there is also no ev- 
idence for the converse. What is clear and what one has to consider if one 
is about to plan a project where the update of the gauge configurations takes 
a considerable part of the computer time, is the considerable overhead of the 
paramctrisation compared to other actions. 

For the future, very accurate scaling tests comparing the classically perfect 
actions to other improved gauge actions are desirable. On one hand, these 
can include large statistics measurements of the critical temperature and the 
hadronic scale, on the other hand, one might perform very careful simulations 
determining the glueball spectrum, including all the limits to be taken and 
excluding all known sources of errors. These simulations should all be performed 
by the same collaboration in order to minimise systematic discrepancies coming 
from the analyses. 

Another promising plan would be to gain control about the adjoint operators 
present in the parametrisation of the classically perfect actions in order to cir- 
cumvent the critical point in the adjoint coupling plane such that the influence 
of the scalar dip is minimised. 
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Appendix A 

The FP Action for Scalar 
Field Theory 



A simple way to get the lattice Fixed Point Action ior d — 2 free scalar theory 
is performing one single renormalisation group transformation (RGT) with an 
infinitely large scale factor starting in the continuum, a procedure called blocking 
out of continuum. The lattice is laid over the continuum field which is averaged 
over certain regions around the future lattice sites to obtain the lattice fields. 
The resulting FP action depends on the form of the averaging (corresponding 
to the blocking in finite scale RGTs) but the procedure presented here will be 
completely general. 

We use the following RGT, accomplishing the task of averaging the contin- 
uum field ip{x) to get to the lattice field </>„, where the form of the averaging is 
described by the function 11 (a;) 



e-^' - W = JVipix)e ' ' "~ , (A.l) 

the parameter k is again allowing the lattice field to fluctuate around the exact 
average of the continuum field (largely for small k). The averaging is normalised 
to / cPxII{x) = 1. To account for an anisotropy, ai ^ 02, we perform the 
calculation using dimensionful quantities, the symbol a is a shorthand notation 
for (ai, 02). 

Our goal is to obtain the FP propagator D^^ on the lattice, we thus calculate 
the lattice two-point function D{n — n') = {(f>n4'n') in this theory: 

1 f ~^Y. p{r)<P^4>^+. ^ 

D{n-n')^- V4>e 0„ (A.2) 



inserting eq. A.l we get 



1 ^d?xa^ip(x)a^Mx)-i^T.('l>n-^<fixn{na-x)if>(x))'^ 

- / V(j)Vipe " 0„ (A.3) 



Z ^ 

The blocking (averaging) kernel is now Fourier transformed in the continuum: 
K{n) = / (fxn{na ~ x)if{x) : (A.4) 
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Appendix A. The FP Action for Scalar Field Theory 



n(na-a;) =y e*f(""-^)n(p), (A.5) 
p 

^{x) = j e"^^ip{k), (A.6) 



thus 



K{n) = Jd^xJ J e*f""e'('=-f)^n(p)(p(fc) = J e''=""n(fc)(p(fc). (A.7) 

p k k 

Now we split up the continuum momentum variables according to: 

ki = qi + —li, k2 = q2 + —l2 (A.8) 
ai a2 

with qi, q2 being lattice momenta in the respective Brillouin zones: 

(A.9) 





-) 


. 92 e 1 


'-1. 


- 




aij 




02/ 



In the following, we will use the notation k = q + for convenience, where 
k = {ki,k2). Inserting this replacement, we get: 



[ dk jf,^^ _ [ dqx f dq2 i(q^ + ^i^)niaij(q2 + ^h)n2a2 

-00 l=-oc_^ _^ 

/oo 
1— 



using the shorthand notation 

TT TV 

ai a2 

dqi f dq2 



(A.ll) 

27r ./ 27r ^ 



9 



Using this, our two-point function reads: 



n q l = — oo f 



(A.12) 
We replace by 

m 



/°° Stt Itz 

^iqna n(g + —l)^{q + — /) (A.13) 
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and get: 



1 



: Jd^a; a^i^(a;)a^i^(a;)-K;^ (0* )^ 



(A.14) 



Performing the integration over 0* , the terms hnear in (/)* disappear. 

If we integrate both sides of eq. A.l with respect to (j)n we get Z on the l.h.s. 
and Z^Z^ on the r.h.s., with 



J d'^xd^ip{x)dfiip{x) 



and 



Using this notation we may write eq. A.14 as: 



Z 



(A.15) 
(A.16) 

(A.17) 
(A.18) 



The first term is just 



Z K K 



Ql°2 / iq{n-n')a 

^ t 



(A.19) 



the second term reads 



Zfj, 

'z 



^ n(fc)^(fc)nt(fc)^t(fc) 



e 2 



(A.20) 



after inserting the content of 6(g). Using eq. A.15, the definition of Zip, we get 

Vipipiq + —l)ipHq + 27r/)e-^ 9^^(:.)a^v W ^ z . Mq + 2ttI), (A.21) 
a 

and thus the second term is just 



Mn-n')a y n(g + + ^;)nt (g + — ; (A.22) 

^ — ' a a a 

1— — 00 



together with the first term: 



D{n^n')= /e*«("-")° - + Tl{q + 2TTl)d{q + 2TTl)\l\q + 2ttI) 



(A.23) 
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Appendix A. The Quadratic Approximation for Gauge Fields 



Thus we get p{q): 



PiQ) 



= E 



U{q + —l)d{q+—l)Il\q + —l) 
a a a 



aia2 

K 



(A.24) 



Setting ai = a2 = 1 we get back the isotropic FP action as expected. 

We are interested in couplings f){q) on the lattice which are dimensionless 
quantities. In order to achieve this, we rescale the above formula by dividing it 
by aia2, this absorbs the product in the «;-term, also, the propagator d{k) will 
become dimensionless: 



'^(^^ a,a2 {qi + + iq2 + T Hpi + 27r/i)2 + ^(^2 + 271/2)^ ' 

(A.25) 

where pi, p2 are dimensionless momenta, G (— tt, tt), and ^ = |^ is the 
anisotropy parameter as defined above. 



Appendix B 

The Quadratic 
Approximation for Gauge 
Fields 



In this Appendix we describe the calculation of the couplings Pfi^{r) in the 
quadratic approximation to the classically perfect action, eq. 



5.ie 



^^"^ = + +0(^3). (B.l) 

^ n,r 

B.l The Anisotropic Recursion Relation 



At the quadratic level, the recursion relation corresponding to eqs. L5, O 
reads in standard matrix notation 



4.2 



-^J2^T{B\kBW{kB)mkB)} 



= min^^j (^^Tr{At(fc)Sp(fc)SA(ft)} 

\ k 

+ 7^ E Tr {(r(fcs) - B{kB)yE^iE{T{kB) - 5(^5))} J , (B.2) 

where Vb = j^V is the volume of the blocked lattice, B{kB) denotes the blocked 
gauge configuration, A^, are the gauge potentials of the configurations U, 



V^inB) = e'^- \ [/^ (n) = e'^- 



and 



/ 1 



(B.3) 



(B.4) 
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Appendix B. The Quadratic Approximation for Gauge Fields 



In the above formula p, k are expressed in units of and the blocking 
function uj{k) is a dimensionless quantity, while 



k'^^a^ and 



n 



u=0 _Jfdf 



(B.5) 



where is a common symbol for the gauge potentials A^, etc. 



B.2 The Calculation of the Couplings 



The blocking used in this work (see Section [4.2| ) corresponds to 



with 



- (1 + e^'O [co'^M. + 6ciT«(fc) + 12c2r(2)(fc) + Scsrl^Jik) 



(B. 



(B.7) 



The tensor uj^^ is fixed by the form of the blocking kernel, where up to this 
order the linear contributions r^'J to the fuzzy link operator have to be taken 
into account. They can be written as 



r(2)( 

r(3)r 



1 

6 
1 

24 
1 

48 



fc^/c* + Sf,^{6 - 7) 



k^ k^j 



7 + fc^fc; + k^kl) + (5^^(24 - 87 + 7^ - ^7 - ikt^kD 
k^kl{2A - 67 + 7' - ?? + 6fc^fc; - 27fc^fc; + 2{k^klf + 6fc,fc: 



27/c^/c* + 2 



+ (5^.(48 - 247 + 67" - 7''^ - 67fc^A:; + 272fc^fc; - 27(A:^fc, 
+ - 2e - 6r] - 2Tjk^k*^) 

The lattice momentum fc^ is defined by 



* \2 



- 1 



(B. 



(B.9) 



and the following definitions have been made for convenience: 



The four parameters of the RGT in eq. B.7 have been optimised in |60| and 
have the values: 



ci = 0.07, C2 = 0.016, C3 = 0.008, k = 8.8. 



(B.IO) 



B.2. The Calculation of the Couplings 
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From eq. B.2, using the relation eq 
the inverse of the couphngs Df^^{k) — p'^^{k) 



we obtain a recursion relation for 



(B.ll) 



One way to find the FP solution for the propagator is to start from the 
anisotropic Wilson propagator and to iterate the propagator to the FP. As we 
have to invert the coupling tensors p^^ and p'^^ we have to introduce a temporary 
gauge fixing. The gauge-fixed anisotropic Wilson action reads 



e\k4 



a- 



{E^\k^\'+e\kAn 



(B.12) 



At the end the gauge-fixing will be removed, which corresponds to the limit 
a "> oo. 

Using the gauge relation 



'-^iiu{k)k,y — 2fc^ 



(B.13) 



one can show that starting from the standard propagator, after an arbitrary 
number of iterations, the propagator D^j_y(k) assumes the form 



Df,,{k)^G^,{k)+af{k)k^kl 



(B.14) 



with G^y{k) and J{k) independent of a. It may now be shown |6(|| that in 
the limit a — )■ oo the function /(fc) (which is the only part differing from the 
isotropic case) drops out and the result for the couplings is 



Pt^y{k) 



G-l{k) 



G-^{k)k,-k*G-^ik) 
k*Gap {k)kp 



(B.15) 



To construct the anisotropic propagator p ^(q) numerically we perform 3 
recursions on a 32'' lattice using the k value obtained optimising the locality in 
the isotropic 8.8. We study the dependence of the number of recursions 

and do not discover any notable changes if we change it from 3 to 4 - the FP 
quadratic propagator is thus reached already after 3 recursions. The following 
steps, inverting p~^{q) to get the couplings in momentum space, removing the 
gauge fixing (analytically) and the subsequent Fourier transformation do not 
change from the isotropic case described in [^|60|. 

A priori there is no reason why the parameter k occurring in the recursion 
relation should be chosen to have the same value for the spatial and the temporal 
part of the relation, i.e. one may think about changing n to the more general 
form 



\ 



Kt / 



124 



Appendix B. The Norm of the Anisotropic APE Action 



There are two "natural" choices, either stick to the (5^,^ form of k and set Kt — 
Kg = 8.8 or change it to Kt — ^'^Kg = ■ 8.8 to compensate for the insertion 
of the two S- matrices in eq. B.ll. Performing numerical studies (in the full 
theory) it indeed turns out that the latter choice is the one yielding actions that 
are well parametrisable. 



Appendix C 

The Norm of the 
Anisotropic APE Action 



In this Appendix we derive the anisotropic normahsation condition which makes 
sure that the action goes over to the continuum form for very smooth fields. 
We first define for a general continuum gauge field a dimension 4 operator 

Rq: 

Ro = -lY.^Ti:Fl). (c.i) 

To work out the form of the conditions for our parametrisation of the ac- 
tion we use constant gauge non-Abelian gauge potentials, df^A^ = 0. In the 
continuum we have 

T^, = [A^,A,]. (C.2) 
The constant non-Abelian gauge-field is defined on the lattice as 

Af^ = ia^ ■ -cTp, (C.3) 

where Cfc, fc=l, 2, 3 are the Pauli matrices and — (cti + (T2 + cr3)/V3- 

To obtain the normalisation condition we expand our anisotropic gauge ac- 



tion, parametrised as described in Section 4.3, in powers of keeping the 
spatial indices k = 1,2,3 and the temporal index 4 separate. This yields two 
operators ~ 0{a'l) and Rq ^ 0{a^a^). Up to quartic order in as, at the 
action can be expressed in i?Q and R^: 

A - ^[2ci2pf,T+Pfo+Po?(l + 2cii + 2r;icn)] 

+ ^ [2%ciip^? + + (1 + 2ci3 -I- mcis + mcu)] , (C.4) 

where all parameters ry, c denote the constant (0*'' order) terms in the poly- 
nomials in x^; the higher orders do not contribute due to — > for very 
smooth fields. 
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Appendix C. APE Action Parameters 



The coefficient of Rq in eq. C.l is required to be 1 and thus a simple com- 
parison yields the two normalisation conditions for our parametrisation of the 
anisotropic gauge action: 



pZ+pZ + 2pZc,i-V2pl^c,2-V2pZci,m = \, (C.5) 

Po™+p'i™ + 2p^Tci3+p^Tci3772 + 2/0^1011^3 +P^Tci2774 = e (C.6) 

These two conditions ensure that the parametrised action goes over to the 
continuum form with correct anisotropy if the action is applied on smooth fields. 
They are required to be fulfilled exactly for all our parameter sets. 



Appendix D 

Action parameters 



In this Appendix we collect the parameters of the FP and classically perfect 
actions that have been used or constructed throughout this work. The isotropic 
FP action presented in Section D.l is the one that has been parametrised and 
examined closely in [Q-^ , results of these studies are collected briefly in Chapter 
0. The isotropic action Af^"^ in Section p.2| is an intermediate parametrisation 



valid on configurations typical for MC simulations that are minimised once, 
parametrised during the cascade process for the isotropic action. It is used in 
the spatial blocking (see Section 6.1) to minimise the coarse ^ = 2 configurations, 
constructed using the ad-hoc anisotropic action presented in Section D.3. The 
resu lting ^ = 2 perfect action for coarse configurations is presented in Section 

To repeat the blocking step, we need a ^ = 2 action which is valid for 
^ = 4 configurations minimised once in a purely spatial blocking step. Such an 
intermediate action has been constructed, see Section D.5. The resulting ^ = 4 
action is presented in Section D.6. 



D.l The Isotropic FP Action 

The parametrisation of the isotropic FP action is described in ||l|-^ as well as 
in Chapter |[ It describes the FP action well in the range of lattice spacings 
0.03 fm < a < 0.3 fm. 

The set of parameters consists of four non-linear parameters r]^^\ cf\ c^"^, 
cf' describing the asymmetrically smeared links W^'^^ and fourteen linear pa- 
rameters pki with < k + 1 < 4. This set approximates reasonably well the true 
FP action in the range of a given above. (For smaller fluctuations — occurring 
e.g. in the intermediate steps of the parametrisation — we use polynomials for 
ri{x) and Ci{x) up to forth order.) 

The optimal non-linear parameters are found to be 

77(") = -0.038445 , cf'> = 0.290643 , = -0.201505 , = 0.084679 . 
The linear parameters are coUected in Table pTl 
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Appendix D. APE Action Parameters 



Pkl 


I = 


I = 1 


I = 2 


/ = 3 


I = 4 


k = 







0.442827 


0.628828 


-0.677790 


0.176159 


k = 


1 


0.051944 


-0.918625 


1.064711 


-0.275300 




k = 


2 


0.864881 


-0.614357 


0.165320 






k = 


3 


-0.094366 


-0.020693 








k = 


4 


0.022283 











Table D.l: The linear parameters pki of the parametrised isotropic FP action. 



D.2 The Intermediate Isotropic Action A 
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The intermediate isotropic action has been parametrised during the construction 
of the isotropic FP action. It is supposed to be valid on configurations that are 
obtained by minimising configurations typical for MC simulations once. It fulfills 
the O(a^) Symanzik conditions (sec ||l|, |^) and uses fluctuation polynomials 
Xfj.{n) (see Section 4.3.1), it thus smoothly approaches the continuum (high /3) 
limit and is expected to interpolate between the rather coarse configurations 
mentioned above and the smooth limit. It is not intended to be used in MC 



simulations, its linear behaviour (see Section 1.4.5) is not checked, in addition 
the set of parameters is not as compact as might be possible. 

This action has been used for the (spatial) minimisation of coarse ^ 



configurations to describe the isotropic configurations on the r.h.s. of eq. 1.5 
The non-linear parameters describing order 3 polynomials are 





= 0.082, 


^(1) 


= 0.292353, 


^(2) 


= 0.115237, 


^(3) 


= 0.011456, 




= 0.282, 


^1 


= -0.302295, 




= -0.302079, 




= -0.052309, 




= 0.054, 




= 0.298882, 




= -0.081365, 




= -0.023762, 


^3 


= -0.201671, 




= 0.022406, 


C3 


= 0.004090, 


„(3) 
C3 


= 0.014886, 


L4 


= -0.008977, 


C-4 


= 0.245363, 


L4 


= 0.140016, 


t-4 


= 0.028783. 



The linear parameters are collected in Table 



Pkl 


I = 


/ = 1 


I = 2 


I = 3 


/ = 4 


k = 







0.629227 


-0.556304 


0.186662 


-0.010110 


k = 


1 


-0.368095 


0.852428 


-0.199034 


0.031614 




k = 


2 


0.389292 


-0.207378 


-0.010898 






k = 


3 


-0.054912 


0.039059 








k = 


4 


-0.000424 











Table D.2: The linear parameters pki of the parametrised intermediate isotropic 
if 



FP action Af^ 



D.3. "Ad-hoc" Anisotropic Actions 
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D.3 "Ad-hoc" Anisotropic Actions 

In the spatial blocking procedure described in Chapter |4| one needs coarse 
anisotropic = 2, 4, 6, . . . ) gauge configurations which are spatially min- 
imised leading to ^' = £,/2 configurations. As the goal of this step is to obtain 
a perfect anisotropic action with anisotropy ^ these coarse anisotropic config- 
urations have to be produced using some other action that is already present. 
This requirement might seem to endanger the whole ansatz, however it is not 
crucial how the coarse configurations exactly look like as the perfectness of the 
resulting coarse action comes from the perfect action on the fine configuration 



as well as from the exactness of the RG transformation (see Section 1.3). Still, 
we try to create coarse configurations that might look similar to future ensem- 
bles produced using the perfect anisotropic action on the coarse level and whose 
minimised configurations appear to have an anisotropy approximately ^' = ^/2 
(i.e. which are isotropic in the case of the construction of the £, = 2 action). 

In order to achieve this, we modify the isotropic FP action adding a term 
(^^ — 1 )pIq (where pf q denotes the simple temporal plaquette) . This modification 
turns the isotropic Wilson action into the Wilson action with bare anisotropy 
^ and is expected to work approximately also for our FP action. The main 
argument of using this coarse action and not e.g. the anisotropic Wilson action 
is that due to the spatial lattice spacing as being larger than the temporal 
spacing at, the 0{a'^) artifacts are also expected to be larger than the 0{af) 
effects. The modification described above should preserve the FP properties 
corresponding to and is thus expected to be considerably better than the 
naive anisotropic Wilson action. 

Using this ad- hoc modification of the isotropic FP action ( "ad-hoc" action) 
in MC simulations shows (as expected) that (at least for small anisotropics 
C ^5 5) its properties are not comparable with the ones of the isotropic FP 
action, however the generated ensembles resemble to the ones generated with 
true perfect anisotropic actions much more than ensembles generated by the 
Wilson action. 

The anisotropics ^ad-hoc that have to be used to generate coarse anisotropic 
configurations turning into minimised configurations with anisotropy ^' are 
Cad-hoc ~ 3.2 for C' = 1 and ^ad-hoc ~ 6 for C' — 2. However, this value of 
Cad-hoc varies considerably if /3 is varied — but as stated at the beginning of this 
section, the exact form of the coarse configurations is not essential. 



D.4 The ^ = 2 Perfect Action 



The £, — 2 perfect action uses the parametrisation described in Sections 4.3.1 



4.3.2. The number of non-zero asymmetry values rjf^^ is 4, the parameters cf^ 
{i = 1, . . . ,3) are splitted into 3 parameters depending on the contribution to 
the smeared plaquette. The linear parameters pki are non-zero for < A: + ^ < 4 
for spatial plaquettes and < fc -I- Z < 3 for temporal plaquettes. 
The non- linear parameters (constants in x^) have the values 
■qi = -0.866007, V2 = -0.884110, vs = 2.212499, r;* = 1.141177; 
cii = 0.399669, cia = 0.519037, cn = -0.071334, 
C21 = -0.076357, C22 = -0.031051, C23 = -0.282800, 
C31 = 0.032396, C32 = -0.015844, C33 = -0.046302. 
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The linear parameters are collected in Table D.3 
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Table D.3: The linear parameters of the ^ 
action. 



2 parametrised classically perfect 



D.5 The ^ = 2 Intermediate Action 

In order to be able to repeat the spatial blocking step constructing a ^ « 4 
action based on the ^ = 2 perfect action we need a parametrisation of the ^ = 2 
action which is valid on « 2) configurations that are obtained by spatially 
minimising coarse £, = i configurations once. To construct such an action, we 
perform a non-linear fit to the derivatives of 5 sets of two configurations each at 
/3 = 6, 10, 20, 50, 100. The non-linear parameters are chosen to be linear in the 
fluctuation parameter x^{n). Having four different parameters 77 and splitting 
up Ci into three parameters (as it is done for all anisotropic parametrisations) , 
this makes 20 non-linear parameters to be fitted whic h is q uite at the edge 
of what is still possible on our computers (see Section 1.4.4| ), that is why we 



restrict the total number of configurations to 10. Rough checks performed on 
a larger number of configurations, with an even larger number of parameters 
show however, that the resulting non-linear parameters are stable and describe 
the data accurately. 

The derivatives and action values of 5 sets of 10 configurations each are 
included in the linear fit (where the relative weight of the action values is chosen 
to be 1.9 • 10^2 for the configurations at /? = 50, 100 and 7.6 • lO"* at /? = 6, 
10, 20). A linear set where the parameters pki are non-zero ior < k + I < 3 
for spatial plaquettes and < fc -I- / < 4 for temporal plaquettes describes the 
full action very well concerning this data. Again, this parametrisation is not 
intended to be used in MC simulations, thus the linear behaviour of the action 
is not checked. 

The non- linear parameters (linear in x^) have the values 
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The linear parameters are collected in Table D.4 
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Table D.4: The linear parameters of the intermediate ^ = 2 parametrised clas- 
sically perfect action. 



D.6 The ^ = 4 Perfect Action 



The C = 4 perfect action uses the parametrisation described in Sections 4.3.1 



4.3.2 . The number of non-zero asymmetry values rlf^ is 4, the parameters cf^ 
(i = 1, . . . ,3) are splitted into 3 parameters depending on the contribution to 
the smeared plaquette. The linear parameters pki are non-zero for < A; + / < 3 
for spatial plaquettes and < fc -I- Z < 2 for temporal plaquettes. 
The non-linear parameters (constants in x^) have the values 
rii = -1.491457, 772 = -1.115141, 773 = 1.510985, 774 = 7.721347; 
cii = 2.014408, C12 = 0.128768, C13 = 0.162296, 
C21 = -0.915620, C22 = 0.134445, C23 = -0.013383, 
C31 = 1.166289, C32 = 0.061278, C33 = 0.000759, 



The linear parameters are collected in Table D.5 
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Table D.5: The linear parameters of the ^ = 4 parametrised classically perfect 
action. 



Appendix E 

Variational Techniques 



Normally, in a Monte Carlo simulation there are several operators that have 
an overlap with a state one would like to measure. In the case of glueballs, 
e.g., we have multiple shapes with combinations of different orientations having 
the same quantum numbers. Another method yielding several gauge operators 



with the same quantum numbers, is smearing (see Section 2.1.2). A priori 
it is unknown which operators have the largest contribution to a given state, 
it is therefore vital to have a method which takes into account the different 
operators according to their overlap with the desired state. In this Appendix 
we describe the variational techniques accomplishing this, used for all mass / 
energy determinations (torelons, static qq-potential, glueballs). 



E.l Obtaining Mass Estimates 



In a Monte Carlo simulation we measure the N x N correlation matrix of oper- 
ators Oa 



C^0{t) = {{)\0^{t)OlmQ) 



(E.l) 



where N is the number of operators having the same quantum numbers as 
the examined state. To determine the coefficients Va of the linear combination 
VaOa which has the largest overlap to the ground state relative to the 
excited states one has to minimise the effective mass given by 

-{v,C{h)v)- 



m{to,ti 



In 



/(h-to). 



(E.2) 



iv,Cito)v)^ 

The vector v is obtained by solving the generalised eigenvalue equation [ 153| , 154|] 

Cih)v = Xito,h)C{to)v, (E.3) 

where < to < ti. 

Assume first that only the lowest lying N states contribute to C(i), i.e.. 



N 



nf3 •> 



(E.4) 
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where Ei < E2 < ■ . ■ < E]\j are the energy levels in the given symmetry channel 
and ipna = (O|0a|n) is the "wave function" of the corresponding state. The 
solution of eq. |E.3|is given by the set of vectors {vn} dual to the wave functions, 



i.e., {vn,'>pm) = Snm- Multiplying eq. E.4 by Vn one obtains 



C{t)vn = e-^"V« = e-^"(*-*«)e-^"*°Vn = e-^"(*-*«)C(io)i;„ ■ (E.5) 



This gives A„(to,^i) = exp(— £'„(ti — to)) for the eigenvalues in eq. |E.3| . Of 
course, contributions from states with n > N and statistical fluctuations distort 



eq. E.4, therefore the stability of eq. E.3 is an important issue 



Observe that eq. |E.3| is well defined only for positive definite C(to)- Because 
of statistical fluctuations, however, the measured correlation matrix C(io) is 
not necessarily positive for tg > 0. This is the reason why one usually considers 
only the to = case in applying the variational method, especially with a large 
number of operators. On the other hand, it is obvious that C(0) is contaminated 
by highly excited states and contains only restricted information on the low 
lying part of the spectrum. Therefore it is desirable to take to > 0. This can be 



achieved in the following way [155|. 
We first diagonalise C{to), 

C{to)ipi ^ Xiipt, Xi>-.->Xn, (E.6) 

and project the correlation matrices to the space of eigenvectors corresponding 
to the M highest eigenvalues, 

C*^(i) = (^,,C(i)^j), z,j-l,...,M. (E.7) 

By choosing the operator space too large we introduce numerical instabilities 
caused by very small (even negative) eigenvalues with large statistical errors 
due to the fact that the chosen operator basis is not sufficiently independent 
on the given MC sample. By choosing M appropriately we can get rid of those 
unstable modes while still keeping all the physical information. In this way we 
render the generalised eigenvalue problem well defined. 

Of course the final result should not depend on the choice of M and one 
has to take care in each case that this is really the case. Our observation is 
that for any acceptable statistics one always finds a plateau in M for which the 
extracted masses are stable under variation of M . 

In a next step we determine the vectors w„, n = 1,... ,Af through the 
generalised eigenvalue equation in the truncated basis: 

C^^(<i)«„ - e-^-^'^-'»^C''ito)vn . (E.8) 

This equation yields the spectrum En- However, the procedure — although it is 



exact for a correlation matrix which has exactly the form of eq. E.4 — is highly 
non-linear, and a small statistical fluctuation can be enhanced by it and cause 
a systematic shift in the energy values obtained, even when the instabilities are 
avoided by the truncation to M < TV. 

In order to avoid this pitfall we use the (approximate) dual vectors Vn ob- 



tained from eq. E.8 to restrict the problem to an even smaller, therefore more 
stable subspace. 

Define the new correlation matrix of size K x K (with K < M) by 

C^^{t)^{v,,C^'{t)vj), i,j = l,... <M. (E.9) 



E.2. Choosing the Right Mass 
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The steps performed until now can be thought of as a preparation for choosing 
the appropriate set of operators, i.e. hnear combinations of original Oa operators 
which effectively eliminate the higher states. The correlation matrix C.^j (t) is 
then considered as a primary, unbiased object. 

The next step is to fit C^j (t) in the range t = t^nin ■ ■ ■ imax using the ansatz 

K 

{t- {i,,E}) = J2 e-'^-'i^mrn, , (E.IO) 

n=l 

where ^ni, En are the free parameters to be fitted. 

Usually we choose K = 1 and 2. For the Af'^ glueball, however, K — 2 and 
3 are chosen since we do not subtract the vacuum contribution (Oq) (Ofj)* from 
the correlators but consider instead the vacuum state together with the glueball 



states in this channel (cf. remarks in Section 2.4.2). 

In the fitting step we use a correlated fit which takes into account the cor- 
relation between (t) and CPy{t'), i.e. using the inverse of the corresponding 
covariance matrix Cov{i, j,t;i\ j' ,t') as a weight in the definition of x^- This 
has the advantage over the uncorrelated that the value of the latter can be 
artificially small if the quantities to be fitted are strongly correlated. Be warned 
however, that (as usually with sophisticated methods) the correlated fit can 



have its own instabilities if the number of data is not sufficiently large [156,157]. 

The result of this fitting process becomes apparent if one compares the initial 
effective mass, e g. |E.2| , to the effective mass obtained using the fitted correlation 



matrix elements C/^ (t) as displayed in Figure E.l for the scalar glueball on the 



coarsest lattice (/? = 3.15) where the plateau is rather short. 



Figure E.2 which displays only the final effective mass plots demonstrates 
that using anisotropic lattices (even at ^ = 2) may greatly improve the quality 
of the results. 



E.2 Choosing the Right Mass 

In the last Section, we described how one obtains estimates of a state's mass 
using sophisticated variational methods, depending on input parameters such 
as the initial input operators, the number M of operators kept after the first 
step or the fitting range t = tmin ■ ■ ■ ^max- The number of degrees of freedom in 
the fit is very large, for typical mass determinations of glueballs it may be of 
0(10000) if one considers all combinations that might turn out to be sensible, 
i.e. that might result in a stable mass plateau and a good fit with small x^. 

The first choice one has to make, is the reduction of parameter sets examined 
to a feasible number, say O(IOO). This can normally be done quite easily, as 
the extracted mass should not depend too much on whether, e.g., one chooses 
M = 13 or M = 14, whether one uses the fitting range t=1...5ori = 1...6 
etc. Figure |E.3| displays the resulting estimates of the full variational method 
with the condition x^/-^df < 2 for the ground state of the A^'^ glueball at /3 = 
3.50 (fixed initial number of operators N = 91). We have included 16 different 
fit ranges as well as 11 different truncations M. In this case it turns out that 
it is unambiguous to pick a mass that is consistent with all the determinations, 
marked by the solid line; the errors of the final result are marked by the dashed 
lines. 
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- a 



O fitted 
□ initial 



Figure E.l: The effective masses, as defined in eq. E.2, before (squares) and after 
(circles) the fitting procedure for the ground state of the scalar glueball on the 
coarse lattice, at ^ = 2, /3 = 3.15. The solid and dashed horizontal lines mark 
the resulting energy and its error, obtained performing a fit on t = 1, . . . ,7. 




Figure E.2: The effective masses, as defined in eq. E.2 after the fitting procedure 
for the separation i? = 1 of the on-axis static quark-antiquark potential at 
(3 = 3. 50. The solid and dashed horizontal lines mark the resulting energy and 
its error, obtained performing a fit on t = 3, . . . ,12 
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Figure E.3: The mass estimates of the scalar gluebah mass on ^ = 2, /3 = 3.50 
configurations. 16 different fit ranges as well as 11 different truncations M are 
included. The finally chosen value and its errors are marked by the solid and 
the dashed lines, respectively. 



Of course, the picture does not look that nice in all the cases investigated. 
Especially for higher masses, where the signal disappears in the noise after just 
a few time slices, one has to take care to study the behaviour of the mass under 
the variation of all the parameters and to pick a region where the mass is stable 
for all sensible choices of parameters. If this is not possible, one might specify 
an average mass estimate with an error that comprises all the mass estimates 
coming from sensible choices of the parameters. However, in this case, the error 
normally gets very large. 



Appendix F 

Representations of 
Two-Glueball States 



In this Appendix we aim at calculating the representation of the cubic group 
Oh to which a pair of glueballs {Tf'~'', T2 ^ ) with opposite momenta p, —p con- 
tributes. These results are used to estimate the minimal energies of glueball pairs 



on the lattice appearing in each representation, see Section 6.5. The method 



has been described by Morningstar and Peardon in and is presented in more 
detail in the following. 



F.l Point Groups Used 



Glue balls at rest transform according to the full cubic group Oh, se e Se ction 
2.4.2| . This group describes the symmetries of a cube (see Figure FT) and 
contains 48 elements in 10 classes: the identity E, 8 three-fold (120°) rotations 
C3, 3 two-fold rotations (180°) C2 around the x,y, z-axes, 6 two-fold rotations 
C2 around the axes a-f and 6 four- fold rotations (90°) C4; the other five classes 
contain the same elements, each preceeded by the spatial inversion /. 




Figure F.l: The rotation axes Oa, Ob, Oc, Od, Oe, Of, Ox, Oy, Oz, Oa, 0(3, O7 
and 05 of the cubic group Oh ■ 
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Momentum breaks full Oh symmetry. Depending on the orientation of 
the spatial momentum, different symmetries are preserved, the so-called lit- 
tle groups, subgroups of Oh , describing these remnant symmetries are listed in 
Table F.l. The group C^v has the symmetries of a regular 4-sided pyramid and 
contains 8 elements; accordingly the groups C2v and Csv have the symmetries 
of 2- and 3- sided pyramids respectively and contain 4 and 6 elements. The lit- 
tle group Cs contains only two operations, namely the identity E and a mirror 
reflection cr; the trivial group Ci includes only the identity 



{nx,ny,nz) little group international notation 



(0,0,0) 


Oh 


— 3— 

m m 


(0,0,n) 


Civ 


Amm 


(0,n,n) 


C2v 


2mm 


(n,n,n) 


C3v 


3m 


(0,m,n) 


Cs 


m 


(m,m,n) 


Cs 


m 


(l,m,n) 


Ci 


1 



Table F.l: The little groups for lattice glueball states having momentum 
p ~ 2n{nx,'ny,nz) / S on a simple cubic spatial lattice with periodic bound- 
ary conditions. Note, that l,m,n ^ 0, I ^ m,m ^ n,l ^ n. The group Cs is 
often denoted as 5*1 . The international notation for the groups is given together 
with the usual Schoenflies notation for convenience. 



F.2 The Method 

In order to work out the representations of Oh a pair of glueballs with net zero 
momentum is contributing to, (1) we subduce the given representations of Oh, 
rf<^ and rf into the little group under consideration (depending on the mo- 
mentum orientation). These representations which in general are reducible are 
then (2) decomposed into direct sums of irreducible representations of the little 
group. Out of these sums, (3) we build pairs of irreducible representations each 
coming from one of the initial glueball representations and build their direct 
product. The direct product is again decomposed into irreducible representa- 
tions of the little group, and finally, (4) the representation is induced back to 
Oh and decomposed into irreducible representations of this group. That this 
is possible follows from the net zero momentum of the glueball pairs. In the 
following, we describe each of these steps in detail. 

^ These groups describe important symmetries occurring in nature, namely the symmetries 
of molecules: e.g., C2v is the symmetry group of molecules like H2O, C^v describes the 
symmetries of NH3. They belong to a large set of groups called molecular point groups. 
A subset of these groups are the 32 crystallographic point groups describing symmetries of 
regular crystal lattices — actually, all the point groups used in this work also belong to this 
subset. 
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First, let us introduce the most important tool that will be used, namely 
characters. The sum of all eigenvalues of a matrix representation D{Ga) of 
the group element Ga, called the character x(Go), is invariant under similarity 
transformations D'(Ga) — AD{Ga)A~^ . Furthermore, all the elements of a 
class Cp of a group have the same character Xpi because 

x{Ga) ^ ^A»(G.) = ^A,(G„G,G-i) 

i i 

= J2D^J{G,n)DJk{Gt)Dk^{G;^') 
ijk 

= Djk{Gb)Dkj{G:^Gra) 

jk 

- Y,D,,{G,)^X{G,). (F.l) 

The characte r s of the cla s ses of all the point groups used here are given in Tables 
fJ, HJ, |F.11| , |F.l4 |F.17[ 



These characters are very useful if one is reducing a representation, i.e. de- 
composing a given (reducible) representation into its irreducible contents 

^ = 0m«T("), (F.2) 

a 

where a runs only over the inequivalent irreducible representations (irreps) T^"^ 
and rria gives the number that this irrep occurs in the reduction. The character 
of the reducible representation D is related to the irreducible characters Xp'^ by 
a similar equation: 

OL 

The vital task of deducing the coefficients rria given the representation D is 
facilitated a lot by the use of characters. Using the orthogonality between 
different irreducible representations, one may show 

- X! '^P^f^*Xp = - X! X! CpXp''^*Xp"^ = ~ X! '^o^a^ap = (F.4) 

a p a a p a a 

where g is the order of the group, Cp denotes the number of elements in a given 
class, p runs over the number of classes whereas a numbers the irreps. The 
coefficient mp is thus just the scalar product of the character x of the reducible 
representation with x^^\ the character of the irrep|^ 

Having prepared the most important tools, we may now start to apply 
the method step by step. As a random example, we study the glueball pair 
(Tj++, Aj^") with momenta p = ±2tt/S{1, 1, 1). 

F.2.1 The Subduction into the Little Group 

First, we subduce the characters of Oh into the little group under consideration. 
This amounts to simply picking the right columns out of Table corresponding 



^Note that all the numbers must come out as positive integer numbers or zero — this 
is a considerable safeguard against errors! 
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to the remnant symmetry operations in the little group under consideration. 



These correspondences are listed in Table F.4. For our example, with the little 



group Csy, we pick the columns denoted by E, 8C3 and 6/C2 corresponding to 
the classes E, 2C3 and Scti, of the little group. 



F.2.2 The Decomposition into Irreps 

The subduced representations, denoted by Oh i Gi, whe re G i denotes the lit- 
tle group, are in general reducible in G;. Applying eq. FA to the subduced 



characters and the respective characters of the little group, we obtain the de- 
composition of all the init i al representation s into the irreps of the little groups, 
listed in Tables |K7[ |F.1C| , |F.13| , | .16| , |F.18[ For our example groups, we build 
the scalar products of the subduced characters of Oh prepared in the last section 



and the characters of C^y listed in Table F.ll obtaining Table F.13 



F.2.3 Building Glueball Pairs 

Up to now, the steps performed referred to single glueballs having momentum p. 
In order to obtain physical objects transforming the same way as our single zero- 
momentum glueballs, we will now combine two glueballs of opposite momentum. 
This can be obviously done in the little group corresponding to the momentum. 
A pair of glueballs transforms according to the direct product of two irreps. 
On the level of characters this amounts to multiplying the characters of the 
irreps corresponding to the same group elements (or class). Having in general 
reducible representations of the initial representations decomposed into irreps, 
there may be several direct products to be taken into account. In our example, 
T^^ I Gsy = A2 & E and i (^Sv = ^2- This leaves us with the direct 

products A2 (8) A2 and E A2. Calculating the characters, one notices easily 
that A2(^A2 = Ai and E(^A2 = E. That the direct products are irreducible is 
true for most of the cases; for the other cases, the reduction can be performed 



exactly the same way as presented in Section F.2 and applied in Section F.2.2 



For convenience, the multiplication tables (together with the decomposition into 
irreps) are given as Tables R3, f\6, |F.9| , F.12| , F.15 (the table for the group Ci 



is completely trivial and thus omitted). 

In the case of both glueballs having momentum 0, i.e. transforming according 
to Oh, the whole method amounts to simply performing the direct products and 
reading off the resulting irreps of Oh in Table F.3. 



F.2. 4 The Induction back to the Cubic Group 

For most of the representations of the glueball pairs in the little groups, the 
induction (getting back onto the level of the full cubic group Oh) is trivial, due 
to the Frobenius reciprocity theorem that states that the number of times a 
representation D'^'^\G) occurs in D'^^'^ ] G, the induced representation, is equal 
to the number of occurrences of the irrep L»(^)(G,) in i Gs, i.e. if the 

representation of the glueball pair in the little group is again an irrep of the 
little group, we can read off the conte nts of the i nduced rep r esent ation in terms 



of the irreps of Oh again from Tables F.10| , |F.13| , |F.16| , |F.18 . 

For the two cases {E®E in C^y and E®E in C^y) where the direct product is 
not an irrep of the little group, we just give the results: {E®E) ] Oh = E~^(BE~ 
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for Civ , whereas E ® E contributes to all representations of Oh for C^^ . The 
way of performing these inductions is employing the left coset decomposition of 
Oh in terms of the little group, see e.g. [158, 159|. 

In our example, the irreps of the little group describing the glueball pair are 
Ai and E. Using the Frobenius theorem with Table F.13 , we notice easily that 
the pair may contribute to the irreps A'i[ , E~^ , , T2 , ,E~, , Tj" of the 
cubic group Oh- 



F.2.5 A Note on Parity and Charge Conjugation 

In the previous sections, the parity quantum number P = ± we have always 
taken into account, whereas the charge conjugation C has been omitted. The 
reason for this is simple: once the momentum of a single glueball is non-zero, 
it is no longer an eigenstate of parity, i.e. that is why we have to consider the 
decomposition into irreps of the little groups for P = + and P — ~ states 
separately. Charge conjugation is completely different: the quantum number 
C2 of the glueball pair is always just the product C • C" of the quantum numbers 
of the initial states, therefore the results for C = C" = + can be trivially 
generalised to the other three cases. 

There is another important issue concerning parity, namely if the two glue- 
balls are indistinguishable (having the same F'^'-^), the final two-glueball state 
has to have positive parity due to the Bose-symmetry. Hypothetic P = — glue- 
ball pairs consisting of two states having the same quantum numbers thus have 
to be excluded. 

For our example, all this means that because of the two representations, 
T^^ and A\ , being different, Bose symmetry does not exclude any states; 
furthermore that the resulting states have charge conjugation quantum number 
C — ~. The two glueballs of our example, having momenta p parallel to the 
space diagonal thus contribute to the representations A'l , E^ , , , 



F.3 Results 

We calculate the energies of all possible glueball pairs on all our lattices for 
momenta up to |p| = 3 • 2-K/Sas. The resulting lowest-lying masses of glueball 
pairs, for all the representations where states can be measured at a given (3 are 
given in Tables p.9| - G.ll| . The spectra are displayed in Figures |F.2| - |F.4| . 
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Table F.2: The character table for the cubic point group Oh- 
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At 




E+ 
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Table F.3: The direct products of the irreducible representations (irreps) of the 
cubic point group O, decomposed back into direct sums of the irreps. For the 
(trivial) extension to the full cubic group Oh, sec text. Due to the obvious 
symmetry, only the upper diagonal part of the table is given. 
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Table F.4: The correspondence of the classes of the cubic point group Oh to the 
classes of the little groups €4^, C2V, Csvi Os,Ci. 
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Table F.5: The character table for the point group 
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Tabic F.6: The direct products of the irreducible representations (irreps) of the 
point group Civ, decomposed back into direct sums of the irreps. Due to the 
obvious symmetry, only the upper diagonal part of the table is given. 
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Table F.7: The decomposition of the subduced representations Oh i Civ in 
terms of the irreducible representations of the little group. 
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Table F.8: The character table for the point group 
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Table F.9: The direct products of the irreducible representations (irreps) of the 
point group C2V, decomposed back into direct sums of the irreps. Due to the 
obvious symmetry, only the upper diagonal part of the table is given. 
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Table F.IO: The decomposition of the subduced representations Oh J, C2v in 
terms of the irreducible representations of the little group. 
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Table F.ll: The character table for the point group C^v 
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Tabic F.12: The direct products of the irreducible representations (irrcps) of 
the point group C^y^ decomposed back into direct sums of the irreps. Due to 
the obvious symmetry, only the upper diagonal part of the table is given. 
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Table F.13: The decomposition of the subduced representations Oh i C^^ in 
terms of the irreducible representations of the little group. 
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Table F.14: The character table for the point group 



Cs Ay A2 

A\ A\ A2 
A2 A2 Ai 



Table F.15: The direct products of the irreducible representations (irreps) of 
the point group Cs, decomposed back into direct sums of the irreps. 
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Table F.16: The decomposition of the subduced representations Oh i Cg in 
terms of the irreducible representations of the little group. 
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Table F.17: The character table for the point group Ci. 
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Table F.18: The decomposition of the subduced representations Oh i C\ in 
terms of the irreducible representations of the little group. 
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Figure F.2: The two-glueball spectrum on the /? = 3.15 lattice. 
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Figure F.3: The two-glueball spectrum on the /3 = 3.30 lattice. 
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Figure F.4: The two-glueball spectrum on the (3 = 3.50 lattice. 
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p2 


M 


fit range 




XV^DF 


3.0 





2 


1-5 


1.368(30) 


0.19 




1 


3 


1-6 


1.372(17) 


0.49 
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3 


1-6 


1.451(17) 


0.28 
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3 


1-8 


1.570(21) 


0.46 
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3 


1-6 


1.610(19) 


0.34 
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1-3 


1.707(29) 


0.05 
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2-5 


0.688(16) 
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1-5 


0.806(5) 


0.06 
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0.902(6) 


0.59 
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1-5 


1.059(8) 


2.35 
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1-5 


1.151(7) 


0.16 
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1-5 


1.310(15) 


0.71 
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1-5 


1.356(23) 


1.58 
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1-3 


1.470(17) 


0.42 
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3 


1-3 


1.621(23) 


0.01 




18 


3 


1-3 


1.762(47) 


0.56 


3.15'' 
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3 


1-5 


1.418(23) 


0.05 




2 


3 


1-5 


1.541(37) 


0.53 
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3 


1-4 


1.751(67) 


1.20 




5 


3 


1-5 


1.922(104) 


1.08 



Table G.l: Collection of results of the torelon measurements using the f = 2 
perfect action. For each /3- value and momentum = +Py "we list the number 
of operators M kept after the first truncation in the variational method, the 
plateau region on which the fit of the correlators to the form Z{p^) exp{—tE{p'^)) 
is performed (fit range), as well as the extracted ener gy t ogether with the 



per degree of freedom (x /Ndf). {Continued in Table GS) 
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p2 


M 


fit range 




XV^DF 


3.3 





3 


1-10 


0.318(6) 


1.39 




1 


3 


1-10 


0.417(3) 


1.20 
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3 


1-5 


0.493(3) 


0.13 
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4 
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0.633(4) 


1.26 
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1-8 
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1-6 
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0.42 
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5 


1-10 


0.921(4) 


0.77 
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5 


1-10 


1.036(6) 


0.30 
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4 


1-6 
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5 


3-8 


0.305(6) 
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5 
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4 
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0.538(4) 
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0.586(3) 
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5 
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0.721(4) 


0.09 
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5 
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0.11 
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4 


1-8 


0.795(4) 
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5 
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0.901(4) 


1.26 
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1.010(16) 
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18 


3 


1-5 


1.215(12) 


0.17 



Table G.2: Collection of results of the torelon measurements using the = 2 
perfect action. For each /3- value and momentum = p^+Py we list the number 
of operators M kept after the first truncation in the variational method, the 
plateau region on which the fit of the correlators to the form Zij?^ exp(— tiJ (p^)) 
is performed (fit range), as well as the extracted energy together with the 
per degree of freedom (x^/-^df). {Continuation from Table 
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atV{f) 




1,0,0) 


1 
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2-7 


0.5053(3) 


0.70 


1,1,0) 


1.414 
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3-5 


0.7322(7) 


0.01 


1,1,1) 


1.732 


3 


3-8 


0.8902(14) 
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2,0,0) 


2 


3 


2-8 


0.9407(10) 


0.89 


2,1,0) 


2.236 


3 


2-8 


1.0642(11) 


0.36 


2,1,1) 


2.449 
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2-5 


1.1676(13) 


0.27 


2,2,0) 


2.828 


3 


2-6 


1.3082(18) 


0.17 


3,0,0) 


3 


3 


2-6 


1.3392(25) 


0.16 


2,2,1) 


3 


3 


2-5 


1.3867(20) 


0.26 


2,2,2) 


3.464 


3 


2-6 


1.5675(55) 


0.84 


4,0,0) 


4 


2 


1-5 


1.7429(24) 


0.14 


4,2,0) 


4.472 


3 


2-6 


1.9328(93) 


1.47 


4,2.2) 


4.<S!)9 


;] 


1 4 


2 1 lOT^.'W) 


0.15 



Table G.3: Collection of results of the off-axis qq measurements using the £, =2 
perfect action at (3 = 3.00. For each separation vector we list the length of the 
vector in spatial lattice units, the number of operators M kept after the first 
truncation in the variational method, the plateau region on which the fit of the 
correlators to the form Z{f) cxp(—tE{f)) is performed (fit range), as well as the 
extracted energy together with the per degree of freedom (x^/A^df)- 
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r 


\r\ 


M 


fit range 


atV{r) 




(1,0,0) 


1 


3 


3 10 


0.3006(1) 


0.20 


(1,1,0) 


1.414 


3 


3 10 


0.3814(2) 


0.21 


(1,1,1) 


1.732 


3 


3-10 


0.4250(3) 


0.65 


(2,0,0) 


2 


3 


6-10 


0.4469(6) 


0.22 


(2,1,0) 


2.236 


3 


4-10 


0.4752(4) 


0.41 


(2,1,1) 


2.449 


3 


4-10 


0.4971(5) 


0.83 


(2,2,0) 


2.828 


3 


4-10 


0.5309(7) 


0.59 


(3,0,0) 


3 


3 


4-10 


0.5445(7) 


1.71 


(2,2,1) 


3 


3 


5-10 


0.5456(8) 


0.62 


(2,2,2) 


3.464 


3 


4-10 


0.5852(9) 


1.53 


(4,0,0) 


4 


3 


4-10 


0.6269(1) 


0.51 


(3,3,0) 


4.243 


3 


5-10 


0.6455(16) 


0.94 


(4,2,0) 


4.472 


3 


5-10 


0.6630(14) 


0.86 


(4,2,2) 


4.899 


3 


5-10 


0.6964(18) 


0.34 


(5,0,0) 


5 


2 


5-10 


0.7044(22) 


1.24 


(3,3,3) 


5.196 


3 


4-10 


0.7215(17) 


1.31 


(4,4,0) 


5.657 


2 


6-10 


0.7509(41) 


1.23 


(4,4,2) 


6 


3 


5-10 


0.7783(25) 


0.77 



Tabic G.4: Collection of results of the off-axis qq measurements using the ^ = 2 
perfect action at /3 = 3.30. For each separation vector we list the length of the 
vector in spatial lattice units, the number of operators M kept after the first 
truncation in the variational method, the plateau region on which the fit of the 
correlators to the form Z{f) exp{—tE{-rj) is performed (fit range), as well as the 
extracted energy together with the per degree of freedom (x^/Nuf)- 
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M 


lit range 


<HV{r) 




3.15 


1 


5 


3-10 


0.3703(2) 


0.456 




2 


5 


3-7 


0.6110(4) 


0.735 




3 


5 


3-10 


0.8019(8) 


1.209 




4 


5 


3-7 


0.9775(16) 


1.238 




5 


•l 


4 !) 


1.12!)(<Sj 


O.-ll!) 


3.50 


1 


5 


3-12 


0.25910(5) 


1.163 




2 


5 


3-10 


0.36490(12) 


0.769 




3 


5 


7-12 


0.42168(38) 


0.664 




4 


5 


5-11 


0.46587(51) 


0.205 




5 


5 


5-12 


0.50388(76) 


1.139 




(i 


•l 


.") 11 


()..')4().")1 (!)()) 


0.8.");-! 



Table G.5: Collection of results of the on- axis qq measurements using the ^ = 2 
perfect action at /3 = 3.15, 3.50. For each separation along the axes we list 
the number of operators M kept after the first truncation in the variational 
method, the plateau region on which the fit of the correlators to the form 
Z{r) exp{—tE{r)) is performed (fit range), as well as the extracted energy to- 
gether with the per degree of freedom (x^/Nuf)- 
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Channel 


iV 


to/h 


M 


fit range 


XV^DF 


energies 


At^ 


87 


1/2 


12 


1-4 


0.69 


0.645(10) 


A++* 


59 


1/2 


10 


1-3 


1.23 


1.365(104) 


E++ 


47 


1/2 


7 


1-3 


0.10 


1.405(29) 


E++* 


47 


1/2 


5 


1-3 


0.17 


1.985(140) 


T++ 


47 


1/2 


10 


1-4 


0.10 


1.416(32) 


T++* 


22 


1/2 


6 


1-3 


1.20 


1.990(131) 


At^ 


21 


1/2 


2 


1-3 


0.03 


1.675(147) 


T++ 


36 


1/2 


4 


1-3 


0.34 


1.881(131) 


T++* 


36 


1/2 


6 


1-3 


0.31 


2.170(182) 




25 


1/2 


3 


1-3 


0.001 


1.588(114) 


E-+ 


25 


1/2 


2 


1-3 


0.01 


1.849(81) 


E-+* 


25 


1/2 


2 


1-3 


2.19 


2.084(268) 


T,-+ 


64 


1/2 


4 


1-3 


0.02 


1.850(83) 




61 


1/2 


4 


1-3 


0.24 


1.798(59) 


r+- 


18 


1/2 


4 


1-3 


0.41 


1.935(177) 



Table G.6: Results from fits to the ^ = 2, /3 = 3.15 glueball correlators on the 
8^ X 16 lattice in units of the temporal lattice spacing at: ta/ti are used in the 
generalised eigenvalue problem, TV is the number of initial operators measured 
and M denotes the number of operators kept after the truncation in C(to). 
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Channel 


to/ti 


M 


fit range 


XV^DF 


energies 


At+ 


1/2 


9 


2-7 


0.90 


0.590(17) 


A++* 


1/2 


11 


1-3 


1.57 


1.133(53) 


E++ 


1/2 


11 


1 - 5 


0.97 


0.983(17) 


E++* 


1/2 


8 


1 - 3 


0.90 


1.453(68) 


T++ 


1/2 


12 


1-4 


1.09 


0.965(15) 




1/2 


10 


1-4 


1.68 


1.386(51) 




1/2 


4 


1 - 3 


0.01 


1.511(94) 


T++ 


1/2 


8 


1-5 


0.84 


1.492(62) 


AI^ 


1/2 


6 


1-4 


4.38 


1.017(39) 


E-+ 


1/2 


3 


1-4 


0.74 


1.365(36) 




1/2 


4 


1-5 


1.16 


1.357(25) 




1/2 


8 


1-4 


0.95 


1.276(33) 



Table G.7: Results from fits to the ^ = 2, /3 = 3.3 glueball correlators on the 
10^ X 20 lattice in units of the temporal lattice spacing a*: to/ti are used in the 
generalised eigenvalue problem, M denotes the number of operators kept after 
the truncation in C{to). 
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4-6 


0.30 


0.720(158) 


E++ 
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23 


2-5 


0.16 


0.675(25) 


E++* 
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14 
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0.70 
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r++ 


53 
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19 


2-5 


0.31 
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53 


1/2 


22 
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15 
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6 


1-3 


0.05 
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15 


1/2 


6 
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T++ 
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0.05 
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A^+* 


15 


1/2 


9 


1-3 


0.61 
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E-+ 


21 


1/2 


14 


1 - 3 


0.27 


0.942(24) 


E-+* 


21 


1/2 


11 


1-3 


1.20 


1.406(80) 


rp- + 


75 


1/2 


9 


2-5 


0.32 


0.952(38) 




75 


1/2 


6 


1-3 


1.70 


1.675(59) 


r+- 


53 


1/2 


14 


1-3 


0.69 


1.013(17) 




11 


1/2 


5 


1 - 3 


0.07 


1.224(64) 


At-* 


11 


1/2 


5 


1 - 3 


0.79 


1.726(147) 


rp+- 


22 


1/2 


9 


2-4 


0.33 


1.142(121) 


rn+-* 
^2 


12 


1/2 


9 


1 - 3 


1.78 


1.663(139) 


E+- 


50 


0/1 


50 


1-3 


1.01 


1.379(34) 


Tf- 


80 


0/1 


80 


2-4 


0.57 


1.200(67) 




80 


0/1 


80 


2-5 


0.20 


1.205(81) 


Ar 


35 


0/1 


35 


1 - 3 


0.16 


1.352(37) 


Ar 


30 


0/1 


30 


1 -4 


0.14 


1.703(99) 



Table G.8: Results from fits to the ^ = 2, /? = 3.5 glueball correlators on the 
12'^ X 24 lattice in units of the temporal lattice spacing at: to/ti are used in the 
generalised eigenvalue problem, TV is the number of initial operators measured 
and M denotes the number of operators kept after the truncation in C{to)- 
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Channel 


Glueballs 




To ■ m2G 


Ar 


{At\ At^) 





5.16 


E++ 


{At\ At^) 


1 


6.07 


T++ 


{At\ 


2 


6.86 




{A++, E++) 


1 


8.89 


T++ 


{At\ 


5 


8.82 




A^+) 





8.94 


E-+ 


iAt+, T^) 


1 


8.93 




{A++, E++) 


1 


8.89 


T+- 


(^++, T+-) 





9.78 


T+- 


{At+, T+-) 





10.33 



Table G.9: The energies of the lowest-lying two-glueball states with zero mo- 
mentum for all the symmetry channels, given in terms of ro, using the masses 
determined in the /3 = 3.15 simulation. The square of the momenta of the 
single glueballs building the pair is given in units of (27r/8as)^ = (1.60/ro)^. 



Channel 


Glueballs 


p' 


ro ■ m2G 


Ar 


{At\ 


AV) 





7.10 


E++ 


(At^, 


A^) 


1 


8.13 


T++ 


iAt\ 


AV) 


2 


9.04 


At^ 




E++) 


1 


10.32 


T++ 


{Af, 


rr++\ 

^2 ) 


1 


10.21 


AI^ 









9.68 


E-+ 


{At\ 


^2 ) 


1 


10.21 




(Ar, 


^2 ) 


1 


10.21 


T+- 









11.24 



Table G.IO: The energies of the lowest-lying two-glueball states with zero mo- 
mentum for all the symmetry channels, given in terms of tq, using the masses 
determined in the (3 = 3.30 simulation. The square of the momenta of the 
single glueballs building the pair is given in units of (27r/10as)^ = (1.98/ro)^. 
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Channel 


Glueballs 




ro • miG 


A++ 







7.30 


E++ 




1 


8.93 


r++ 







9.78 




£;++) 


1 


11.06 


T++ 




1 


11.11 









10.44 


E-+ 




1 


11.11 


T,^ 


£;++) 


1 


11.06 


T+- 







12.77 




At-) 





14.67 


rp+- 







13.94 


E+- 


K+, r+-) 


2 


14.97 


Tf- 


r+-) 


1 


13.94 




T+-) 


1 


13.94 




T+-) 


2 


14.97 


A~ 




1 


13.94 



Table G.ll: The energies of the lowest-lying two-glueball states with zero mo- 
mentum for all the symmetry channels, given in terms of ro, using the masses 
determined in the /? = 3.50 simulation. The square of the momenta of the 
single glueballs building the pair is given in units of (27r/12as)^ = (2.57/ro)^. 
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Ratio 


P = 3.15 


P = 3.30 


f3 = 3.50 


m{A++*)/m{A++) 


2.12(19) 


1.92(15) 


1.78(45) 


m{E++)/m{A++) 


2.18(8) 


1.67(8) 


1.67(12) 


m(S++*)/m(A++) 


3.08(26) 


2.46(19) 


2.92(23) 


miT++)/m{At+) 


2.20(8) 


1.64(7) 


1.68(8) 


m{T++*)/m{A++) 


3.09(25) 


2.35(15) 


2.79(15) 


m{A++)/m{A++) 


2.60(27) 


2.56(23) 


2.99(20) 


miA++*)/m{A++) 






3.87(45) 


m(r++)/m(A++) 


2.92(25) 


2.53(18) 


3.07(16) 


miA^+)/m{A++) 


2.46(21) 


1.72(12) 


1.86(12) 


m{A-+*)/m{A++) 






2.81(37) 


m{E-+)/m{A++) 


2.87(17) 


2.31(13) 


2.33(13) 


m{E-+*)/m{At+) 


3.23(47) 




3.47(31) 


miT,+)/m{A++) 


2.87(17) 


2.30(11) 


2.35(17) 


m{T^-+*)/m{A++) 






4.14(28) 


m{T+-)/m,{A++) 


2.79(13) 


2.16(12) 


2.50(12) 


m{A+-)/m,{A++) 






3.02(26) 


m{T+-)/m{A++) 






2.82(39) 



Table G.12: Glueball mass ratios with the scalar glueball mass, mG/Tn{Af~^). 
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Ratio 


/? = 3.15 


13 = 3.30 


P = 3.50 


m(A++)/m(r++) 


0.46(2) 


0.61(3) 


0.59(3) 


miAt+*)/m{T++) 


0.96(10) 


1.17(7) 


1.06(25) 


m{E++)/m{T++) 


0.99(4) 


1.02(3) 


0.99(5) 


m{E++*)/m{T++) 


1.40(13) 


1.51(9) 


1.74(11) 


m{T++*)/m{T++) 


1.41(19) 


1.44(8) 


1.66(6) 


m{At+)/m{T++) 


1.18(13) 


1.57(12) 


1.78(9) 


m{A++'*)/m{T++) 






2.30(23) 


m{T++)/m{T++) 


1.33(12) 


1.55(9) 


1.83(7) 


m{A^+)/m{T++) 


1.12(11) 


1.05(6) 


1.11(5) 


miA^+*)/m{T++) 






1.67(19) 


m{E'+)/m{T++) 


1.31(9) 


1.41(6) 


1.38(6) 


m{E-+*)/m{T++) 


1.47(22) 




2.06(15) 


m{T^+)/m{T++) 


1.31(9) 


1.41(5) 


1.40(8) 


m{T^+*)/m{T++) 






2.46(13) 


m{T+-)/m{T++) 


1.27(7) 


1.32(5) 


1.49(5) 


m(At')lm{T++) 






1.80(13) 


m{T+-)/m{T++) 






1.68(21) 



Table G.13: Glueball mass ratios with the tensor glueball mass, mG/m{T2"^). 
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Channel 


J 


terms in the fit 


XV^DF 




At+ 





C, C6 


0.99 


3.64(10) 


A++* 





C, C6 


1.06 


6.90(38) 


E++ 


2 


C, C2 


1.06 


6.16(24) 


E++* 


2 


C, C2 


3.85 


10.59(63) 


T++ 


2 


C, C2 


1.98 


6.14(18) 




3 


C, C2 


1.65 


11.65(57) 







c 


2.50 


6.46(18) 


E-+ 


2 


C, C2 


1.06 


8.73(34) 


m-+ 
-'-2 


2 


C, C2 


1.07 


8.77(39) 



Table G.14: Results of the continuum extrapolations of selected glueball repre- 
sentations in terms of the hadronic scale tq. The continuum spin assignment 
J, the terms in the fit, constant (c), (a^/ro)^ (C2) and (fls/ro)^ {c%) and the 
goodness of the fit, x^/-^df, are also given. 





j 


terms in the fit 




mppc /mj,++ 


romc 


E++ 


2 


c 


1.20 


1.01(2) 


6.20(30) 


E++* 


2 


C, C2 


2.20 


1.74(12) 


10.68(105) 




2 


C, C2 


3.01 


1.70(8) 


10.44(80) 




3 


C, C2 


1.11 


1.90(11) 


11.66(102) 


T++ 


3 


C, C2 


2.60 


1.90(8) 


11.66(83) 


AI^ 





c 


1.39 


1.09(4) 


6.69(44) 


T2-+ 


2 


C, C6 


1.00 


1.41(4) 


8.66(50) 


r+- 


1 


C, C2 


2.29 


1.54(7) 


9.45(71) 



Table G.15: Results of the continuum extrapolations of selected glueball ratios, 

mQ/mrp++. The continuum spin assignment J, the terms in the fit, constant 

(c), {as/roY (C2) and (ag/ro)^ (ce) and the goodness of the fit, x^/A^df, are 
also given. The last column lists the masses converted to units of Tq ^ using the 
continuum rc;siilt for the mass of the tensor glueball T2"'', values in bold face 
will be used further on. 
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J p- 


AI 


lit range 


aiE(p-} 




3.0 1 


5 


3-10 


0.317(5) 


1.34 


2 


3 


4-10 


0.375(8) 


0.31 


4 


3 


3-11 


0.487(7) 


0.65 


5 


5 


4-11 


0.520(11) 


0.35 


8 


5 


3-7 


0.663(17) 


0.09 


9 


3 


4-11 


0.754(37) 


0.87 


10 


5 


4 11 


().()7<S(;-',2) 


1.56 


Table G.16: ' 


Collection of 


results of the torelon 


measurement at [3 = 


3.0, using 



the £, = 4: perfect action. For each momentum = + Py we list the number 
of operators M kept after the first truncation in the variational method, the 
plateau region on which the fit of the correlators to the form Z{p^) exp{—tE{p'^)) 
is performed (fit range), as well as the extracted energy together with the 
per degree of freedom (x^/A^df)- 
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